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ABSTRACT 
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<^> We review the current status of nonperturbative studies of gauge field theory using the Dyson-Schwinger 

equation formalism and its application to hadronic physics. We begin with an introduction to the 

formalism and a discussion of renormalisation in this approach. We then review the current status of 

studies of Abelian gauge theories [e.g., strong coupling quantum electrodynamics] before turning our 

^\ attention to the non- Abelian gauge theory of the strong interaction, quantum chromodynamics. We 

| discuss confinement, dynamical chiral symmetry breaking and the application and contribution of these 

^ ■ techniques to our understanding of the strong interactions. 
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1 Introduction 

As computer technology continues to improve, lattice gauge theory [LGT] will become an increasingly 
useful means of studying hadronic physics through investigations of discretised quantum chromody- 
namics [QCD]. [For a recent review of LGT with numerous references see Rothe (1992).] However, it 
is equally important to develop other complementary nonperturbative methods based on continuum 
descriptions. In particular, with the advent of new accelerators such as CEBAF and RHIC, there is a 
need for the development of approximation techniques and models which bridge the gap between short- 
distance, perturbative QCD and the extensive amount of low- and intermediate-energy phenomenology 
in a single covariant framework. Cross-fertilisation between LGT studies and continuum techniques 
provides a particularly useful means of developing a detailed understanding of nonperturbative QCD. 

One such continuum approach is based on the infinite tower of Dyson-Schwinger equations [DSEs] . The 
DSEs are coupled integral equations which relate the Green's functions of a field theory to each other. 
Solving these equations provides a solution of the theory; a field theory being completely defined when 
all of its n-point Green's functions are known. The DSEs include, for example, the Bethe-Salpeter 
equation [BSE] which is needed for the description of relativistic two-body scattering and bound states. 
Quantitative studies of a field theory must be based on one [or more] systematic approximation schemes: 
for example, perturbation theory for weak coupling; lattice studies using increasing numbers of sites 
and /^-values in order to represent a larger volume of the spacetime-continuum; large N expansions for 
SU(N) gauge theories; h expansions for semiclassical approximations; etc. For studies based on DSEs it 
is unavoidable that the infinite tower of coupled equations must be truncated at some point. This means 
that the tower of equations must be limited to some n, where n is the maximum number of legs on any 
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Green's function included in the self-consistent solution of the equations. Because of the complex issues 
involved and the computational effort required, current efforts are necessarily modest; for example, we 
know of no attempts as yet to study coupled DSEs for n > 3. The penalty incurred by truncation is the 
need to employ an Ansatz for the omitted function(s). However, much can be achieved by exploiting 
the need to maintain certain properties of the theory, including the various global and local symmetries, 
multiplicative renormalisability, analyticity, known perturbative behaviour in the weak coupling limit, 
etc. These can provide stringent constraints on the Ansatze. In addition, there is the hope that future 
lattice gauge theory simulations will be able to provide additional insight into the form of these higher 
n-point functions. 

In Sec. [2] we introduce the DSE formalism, beginning with a discussion of Abelian gauge theories 
[quantum electrodynamics in three- and four- dimensions, QED3 and QED4 ] and then extend this 
discussion to QCD. In Sees. |3| and f| we review in some detail the current status of DSE-based studies 
of QED3 and QED4 with particular emphasis on gauge covariance and multiplicative renormalisability. 
Such studies are extremely useful as a guide to the more complicated case of QCD. We devote Sec. |] to 
a study of the gauge boson sector of QCD and focus on studies of the infrared behaviour of the gluon 
propagator, since this is thought to be crucial to confinement in QCD. In Sec. |6] we examine the quark 
sector of QCD, discussing the crucial issues of dynamical chiral symmetry breaking [DCSB] and quark 
confinement. In Sec. [7| we review applications of the ideas discussed in the preceding sections to studies 
of hadronic structure. Finally, in Sec. |8], we summarise and discuss possible future extensions of these 
studies. 



2 Dyson-Schwinger Equation Formalism 

It has been known for quite some time that, from the field equations of a quantum field theory, one 
can derive a system of coupled integral equations relating the Green's functions for the theory to each 
other (Dyson, 1949; Schwinger, 1951). This infinite tower of equations is sometimes referred to as the 
complex of Dyson-Schwinger equations. An introduction to the formal details of DSEs, including their 
derivation etc., can be found in a number of text books; Bjorken and Drell (1965, pp. 283-376) and 
Itzykson and Zuber (1980, pp. 475-481) are two examples. For completeness we will summarise those 
aspects of the formalism which are relevant to the gauge field theories of quantum electrodynamics 
[QED] and quantum chromodynamics [QCD]. These are the areas where the DSE approach to the 
solution of a quantum field theory [QFT] has been applied most widely. Our exposition employs the 
functional integral formulation of these field theories; useful introductions to which can be found in 
Itzykson and Zuber (1980, pp. 425-474) and Rivers (1987). The separate but related question of 
whether the functional integral formulation of these field theories can be rigorously defined is as yet 
unresolved. A discussion of this problem can be found, for example, in Seiler (1982). The Euclidean- 
space, discretised, lattice theory is well-defined, but there the question is simply deferred and becomes 
a question of rigorously establishing the existence of the continuum limit. This limit is a critical point 
of the lattice theory, since, in units of the lattice spacing, physical correlation lengths must diverge. 



2.1 Quantum Electrodynamics [QED] 

In the following we will discuss quantum electrodynamics in D = 2 and D = 3 space dimensions, referred 
to as QED 3 and QED 4 respectively. We denote the number of space-time dimensions as d = D + 1 and 
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write the action for QEDd in Minkowski metric, with signature goo = 1 and gu = — 1 for i — 1 . . . D: 



d d x 



N 



/=i 4 



(2.1) 



[We have set % = 1 = c] In Eq. Q2.1p the superscript / is a flavour label; / = 1 . . . Nf for a theory with 
Nf distinct types or flavours of electrically active fermions, each represented by the field ijj^{x)\ ml 
and e{ are, respectively, the bare mass and charge of each of these fermions; A^{x) is the gauge-boson 
[photon] field; and 

(2.2) 



d v A, 



is the Abelian gauge-boson field strength tensor. We use standard notation and conventions, where 

etc. Note that for an electron we would have 



27°7 1 7 2 7 3 , 



4 = A^r = g^A^, {Y,r} = 2<r, 75 

for the [physical] charge e* = — e, where by definition e = |e| is the magnitude of the electron charge. 



The quantum field theory associated with Eq. ( |2.1| ) is defined by the generating functional 



dfi(ijj, if), A) exp iS[ip, ip, Afj] +i d d x 



J2(y f r} f + rj f iP f )+A^ 
f 



(2.3) 



where 77*, rj^ and J M are, respectively, source fields for the fermion, antifermion and gauge boson and 
where we have defined 

dfi@, 1>,A) = H Vi> f V^ J] VA, . (2.4) 
/ n 

[A constant normalisation factor is understood which ensures that £[0,0,0] = 1.] This generating 
functional is the field theoretic analogue of the partition function of a statistical mechanical system and 
serves the same purpose; i.e., all the physical quantities of the theory can be obtained from Z\fj,rj, J M ]. 



Operationally, the integral in Eq. (|2.3| ) represents "an integral over all possible values of each of the 
fields at all spacetime points". Detailed discussions of the functional integral can be found in Rivers 
(1987), Pascual and Tarrach (1984, pp. 198-226) and Seiler (1982). For our purposes, however, we may 
proceed upon noting that the fermion fields {ip, ijj} and their sources {r], 77} are elements of a Grassmann 
algebra with involution which entails that all of these fields anticommute with each other [this ensures 
that Fermi-statistics are obeyed by the fermion fields]; and that A^ and its source J M are c-number 
fields. 



In order to complete the operational definition of QEDd we note that the action in Eq. Q2.1p is invariant 
under the local Abelian gauge transformation 

i{>(x) -> ty(z) = e- ieoX{x) ^j{x) , tf(x) -> x ^(x) = e ieoX{x) ^{x) , 
A^x) -> % t (x) = A^x) - d^\(x) (2.5) 

where A (2) is an arbitrary scalar function and, for now, we have suppressed the flavor indices. This 
entails that the generating functional as it has been defined so far is meaningless [even assuming that 
our above caveats have been properly taken into account]. This is because gauge invariance ensures that 
for any field configuration, {ip(x), ip(x), A^x)}, there are [uncountably] many related configurations, 
{\p(x), x ip(x), x A^(x)}, which have the same action; i.e., 

Sffl^AJ = S[\p x ip x A»] ■ (2.6) 
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The Grassmann integration over ip and ip gives the same result independent of X(x) since the corre- 
sponding Grassmannian Jacobian for this transformation is unity. Hence there is an overall volume 
divergence in the functional integration over the gauge field in Eq. (|2.3| ). This is analogous to the 
divergence one obtains in the spacetime integral of a translationally invariant function. The proper 
definition of the measure in Eq. fl2.4|) must ensure that the gauge field integration extends only over 
gauge-inequivalent configurations (Bailin and Love, 1986, pp. 116-119). 

This problem can be resolved via the introduction of the Faddeev-Popov determinant (Popov, 1983). 
[Historically, this first arose in connection with the correct quantisation of non-Abelian gauge theories 
(Faddeev and Popov, 1967).] The net effect of this procedure in QED d is simply to introduce a gauge 
fixing term into the action of Eq. ( |2.3j ) (Rivers, 1987, pp. 182-185). A common choice for this is the 
covariant gauge fixing term, which we use here, and which gives 

sfatM - ss,^] = sfatM - ^- I d d x {d^f , (2.7) 

where £o is the bare gauge fixing parameter. It should be noted that the difficulty with gauge fixing 
involving Gribov ambiguities does not appear in Abelian gauge theories, [unless one adopts unusual 
nonlinear gauge fixing prescriptions]. It can be shown that physical observables are independent of the 
[renormalised] gauge parameter, £, and, in fact, that they are independent of the actual form of the 
gauge fixing term (Bailin and Love, 1986, pp. 330-334 and references therein). The unrenormalised 
quantum field theory of electrodynamics is then defined by Eq. ( |2.3|) with the action Sj of Eq. (|2.7f ) 
and we may now proceed to obtain the unrenormalised DSEs. 



Unrenormalised Dyson-Schwinger equation for the Photon Polarisation Tensor. The example of the 
derivation of the DSE for the photon polarisation tensor in QED 4 is given in Itzykson and Zuber (1980, 
pp. 476-477) and we summarise these arguments here to illustrate the technique and to emphasise the 
fact that the DSEs are simply the Euler-Lagrange equations of quantum field theory. It should also be 
noted that independent of this functional analysis one can directly obtain these equations by simply 
decomposing the infinite sums of Feynman diagrams for Green's functions in terms of other proper and 
full Green's functions. This Feynman diagram approach is explained in some detail in Bjorken and 
Drell (1965, pp. 283-376). 

Consider the generating functional of Eq. ( |2.3| ) with the action of Eq. (|2.7| ). It can be shown [e.g., 
Itzykson and Zuber (1980, pp. 211-212)] that the generating functional of connected Green's functions, 
G\rJ, r), J^}, is defined via 

z % V, M = exp {Q% rj, J M ]) . (2.8) 

In order to obtain the DSEs one must simply note that, in analogy with the case of standard calculus, 
the functional integral of a total functional derivative is zero given appropriate boundary conditions, 
[see, e.g., Collins (1984, pp. 13-15)]. Hence, for example, 







dfi(?p, ip, A) 
d(Jt(il>,ip,A) 



S 



SA^x) 



expjz (S^,^,AJ + j d<ix [V»V + + A tl J> x 
SS i , t f„\\ ^ /„• (q CT, „/. a l , f jd„ \~J„f , »/,/./ 



SA^x) 



+ J^x) exp Ui SM, 1>, AJ + / d x V ^ + + 



5S e 



SA^x) 



5 5^5 
i5J' ^Sr] , iSr] 



+ J n( x ) \ z [v,V,Jfi} ■ 



(2.9) 
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Differentiating Eq. (|2.7j) immediately gives 



SS e 



SA^x) 



d p d p g ilv - f 1 - — j d,,,d v 



f 



(2.10) 



from which it follows that after dividing through by Z we can write Eq. (|2.9| ) as 




SQ 



iSJ u (x) 



Srj-f \x)^^ 5rjf(x) ' 8r]f(x) 



SQ 



S 



sg 

7 ^r//( 



= -J M (x) . 

(2.11) 

Equation. ( 2.11 ) represents a compact form of the nonperturbative equivalent of Maxwell's equations. 
To illustrate this we use it to obtain an expression for the photon vacuum polarisation. We can now 
perform a Legendre transformation and introduce the generating functional for the connected, one- 
particle irreducible [1-PI] Green's functions, T[ip, ip, A^]: 



g[rj, 7], = ir#, ip, + i / rf d x U V + ^ + a^ 



:2.12) 



An explanation of the fact that T[ip, ip, A M ] generates the 1-PI Green's functions is given in Itzykson 
and Zuber (1980, pp. 289-294). The 1-PI Green's functions are also frequently referred to as the 
proper Green's functions. It follows from the rules of Grassmannian integration that Z [77, 77, J^] and 
hence Q\rj, rj, J M ] depend only on powers of pairs of Tj and 77, which implies that setting Tj = 77 = 
after differentiating Q [or Z] will only give a nonzero result for equal numbers of Tj and 77 derivatives. 
Similarly, in the absence of fermion derivatives it can be seen that only even numbers of derivatives of 
Z and Q with respect to J survive when we set J = 0, which immediately leads to Furry 's theorem 
for QED, [see, e.g., Itzykson and Zuber (1980, pg. 276)]. Throughout this work we diagrammatically 
represent proper Green's functions by unfilled circles and full [connected] Green's functions by shaded 
circles. See, e.g., Fig. [2.1| , where the use of the proper electron-photon vertex, T^, is necessary to 
avoid the double-counting of Feynman diagrams using the arguments of Bjorken and Drell op cit. From 
Eq. ( |2.12j) the following relations immediately follow: 



SG 



i5J»(x) 
ST 



x) 



SQ 



rjf(x) 



iSriHx) 
ST 

Sip (x) 



ip (x) 



rf{x) 



sg 



iSrjf (x) ' 
ST 



(2.13) 



Sipf(x) 



From Eq. ( |2.13| ) we see that we now have expressions for ip, ip, and A^ in terms of 7/, 77, and J M and vice 

versa; i.e., ip a {x) = ip a [Tj,rj, J^{x)] = iSQ\r],r], J^/SrjKx) with the spinor index a now explicitly shown. 
It is now easy to see that setting ip = ip = after differentiating T gives a nonzero result only when 
there are equal numbers of ip and ip derivatives in analogy to the case for Q. Using Eq. Q2.13] ) and the 
fact that Sipa/Sipp = S a /35f g 5 d (x — y) we find that 



I 



i / d d z 



s 2 g 



S 2 T 



5r)l(x)rf(z) Sip^(z)ipp(y) 



7] = 7] = 

ip = ^p = o 



S rs S fg S d (x - y) . 



(2.14) 



Hence it follows that, when the fermion sources (77,7/) vanish, Eq. ( |2.11| ) can be written as 



ST 



SA^(x) 



d p d"g 



[IV 




A>{x) -iDeftr^i,^]) 
/ 



(2.15) 



7 



where we have made the identification 



S f (x,y, [A,]) 



I 



5 2 r 



Sip (x)5i/jf (y) 



(2.16) 



which is the propagator of the fermion of flavour / in an external electromagnetic field A^. This iden- 
tification follows from Eq. fl2.14| ) and the fact that the second derivative of Q is the full 2-point fermion 
Green's function; i.e., the fermion propagator in the presence of a background field, Sf(x, y, [A^]). 
Clearly, the full fermion Green's function (x,y) follows from setting A^ = in Eq. ( |2.16| ). 

To obtain the DSE for the photon polarisation tensor it remains only to act with 8/8A u (y) on Eq. ( [2.15| ) 
and set J M (x) = 0. It can be shown that 
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8 2 T 



6A»{x) ty f (y)8^jf(z) 



;2.17) 



Q=A ll =i>=il> 



is the proper [i.e., 1-PI] fermion-gauge-boson vertex [not to be confused with the generating functional 
r]. Similarly to the fermion case it can be shown that the second derivative of V with respect to A gives 



the inverse photon propagator, (D 
photon propagator 

(D~ l y v (x, y) 



x 



y). Thus from Eq. ( |2.15 ) we obtain the DSE for the inverse 



8 2 r 



8A^(x)8A v {y) 



=i/j =q p=0 



1 

T. 



\x 



y)+U fll/ {x,y) (2.18) 



where we have defined the photon polarisation tensor, II 



f 



If 



d A zi d d z 2 tr 7 At S' / (x, z 1 )Y f v {y\ z u z 2 )S f (z 2 , x) 



(2.19) 



Using translational invariance we can write the photon propagator in momentum space as 



D^(q) 



q 2 + ie l + n(g 2 ) su (g 2 + ze) 2 ' 

where we have made use of the Ward-Takahashi identity [WTI] for the photon propagator, 

^n^(g) = , 



(2.20) 



:2.2ii 



to define Il(g 2 ) by W^q) = (-g^q 2 + q^ 1 q v )Ii(q 2 ) . Note that Il(g 2 ) is independent of the gauge 
parameter £o in QED as a result of current conservation. A discussion of WTIs is given, for example, in 
Itzykson and Zuber (1980, pp. 407-411) and in Bjorken and Drell (1965, pp. 299-303). There is a whole 
family of such identities, which follow from the fact that gauge invariance leads to the conservation of 
the electromagnetic [e.m.] current. Because renormalisation does not affect e.m. current conservation 
the WTIs also apply to the renormalised Green's functions of the theory provided that the regularisation 
of infinities is performed in a gauge-invariant way. That H.^ is transverse leads to the fact that the 
photon remains massless, [see, e.g., Itzykson and Zuber (1980, pp. 318-329]. The cases £ = and 1 are 
referred to as Landau and Feynman gauges respectively [and similarly for £ after renormalisation]. The 
normalisations in these equations are such that at lowest order in perturbation theory we have II = 0, 
as well as 



Fl(y,z 1 ,z 2 ) =l v 8 d {y -z 1 )8 d (y- z 2 ) and (z - mjf) S f (x, y) = 8 d (x - y) . 



(2.22) 
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in 




bO'W'Q^VA, ZZ VWV + VW\Qw0AA/ 

zD iD ° iD o /n i'D 

Fig. 2.1. The Dyson-Schwinger equation for the photon propagator. 

It should be noted that, once eg is factored out, there is no explicit flavour dependence in the proper 
vertex since the same sum of Feynman diagrams will contribute in each case. We will later see that 
once we define the renormalised charge through appropriate boundary conditions a flavour-dependence 
still arises for nondegenerate fermion masses. 

We have seen that second derivatives of the generating functional, T[tp, if), A^], give the inverse fermion 
and photon propagators and that the third derivative gave the proper photon-fermion vertex. In general, 
all derivatives of the generating functional, T[ifj,ip, A^], higher than two produce the corresponding 
proper Green's functions, where the number and type of derivatives give the number and type of proper 
Green's function legs. 



The DSE for the photon propagator is represented diagrammatically in Fig. [2.1| . Part a) of the figure 
represents Eq. ( 2.19| ) for the photon polarisation tensor and part b) shows the photon propagator, 



D^ v (x, y), defined in terms of IF"" (a;, y) and the bare photon propagator, Dq U (x, y) [given by Eq. ( ]2.20| ) 
with II = 0] and corresponds to the inverse of Eq. ( |2.18| ). The momentum-space representation of the 
DSEs is readily obtained by either Fourier transforming the coordinate-space form or, more readily, by 
using the standard rules for Feynman diagrams based on the lowest-order perturbative contribution to 
the nonperturbative quantities. For example, for the photon polarisation tensor we obtain 

= (-l)E(^) 2 / 7^rtr[(z 7At )(z^(£))(zr/(£,£ + g))(z^(£ + g))] , (2.23) 



where the factor of (—1) arises from the fermion loop in the usual way. We have elected to not explicitly 
show such factors of (—1) in the diagrammatic representation. In momentum space, Fig. |2.1| b) corre- 
sponds to iD pl/ (q) = iDQ T (q){5" + iU Tp (q)iD pu (q)], which can be obtained from the Fourier transform 
of Eq. (|2.20|) . Throughout this work we use figures like Fig. ^j] and the associated notation so that the 



momentum space equations are easily written down by inspection. 



Unrenormalised Dyson-Schwinger equation for the Fermion Self Energy. Following a procedure sim- 
ilar to that above (Itzykson and Zuber, 1980, pp. 478-479) one may derive an integral equation for the 
fermion propagator starting from 







d^,^A)-^expl [ i(s^,'i/j,A tl ) + J d d x [/r/ + r7V + A^])} 
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iS 

Fig. 2.2. The Dyson-Schwinger equation for the electron propagator. 



f SSt 



[ Sip(x) 



5 _6 5_ 

iSJ' iSrj' i8r] 



+ r}(x) \ Z[r), 77, J M ] . 



(2.24) 



After differentiating with respect to 77 and setting all sources to zero [77 = 77 = J = 0] we can rewrite 
Eq. flpD as 



(7 0-m f o )S f (x,y)-i(e f o f j d d Zl d d z 2 d d z 3lfl D^(x, Zl )S f \x, z 2 )Tl(z 1 ; z 2 , z 3 )S f \z 3 ,y) = 5 d (x - y) , 

(2.25) 

where D^x, y) is the photon propagator which couples Eq. ( |2.25| ) to Eq. ( |2.19| ). So, one sees that the 
equations for the 2-point functions are coupled to each other and that both also depend on the 3-point 
function, T^. This is the first indication of the general rule that the DSE for an n-point function is 
coupled to other functions of lesser and the same order and to functions of order (n+1) and (n+2). 

The structure of Eq. ( |2.25| ) allows one to rewrite it in terms of the fermion self energy, —iY}'(x,y), 
defined such that 



(i - m{) S*{x, y)- [ d d z, &(x, Zl ) S^z u y) = 5 d (x - y) 



and hence satisfying 



-i&(x,y) = {e{f J d d z 1 d d z 2%l D^(x,z 1 )Sf(x,z 2 )rl(z 1 ;z 2 ,y) 



(2.26) 



;2.27) 



The equation for the fermion self-energy is represented diagrammatically in Fig. |2.2| a) while part b) of 
this figure shows the definition of the fermion self-energy (— i£* (p)) in terms of the fermion propagator 
Sf(p) with S$(p) = l/{tf — 777,0) the bare fermion propagator. Again, the momentum-space form for the 
proper fermion self-energy (— i"E->) is easily obtained from Fig. |2.2| a) using the usual Feynman rules [or 
equivalently from Fourier transforming Eq. ( |2.27| )1 and can be written as 



d d l 



(2tt)' 



-(i lfl )(iSf(£))(iD^(p-£)){irl(£,p)) . 



(2.28) 



From Fig. |2.2|b) or, equivalently, from Eq. ( |2.26| ), we can solve for S^(p) to give S^(p) = 1/[(Sq) 1 
Zf(p)] = l/y-m f -Zf(p)}. 



Unrenormalised Dyson-Schwinger equation for the Fermion-Photon Vertex. This equation can be de- 
rived in a similar way. For completeness, we present it here in momentum space where it is most concisely 
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iS 




Fig. 2.3. The Dyson-Schwinger equation for the electron-photon proper 
vertex. 

written (Bjorken and Drell, 1965, pp. 291-293): 

*r£(p',p) =<7m + E/ ^^(p' + im^ip' + ^p + imtS^p + i^K^ip + ^p' + ej) , (2.29) 

where X is referred to as the fermion-antifermion scattering kernel. The diagrammatic representation of 
this is given in Fig. 2.3. Clearly, T M is coupled to the fermion 2-point function [the fermion propagator S] 
and also to a fermion-antifermion scattering amplitude denoted by M, a 4-point function, which itself 
satisfies an integral equation (Bjorken and Drell, 1965, pp. 293-298) and this illustrates the general 
rule again. The amplitude M is 1-PI with respect to the fermion lines and does not contain any 
fermion-antifermion annihilation contributions [i.e., no intermediate single photon state] since these 
would not be 1-PI contributions to with respect to the photon line. In an obvious shorthand 
notation we see from Fig. 2.3 that M = K + K{iS) 2 K + K(iS) 2 K(iSfK + ■■■ = K + K{iSfM. 
Hence it follows that the fermion-antifermion scattering kernel K has no annihilation contributions 
and is 2-PI with respect to the fermion-antifermion pair of lines. The latter requirement ensures that 
there is no double counting when (iS) 2 K is iterated to form M. It then follows that to lowest order in 
perturbation theory M? 9 = K^ 9 = 5^ 6 '(eg) 2 (27 il ){iDQ V ){i^j u ), where D is the bare photon propagator. 
The simplest so-called ladder approximation consists of approximating M by iterating the lowest order 
perturbative contribution for the kernel K together with the replacement of the fermion propagators 
by their perturbative form; i.e., Sq(p) = l/[^ — m{,]. 

Just as there was a WTI for photon vacuum polarisation [i.e., g M IP^(g) = 0] there is a WTI which 
relates the fermion propagator to the proper fermion-photon vertex [see, for example, Itzykson and 
Zuber (1980, pp. 407-411)] 

k,T^(p + k,p) = (Sfy\p + k) - (Sfy\p) . (2.30) 



Renormalisation of the Equations. Once the technique for deriving the unrenormalised equations is 
known then the renormalised equations follow. Essentially, one need only modify the action of the 
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theory to include the necessary counterterms and repeat the above procedure. For the purposes of 
this discussion we will temporarily adopt a notation similar to that used in Bjorken and Drell (1965, 
pp. 283-376), where the renormalised action and the quantities derived from it are denoted by a tilde 
to distinguish these from the unrenormalised quantities. 



In Nf flavour QEDd one has the following [infinite] renormalisation constants (Dyson, 1949): the fermion 
wave function renormalisation, Z2, which relates the bare fermion field to the renormalised one; the 
photon wave function renormalisation, Z 3 ; and the vertex renormalisation, Z±, which relates the bare 
and renormalised charge: 



f 



'Z«An ; e 



Z{, 



J 



2.31 



A fundamental requirement and consequence of gauge invariant regularisation and renormalisation of 
the theory is the Ward Identity (Ward, 1950; Takahashi, 1957), which leads to 

Z{ = z{ . (2.32) 

This identity together with the multiplicative nature of renormalisation ensures that the WTIs are pre- 
served in the renormalised theory. In addition it ensures charge universality; i.e., that different species of 
fermion have their charge renormalised by the same multiplicative factor. The QEDd renormalised action 
[i.e., including the counterterms] is obtained from Eq. ( |2.1[ ) by identifying S^[ip,ip, A 11 ] = S^[ip , ipo, Aq] 
and making use of Eq. (|2.31|) . Hence the renormalised action is given by 



S^,ip,A^ 



N f 



I 



d d x 



N f 



+£[(4 
/ 



/=i 



/=1 



+ e f 4) V 



4> f 



2£ 



1) 



(2.33) 



where the four terms in the last line of this equation are called counterterms and depend on both the 



regularisation parameter and the renormalisation point. We see that 5m^ = m B - 
and Se? = e B — = [Z{ — l)eA There is no gauge-parameter counterterm since £0 = Z^, which 



^ = Zlml — m/ 



follows from Eq. (|2.31 ) and from the WTI, q fJi H' J ' u = 0, which ensures that vacuum polarisation only 
affects the transverse part of the photon propagator. We will refer to ip ,ipQ,Ao,ml,eQ, and £0 as the 
bare quantities, although occasionally in the literature this term is also used to describe m B and e B . 
This is discussed for example in Collins (1984, pp. 10-11). The renormalised [i.e., physical ] fermion 
masses and charges are vrJ and e' , respectively, and £ is the renormalised gauge parameter. 



The renormalised action, S%, will lead immediately to the corresponding renormalised generating func- 
tional Z, Q, and T, from Eqs. (|2.3|) , Q2.8|) , and fl2.12|) respectively. Using the same steps as before 
these in turn lead to the renormalised Green's functions D^ v , S r * , V* , K$ 9 , etc. The renormalisation is 
carried out subject to the following boundary conditions on renormalised quantities: 



(s f r\ P ) 



--(mf) 2 



— m J 



r 



n(o) = 



(2.34) 



together with the definition of the renormalised fermion and photon propagators: 
{S f )-\p) = p'-m 1 -t f {p) , 

-cT + (<f<f/V + ie)) 1 t q»q v 



D^(q) 



q 2 + ie 



1 + n(g 2 ) % 2 + ie) 



;2.35) 
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The first condition in Eq. (|2.34j ) specifies that the renormalised fermion propagator, S^(p), has a pole of 
residue one at the physical fermion mass m?; i.e., that E'(p) = at p 2 = (mf) 2 . The second condition 
ensures that when an on-shell fermion of flavour / is probed with a zero-momentum photon we measure 
the physical charge, e*. The notation p 1 = uJ means that all factors of p 1 are replaced by m*, which 
is equivalent to having act on a free spinor at the fermion mass pole, p 2 = (m^) 2 . Note that for 
nondegenerate fermion masses [i.e., m* 7^ vn/\ this boundary condition introduces a subtle distinction 
between and , since otherwise these vertices would be identical. The last boundary condition 
ensures that an on-shell photon has a pole at p 2 = with unit residue after vacuum polarisation 
corrections. 



There are some subtle difficulties with these boundary conditions due to infrared divergences but these 
are readily taken care of by introducing a small photon mass as discussed, e.g., in Itzykson and Zuber 
(1980, pp. 413-414). It should be noted that this choice of boundary conditions corresponds to the 
particular choice of the on-shell renormalisation point, which is the usual choice for QED. We know, 
however, from renormalisation group arguments, that this choice is in fact arbitrary. Since QCD is 
believed to be confining, it has no corresponding natural choice of renormalisation point. 



Following these arguments gives, in place of Eq. ( p. 17] ), 

S 5 2 t 



V 1 z ) 5A^{x)s^ f ( y )5ipf(z) 
Similarly, in place of Eqs. ( |2.18|) and ( |2.19| ), we find 



(2.36) 



0=A Al =V=V 



5 2 T 



5A^(x)5A u (y) 



6> 



6 d (x-y)+W fMU {x,y) 



Z 3 [d p d<>g^ - S d (x -y) + -d^dj^x -y) + II' (x, y) 



where we have [c.f., Eq. (|2.19|) 1 



n^(a?,2/) = i^2z((e f ) 2 / d d z 1 d d z 2 tr 7^ \x, z x )Y v (y] z u z 2 )S j \z 2 , 
f 



x , 



(2.37) 



(2.38) 



In momentum space we then define II' (q) = (—g^uQ + QhQv)H-'(q ) an d hence it follows from the 
momentum space form of Eq. Q2.37 ) and from Eq. (|2.35 ) that II(g 2 ) = (Z% — 1) + U'(q 2 ). From the 
boundary condition 11(0) =0 we find that 



z 3 = 1 - n'(o) 

and so finally that 

ifl^(q) = i(-g, v q 2 + q^Mq 2 ) = l(-g, v q 2 + q»q„W(q 2 ) - n'(0)] . 



(2.39) 



(2.40) 



Note that q^H'^ = [assuming a covariant, gauge-invariant regular isat ion scheme] just as before, 
and hence we have seen that the WTI for the vacuum polarisation has survived renormalisation; i.e., 
g A1 IP iy = 0, as expected. 



Repeating the same arguments for the fermion self-energy we find Eq. (|2.25 ) becomes 
Zl(i0-m f Q )Sf(x,y)- [ d d z 1 ^'(x,z 1 )S^z 1 ,y) = 5 d (x-y) , 



(2.41) 
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where 



- iZf'{x, y) = Z{(ef ) 2 J d d z x d d z 2 ^D^(x, Zl )&(x, z 2 )f z 2 , y) 



[2A2) 



E*(p). Using Lorentz 



In momentum space we have (S?) l (p) = Z 2 {p" — ml) — S^'(p) = i> — m) 

invariance allows us to decompose YJ 1 into Dirac-vector and -scalar pieces: Yj'{v) = Y a ' (p 2 ) $ + Y?' (p 
and similarly for &. Hence we find Y f d (p 2 ) = Y f d '(p 2 ) - [Z{ - 1) and Y{{p 2 ) = Ef (p 2 ) + (Z 
Hence the boundary condition E^(p) = at p 2 



m J 



\m*) 2 gives 



4 = 1 + E^((m') 2 ) 



f 

m 



Then we have Y^(p) = Tj d (p 2 )p'+ E{(p 2 ), where 



U{P 2 



-E{'((m 



A2\ 



/Z 



E?((m 



(2.43) 



(2.44) 



If we introduce the notation {S f )~ 1 {p) = A f {p 2 )i> - B f (p 2 ) = (Z f y l {p 2 )[p' - M f (p 2 )}, then this leads 
to the renormalised DSE for the fermion propagator in the form that we use in later sections: 



Af(p 2 ) = (Zf)-\p 2 ) = 1 - t{{p 2 ) , Bf(p 2 ) = {Zf)-\p 2 )Mf{p 2 ) = m f + t{{p 2 ) 



(2.45) 



For the fermion-photon proper vertex we have, in place of Eq. (|2.29|) , 



Ktp'ip) 



J TT^d^V + ?)rj(p' + Q,P + l)S 9 (P + q)K f9 (p + q,p' + q, q) 



z ' 7 " J {2*)* 

z[ ltl + Al'i P , , P )= 1 , + Al( P ',p) 



(2.46) 



which also defines both A^' and A^. Hence, we have A^(p',p) = A^{p\p) + (Z{ — 1)7^ and from the 
boundary condition for the vertex in Eq. ( p.34j ) we find 



which gives finally that 



(Zf-1)7 M = -Ai'(p,p) 



Al(p',p) = A{:(p',p)-Af'(p,p) 



(2.47) 



(2.48) 



Since (e^TjJ, (D 1 ) M ^, and (S-f) 1 are all obtained by differentiating T[ip, ip, A M ], then it follows, for 
example, from Eqs. ( gT7||2T3l , |2T33| , |2T36|) that (e£r£) = (Z 2 / v ^)- 1 (e / f^). Similarly, S f = Z S 2 S S and 



D^y = Z^D^ V . Thus from the unrenormalised WTI of Eq. ( p. 30 ) we find for the renormalised quantities 



{Z{/Z{)k^T^(p+k,p) = (S / )- 1 (p+A;)-(5 / )- 1 (p), which implies that since the theory is renormalisable 
then Z( and Z 2 can only differ by a finite amount. In fact from the boundary conditions in Eq. ( [2.34j ) 
we see immediately that we must have Z{ = Z 2 , as was stated earlier. The WTI implies the Ward 
identity, which can be written as r^(p, p) = d(S^)~ 1 (p)/dp 11 . [For ease of comparison with Bjorken and 

Drell (1965), pp. 311-312, note that, e.g., A and n used here correspond to A c and e 2 n c , respectively, 
in that work.l 



In repeating the previous arguments and deriving the above renormalised equations it soon becomes 
apparent that Z-factors are only introduced when bare propagators or vertices appear in the DSE. 
We have not discussed details of the various regularisation schemes but it should be clear that all of 
the regularisation-parameter dependence has been absorbed into these factors. Let A denote a generic 
regularisation parameter such that A — > 00 as the regularisation is removed, e.g., A might denote a 
simple momentum cut-off or using dimensional regularisation, where d = 4 — e, we would identify 
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A = l/e, etc. Then the Z-factors are functions of the regularisation parameter and the renormalisation 
point; i.e., Z 2 = Z 2 (A, fi), etc. For m + n > 3, where m is the [even] number of external fermion legs and 
n the number of external photon legs, it can be shown that the unrenormalised Greens functions can be 
decomposed into an infinite sum of skeleton diagrams each of which is made up only of combinations of 
S?, D^ u , and eor^j. However, we have seen that in the renormalised theory the quantities e^T^, D^, and 
are simply replaced by the renormalised quantities e^T^j, D^, and S*. Hence, for Green's functions 
with m + n > 3 we simply replace the skeleton expansion for the unrenormalised Greens function by 
the same skeleton expansion with the renormalised quantities. This is the meaning of K in Eq. (|2.46[ ), 
for example. This completes the renormalisation program for QED. A more detailed discussion of the 
skeleton expansion and associated proof of renormalisability can be found in Bjorken and Drell (1965). 



Review. We have illustrated the general technique for obtaining the DSEs for a field theory and 
have demonstrated by example the nature of the hierarchy of equations; i.e, that the equations couple 
a given n-point function to adjacent n-point functions in an infinite tower of coupled equations. We 
have also discussed in some detail the renormalisation of these equations. As was already noted, choices 
of subtraction point [i.e., renormalisation point] other than the on-shell subtraction point of Eq. ( |2.34| ) 
are possible. From the previous discussion it is clear that such a change can simply be absorbed as a 
finite rescaling of the renormalised charges and masses whilst leaving the physical content of the theory 
unchanged. This observation is the basis of renormalisation group studies of field theories. We note 
in passing that with a suitably unconventional choice of subtraction scheme one can find Z{ + Z£, 
although this is not the case in standard approaches [of course, this has no physical significance]. The 
renormalised masses and charges in QED 4 only correspond to the physical ones when the on-shell 
renormalisation point is used. There is no corresponding on-shell renormalisation in QCD, since QCD 
is believed to be a confining theory. We have now introduced all of the concepts that are necessary 
for the discussion of DSE studies of QED d in Sees. |] and |. For our discussion of QCD in Sec. pT2 
simply point out the added complications due to the non-Abelian nature of the theory. 



we 



2.2 Quantum Chromodynamics [QCD] 

QCD is a non-Abelian gauge theory [based on the fundamental representation of the group 577(3)], 
which leads to significantly different behaviour than that shown by the Abelian theory of QED [based 
on the group £7(1)]. In simple terms this arises because the non-Abelian gauge fields are self-interacting. 
The basic Lagrangian of QCD is given by 



S[i>,ijj,A^ 



d d x 



2 / (ty-ml) tf- - 



N 



- 1 fiua- 1 a 



[2A9) 



where 



K v = d»A:-d»AZ + g f abc AZA» 
D» = d»-ig AH a = d»-ig A< 1 



(2.50) 
(2.51) 



N 



f 



and where A 11 = A%t a . The fields A% are the gluon fields and tp* are the quark fields, where / = 1, 
is the quark flavour index. The indices a = 1, • • • , 8 are colour indices and summation over repeated 
indices is to be understood. We use standard conventions [e.g., Yndurain (1993), Muta (1987)] where 
t a = A a /2 and the hermitian [Gell-Mann] matrices A a are the generators of SU(3). The A a matrices 
satisfy tr(\ a \ b ) = 25 ah and [A a , X b ] = 2if abc X c , which leads to tr(t a t b ) = T F 5 ab = 5 ab /2 and [t a ,t b ] = 
fabctc = —C abc t c . [Note that Itzykson and Zuber (1980, pp. 563-567) define instead antihermitian 
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matrices t a = iX a /2]. We can write, e.g 
on the values Ca = N, Cp 



(l/4)FjTF t 



(l/2)tr(F' iV F /iI/ ), where F^ = t a F£ v 



(if g )[D^, D u \. For SU(N), the Casimir invariants Ca,Cf, and T F take 
(N 2 - l)/2N, and T F = 1/2. 



As for QED, it is necessary to remove the infinite gauge- volume problem due to integrations over gauge- 
equivalent field configurations. This is remedied by the prescription of Faddeev and Popov (1967) and 
has the added complication that in order to maintain gauge invariance and unitarity in covariant 
gauges it is necessary to introduce unphysical auxiliary fields (uJ a , ui a ). These are anticommuting spin- 
zero fields called ghost fields and it is also necessary to integrate over these fields in order to define the 
generation functional Z for QCD. A technical point, which we shall neglect henceforth, is the fact that 
the Faddeev-Popov gauge-fixing procedure only uniquely specifies the gauge field configuration with 
respect to infinitesimal local gauge transformations. Under finite transformations gauge-equivalent field 
configurations remain and are referred to as Gribov copies [see, e.g., Itzykson and Zuber (1980, pp. 
574-582) and Gribov (1979)]. It is possible to choose ghost-free gauge-fixing prescriptions such as the 
axial gauge, where n ■ A = for some spacelike four- vector n M , (n 2 < 0). The cases n 2 > and n 2 = 
are referred to as the temporal and light-cone gauges, respectively. The disadvantages of such a choice 
are that Lorentz invariance is now broken explicitly and unphysical singularities are present in the 
denominators of momentum-space propagators; i.e., when p ■ n — 0. 



The quantum field theory associated with Eq. ( p.49|) is defined for covariant gauges by the generating 
functional 



Z[r],V,Jn,t,£,} = 

J dfi(ij),ip, A,uJ,u>) exp (iS^[ijj, ip, A M , 57, u>] + i / <7".r 

where 



(2.52) 



S^ip, ip, u, u] = S$, ip, A h 



d d x 



(d^ a )(6 ab d» - gofabcAfiub - — (d,A») 2 



(2.53) 



where ^ is the bare gauge fixing parameter. The Faddeev-Popov technique for non-Abelian gauge fields 
[as well as the gauge fixing procedure for QED] can be understood in terms of the more general class of 
Becchi-Rouet-Stora [BRS] transformations. These are a generalisation of the local gauge transforma- 
tions to include the ghost and gauge-fixing terms. The Slavnov- Taylor identities [STIs] can be derived 
from the BRS invariance of QCD and correspond to the WTIs of QED. It can be shown that requiring 
BRS invariance of a gauge theory automatically generates both the ghosts [where appropriate] and the 
gauge fixing as well as ensuring the gauge-invariance of physical observables, [see e.g., Pokorski (1987) 
pp. 64-85 and Pascual and Tarrach (1984) and references therein]. 



Just as for QED we introduce appropriate counterterms into the QCD action in Eq. ( |2.53 ) in order to 
define the QFT and then apply appropriate boundary conditions at the chosen renormalisation point 
[a momentum scale that we refer to as Since the renormalisation point, fi, is arbitrary, physical 
observables must be independent of its choice. In other words, hadronic masses and cross-sections 
must be renormalisation-point independent even though individual Green's functions [e.g., the quark 
and gluon propagators] may not be. In analogy with QED we identify the renormalised action using 
S>£ [if), ip, 57, uj] = S^[ip ,ipo, Aq,uj ,u! } together with the relations [see, e.g., Yndurain (1993) pp. 45- 
83] 

i)l = y/Zp'ip 1 ; A% = y/Z^A* ; u a0 = \[Z^u) a ; g = Z g g . (2.54) 

As a consequence of gauge invariant regularisation and renormalisation of the theory the STIs are 
maintained, which leads to the fact that all quark colours have the same Z F and that all gluons and 
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Fig. 2.4. The Dyson-Schwinger equation for the quark self-energy. 




Fig. 2.5. The Dyson-Schwinger equation for the gluon propagator. 
[Here and below the broken line represents the propagator for the ghost 
field.] 

ghosts have the same Z B and Z u respectively. In analogy with QED we see that £ = Zb£, since 
only the transverse part of the gluon propagator is modified by vacuum polarisation [i.e., q fJ ,Tl fJ ' 1 ' = 0]. 
Here, Z g plays the role of the QED combination Z[j 'Z^y/Z^, so that in place of Z{ in QED we have 
Z[ = ZgZpy/Zg in QCD. The STIs also imply that the same renormalisation constant Z g applies to 
the quark-gluon, ghost-gluon, three-gluon, and four-gluon vertices. The STIs are the reason that we do 
not need other independent renormalisation constants for these couplings, [see, e.g., Itzykson and Zuber 
(1980, pp. 593-594) and Muta (1987, pp. 158-179)]. We can define 5m f = m s B - m f = Z s F m{ - m f = 
(Z F Zf m — 1)772/ as for QED, where the last result follows from the definition of Z^ by ml = Z^m) '. 



The derivation of the unrenormalised and renormalised DSEs proceeds in an analogous way to that for 
QED and, as already stated, these are a direct result of the BRS invariance of the theory. In Fig. |2]4] we 
show the DSE for the quark self-energy. The graphical representation of the quark propagator DSE is 
the same as that for the electron given in Fig. |2.2| b). Figure 2.5 specifies the Dyson-Schwinger equation 
for the gluon propagator and can be compared with the photon DSE in Fig. |2.1| . The symmetrisation 
factors of 1/2 and 1/6 arise from the usual Feynman rules, which also require a negative sign [unshown] 
to be included for every fermion and ghost loop. In Fig. |2l^ we show the quark-gluon proper vertex 
DSE. The notation is analogous to that used for the QED fermion-photon vertex in Fig. 2.3; i.e., the 
amplitudes M, M' , M" , and M'" are 1-PI with respect to all external legs and do not contain any 
single-gluon intermediate states. 



The difficult challenge is to find an Ansatz for the unknown renormalised propagators, vertices, etc., of 
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Fig. 2.6. The quark-gluon vertex Dyson-Schwinger equation. 



QCD which satisfy the DSEs and respect the STIs of the theory. Progress can be made by satisfying a 
subset of the DSEs and STIs and then supplementing these with information gleaned from lattice gauge 
theory, phenomenology, etc. The test of such a scheme is to enlarge the set of DSEs considered and look 
for a decreased dependence on the phenomenological input. Since the most significant difficulties in 
QCD arise in the gluon sector this is the natural place to appeal to phenomenology while attempting to 
construct explicit solutions and maintain important symmetries in the quark sector. Two symmetries 
of particular importance are gauge invariance and chiral symmetry. 



Note: Having completed the discussion of renormalisation, we now dispense with the tilde notation 
for renormalised Green's functions, self-energies, vacuum polarisations, etc. In analogy with the QED 
case we then denote the renormalised quark and gluon propagators [in momentum space] by S'{jp) and 
D^{p) = SabD^ij)) respectively, (a, b = 1,...,8 are SU(3) colour indices. For a covariant gauge the 
renormalised gluon propagator and quark propagators have forms identical to those in Eq. ( |2.35| ). We 
can write for the proper quark-gluon vertex T^(p',p) = t a T' J '(p',p). 

The renormalised quantities S?, £•*, D, II, T*, £, g, m' , etc. have a dependence on the both the 
renormalisation point fi and the gauge parameter £ in general, but we do not explicitly indicate these 
for notational brevity. The renormalised coupling g is related to the running coupling constant a s (Q 2 ) = 
g 2 (Q 2 )/4nby 



g 2 = g 2 (Q" 



4vr a s (Q 2 



(2.55) 



where q 2 = —Q 2 > 0. For Q 2 3> Aq CD the running coupling constant is given by the leading-log [i.e., 
one-loop] result 

asW > (33 - 2JV» HQVKcd) ~ HQ 2 1 Ken) 1 ' 

where Nf is the number of quark flavors and Aqcd is the scale parameter of QCD. We will see that 
oIm = 12/(33 — 2Nf) is the anomalous dimension of the mass. It was shown by Symanzik (1973) 
and Appelquist and Carrazzone (1975) that when Q 2 <^ (ttt/) 2 these quarks can be neglected since 
they only give rise to effects of 0(Q 2 jvrJ 2 \ This is sometimes referred to as the decoupling theorem 
[Yndurain (1993), pp. 75-79]. Note that Eq. (|2.56|) is renormalisation-scheme and gauge independent, 
as is its second-order [i.e., two- loop] extension. The [renormalised] gauge-parameter £ is chosen at the 
renormalisation point and, in general, it also runs as the renormalisation scale \x is varied. 



At the renormalisation point \i we impose the boundary conditions 



D™{p) 



-9 av + (1 - 



p°p v 

p2 



(s f r\ P ) 



i>-m s 
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tk 2 C 



r 



(2.57) 



where = (p' — p)^. These are to be understood to imply that U(p 2 ~ —pi 2 ) — 0, Y,f(p 2 ~ —pi 2 ) — 0, 
Zf(p 2 ~ —p 2 ) — 1, and M^(p 2 ~ — /i 2 ) ~ using the (Z,M) representation of the fermion prop- 
agator in Eq. (|2.45|) . Here we have introduced the notation to denote the "running" mass at the 
renormalisation point p. In a nonconfining theory such as QED the mass would only correspond 
to the physical mass when (mjj) 2 = —p 2 such as when the on-shell renormalisation scheme is used. 
The boundary condition on A^ u (p',p) should obviously be chosen so as not to violate the STI relating 
the quark propagator and the quark-gluon vertex [see later] but since QCD is an asymptotically free 
theory then the boundary condition given above is certainly approximately correct for p sufficiently 
large. A^ u (p',p) ~ at p' 2 ~ p 2 ~ k 2 ~ — p 2 . Note that we have used approximation symbols in 
specifying Eq. (|2.57|) since the exact boundary conditions depend on the detailed choice of renormali- 
sation scheme, e.g., MS or MS if we use dimensional regularisation. In general there is some residual 
renormalisation-scheme dependence as well as the obvious renormalisation-point dependence for the 
renormalised quantities. Within a particular renormalisation scheme the renormalisation-point inde- 
pendence of physical observables gives rise to the renormalisation group equations, [see, e.g., Muta 
(1987) and Yndurain (1993)]. 



We will now also omit the explicit flavour index on the quarks and restore it only when necessary. The 
STI for the gluon propagator can be written as q^JX^^q) = where U^l(q) = 5 a b^ v {q) and IP^g) = 
[ — g^ u q 2 + g M g 1/ ]II(g 2 ), just as was the case for the photon. Equivalently, we can write q^D^ (q) = £. 
The STI for the quark-gluon vertex can be written as (Marciano and Pagels, 1978, p. 172) 

k a r°(p',p)[l + b(k 2 )] = [1 - B(k,p)]S-\p') - S-\p)[l - B(k,p)} (2.58) 

where p' = p + k. The ghost self-energy and ghost-quark scattering kernel are denoted b(k 2 ) and B(k,p), 
respectively. Note that in Landau gauge B(k,p) ~ at the renormalisation point (p 12 ~ p 2 ~ k 2 ~ —p 2 ) 
to lowest order in perturbation theory. Without ghosts Eq. (|2.58| ) is identical to the corresponding WTI 
of QED. In general, other symmetries will also give rise to STIs. In the limit that there are no explicit 
chiral symmetry breaking [ECSB] quark masses, QCD is a chirally symmetric theory and has a STI for 
the proper vertex for an isovector axial- vector current [e.g., a W-boson] coupling to a quark, 

KT£(p',p) = T - [S-\p') l5 + 75 S-^p)} . (2.59) 

Since ghost terms do not contribute this actually has the form of a WTI. In the presence of ECSB quark 
masses the chiral symmetry WTI is replaced, for flavour /, by [see, e.g., Itzykson and Zuber (1980, 
pg. 557)] 

k a rt(p>,p) = [(S^ V)7s + l^r'ip)] + W T f 5 (p',p) . (2.60) 

The additional term contains the ECSB quark mass m* and the quark pseudoscalar isovector vertex 
f 7 

QCD is well known to be an asymptotically free theory, which means that if the characteristic momenta 
in some physical process are sufficiently large and space-like then QCD behaves as a free theory with 
logarithmic corrections. Let Ti(q,g;n) represent some renormalised Green's function, where g is the 
renormalised coupling at the renormalisation point p, q u is a large spacelike four-momentum character- 
istic of the momenta of the external legs, and where we define, as usual, Q 2 = —q 2 > 0. In order to use 
perturbation theory effectively we also need Q 2 ~ p 2 3> Aq CD and so, as stated earlier, g is as defined 
in Eq. (|2.55|) . If we scale all of the external momenta by the factor A then the characteristic momentum 
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will also scale; i.e., q v — ► q w = Ag y and Q — > Q = — q = X Q . The meaning of asymptotic freedom 
is that to leading order in a s we have (for Q' 2 ^> Q 2 ^> Aq CD ) 



Ti(Xq,g;fi) 



(AfT^, (-In A 



(2.61; 



where we have defined g = g(/i 2 ) and g' = g(X 2 fi 2 ) and where g 2 {Q 2 ) was defined in Eq. ( |2.55| ). The 
exponent n, is the naive [canonical] dimension of Tj and the exponent di is defined as its anomalous 
dimension, [see, e.g., Yndurain (1993), pp. 66-75, and also Altarelli (1982) and Reya (1981)]. For 
example, for Q 2 = —q 2 ^> Aq CD we have the leading-log result for the running mass 



M(-Q 2 



m 



l\n(Q 2 /A 2 QCD 



d M 



(2.62) 



where rh is the renormalisation group invariant mass parameter and du = 12/(33 — 2Nf) is the anoma- 
lous dimension of the mass. Note that the asymptotic behaviour of M(p 2 ) when m / is gauge 
independent. The parameter rh sets the scale of the ECSB in QCD and is analogous to Aq C d- To 
lowest order the renormalised ECSB mass, m M , is related to rh using Eq. ( |2.62| ) by 

/ r l i du 

77v = M(-/i 2 ) = rh/ - ln(/i 2 /Aj co )j , (2.63) 

-q 2 and Q 2 > fi 2 > 

L Jn(Q7A^ D )j , (2.64) 

which gives Z(-Q 2 ) = 1/ 1 1 m(Q 2 / A 2 QOD )] S , where d s = -2f/(33 - 2N f ) is the [gauge-dependent] 
anomalous quark dimension. Note that Z(—Q 2 ) = 1 to leading order in Landau gauge [£= 0]. Similarly, 
the asymptotic behaviour of the transverse part of the gluon propagator at leading order is 

-dr, 



where the renormalisation-point dependence is explicitly indicated. For Q 2 
Aq CD the inverse quark propagator is, to lowest order, 



S-\q)= Z-\q 2 )i =i 



1 



q^q v ~ 


1 ■ 




q 2 . 



-HQ 2 /A 2 QCD ) 



(2.65) 



where do = (39 — 9£ — 4Nf)/ [4(33 — 2Nf)] is the gluon anomalous dimension. The quark-gluon vertex 
for Q 2 ee — p 2 and Q 2 > /i 2 > A 2 QCD is given by 



r CT (p',p) 



7 



l -HQ 2 /A 2 QCD f 



(2.66) 



where ci r = -(27 + 250/ [ 8 ( 33 - 2N f)}- lt is straightforward to verify that 2d s - 2d r + d D = 1/2 
which, since Z g = Z^/Zf^Zb, implies, for example, that the /3-function is independent of the choice 
of gauge to leading order, as already stated. It is also possible to similarly define a running gauge 
parameter £(Q 2 ), where £ ee ^{Q 2 = fi 2 ). These results all follow from asymptotic freedom and the 
renormalisation group equations. It should be noted that the asymptotic spacelike behaviour of the 
various quantities [e.g., S(p), Z(p 2 ), M(p 2 ), D ua (q), II(g 2 ), and T v (p',p)] has no explicit dependence on 
the renormalisation point /i, although there is an implicit dependence through the renormalisation-point 
dependence of the gauge parameter. 

If we consider the electromagnetic couplings of the quarks, then it follows from electromagnetic gauge 
invariance that the photon-quark vertex satisfies 



k a ri m Xp',p) = s-\p')-S'\ P ) . 



(2.67) 
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This is the same as the corresponding WTI of QED and has no ghost complications. 

In momentum space the renormalised inverse quark propagator is 

S-\p) = ^-m-E(p) = Z-\ P 2 )^-M{ P 2 )] = A(p 2 )^-B(p 2 ) , E(p) ee E'(p)-E'(p) , (2.68) 

where 

-&(p) = \z v g 2 J ^(^(iSmiD^ip-mU^P)) • (2-69) 
We have used J2 A a A a /4 = 4/3. The signal that DCSB has occurred is that the renormalised self-energy 

a 

develops a nonzero Dirac-scalar part; i.e., when £ s Qo 2 ) ^ 0. This leads to a nonzero value for the quark 
condensate and, in the limit of exact chiral symmetry, leads to the pion becoming a massless Goldstone 
boson. 

In general, Aq C d and rh are renormalisation-scheme dependent and so take different values in, for 
example, the modified minimal subtraction [MS] scheme, the minimal subtraction [MS] scheme, and 
the momentum subtraction [MOM] scheme. However, the difference does not appear in leading order 
but only at second-order and higher in a s . If we retain only leading-order terms in the asymptotic 
region we do not need to concern ourselves with this. 



2.3 Euclidean Space Formulation 

Up to this point we have discussed the DSEs in Minkowski space. To an audience unfamiliar with 
studies of lattice gauge theory and the numerical solution of DSEs this would presumably seem natural. 
However, both of these approaches are typically formulated in Euclidean space. 

Our Euclidean space conventions are as follows: 1) we use a non- negative metric for Euclidean four- 
vectors 

4 

a ■ b = 5^ u a^b u = J2 (2.70) 
i=i 

where 5^ v is the Kronecker delta; 2) our Dirac matrices are hermitian and satisfy the algebra 

{7 tl ,7,} = 2^; (2.71) 

and we have 

75 = -7i72 73 74 (2.72) 

so that 

tr[7 5 7a 7^ lu lp\ = ~ 4 e x ^ p (2.73) 

where ex^p is the completely antisymmetric Levi-Civita tensor in d = 4 dimensions. One realisation of 
this algebra is 

7 f = 7° and 7 f = -il 3 , 3 = 1, 2, 3, (2.74) 

where 7 and 7^ can be any one of the commonly used Minkowski space representations of the usual 
Dirac algebra. We note that with these conventions a spacelike four-vector, p M , has p 2 > 0. 
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A straightforward transcription procedure can be employed to determine the action for the Euclidean 
field theory that corresponds to one formulated in Minkowski space: 



M 



4_M 



d x 



4 



-i I d x 

i 7 E ■ d E , 
-ij E ■ A E 
—A E ■ B E . 



(2.75) 

(2.76) 
(2.77) 
(2.78) 



where, as usual, 4 represents Q i w r fuA v M . These transcription rules can be used as a blind implementation 
of an analytic continuation in the time variable, x°: x° — > —2x4 with x M 



x E , etc. 



Employing these rules in Eqs. (|2.1|) and (|2.3|) it is easy to verify that the generating functional for 
Euclidean QED d is 



with the action 

s E ® E ,r,A E ] 



r N f 1 

{ d d x E J2 ^ E {l E ■d E + m f + ie f ol E ■ A E ) ^ E + ~F E V F» V1 



(2.79) 



(2.80) 



[The analogue of this result in QCD is obvious.] 



Equations ( |2.79| ) and ( |2.80| ) are the foundation of Lagrangian lattice studies of QEDd ; as are their 
analogues for other field theories. When formulated on a discrete lattice, the measure in Eq. ( [2.79 ) is 
non-negative and this makes a direct numerical simulation of the theory possible using the standard 
tools of statistical mechanics, which is an important way in which lattice gauge theory benefits from 
being formulated in Euclidean space. 

In addition, the discrete lattice formulation in Euclidean space has allowed some progress to be made 
in attempting to answer existence questions for interacting gauge field theories (Seiler, 1982). In fact, it 
is possible to view the Euclidean formulation of a field theory as definitive [see, for example, Symanzik 
(1969)]. The moments of the Euclidean measure, which can be obtained operationally via functional 
differentiation of Eq. (|2.79|) and setting the sources to zero, are the Schwinger functions: 



S n (x\... 



X 



(2.81; 



which are sometimes called Euclidean space Green functions. Given a measure [i.e., the piece of Eq. ( |2.79| ) 
that does not involve source terms] and given that it satisfies certain conditions [i.e., the Wightman and 
Haag-Kastler axioms], then it can be shown that the Wightman functions, W™(xi, . . . , x n ), can be 
obtained from the Schwinger functions by analytic continuation in each of the time coordinates: 



W n (x 1 ,...,x n ) = lim S n {[x\xl + ix° l },...,[x n ,xl + ix° n }) 



(2.82) 



with x\ < x\ < ■ ■ ■ < x° n . These Wightman functions are simply the vacuum expectation values of 
products of field operators from which the Green functions [i.e., the Minkowski space propagators] are 
obtained through the inclusion of step [8] functions in order to obtain the appropriate time ordering. 
[This is described in some detail in Streater and Wightman (1980, Appendix); Seiler (1982); Glimm 
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and Jaffe, (1987, Chap. 6).] Thus the Schwinger functions contain all of the information necessary 
to calculate physical observables; this notion is used directly in obtaining masses and charge radii in 
lattice simulations of QCD. The same notion can be employed in the DSE approach since the Euclidean 
space DSEs in QED<j can all be obtained from Eq. Q2.79| ) [and the DSEs for other theories obtained 
from their respective Euclidean generating functionals] and the solutions of these equations are the 
Schwinger functions. 



Euclidean — » Minkowski Continuation. In studying DSEs it has been commonplace to obtain DSEs 
from the Minkowski space generating functional and then use transcription rules in momentum space 
[i.e., k° -> ik d and k M -> -k E ] 



M r E 



(fk M -> i J <fk* , (2.83) 

# -> i-f E -k E , (2.84) 
~k E -q E , (2.85) 
k^ -> k E ■ x E , (2.86) 

[the last of which appears in four-dimensional Fourier transforms and follows since x° — > —ix^ and 
x M — > x E ] which are analogues of Eqs. (|2.75 )-( 2.78|) , to obtain equations in Euclidean metric which 



are assumed to be the Euclidean space counterparts of the original equations. This approach is often 
argued to be connected with the "Wick Rotation" (Wick, 1954) and is discussed in most text books in 
association with dimensional regularisation. [In the following we take (k E — > — k E ) so that the inverse 
of the free-fermion propagator is (ry • k + m).] As remarked in Itzykson and Zuber (1980, pg. 485), 
this procedure is easy to justify in perturbation theory when one neglects the possibility of dynamically 
generated singularities in the first and third quadrants of the complex p° plane, however, in connection 
with nonperturbative studies this simple assumption is highly nontrivial and in model studies is often 
incorrect. 



Indeed, in a study in Euclidean space of an approximate DSE for the electron in QED 4 (Atkinson 
and Blatt, 1979) it was found that, instead of a single physical branch point on the timelike p 2 axis, 
the electron propagator had complex conjugate branch points whose position depended on the gauge 
coupling. A straightforward application of transcription rules is invalid in this case. Subsequently 
there have been a number of additional studies of approximate and/or model DSEs for the electron 
propagator in QED 4 (Maris 1993) and the quark propagator in QCD (Stainsby and Cahill, 1990; Maris 
and Holties, 1992; Burden et al, 1992b; Stainsby and Cahill, 1992; Stainsby, 1993; Maris, 1993) which 
have addressed this question; i.e., whether singularities are generated dynamically in the complex p° 
plane which complicate or even prohibit the Wick Rotation. In each of the cases studied so far the 
straightforward transcription has been found to be invalid. Indeed, in one of these studies (Burden et 
ai, 1992b) the quark propagator had an essential singularity in p 2 at timelike infinity which completely 
prohibits the rotation of the p° contour. 

In each of these cases the solution of that Minkowski space equation which is obtained via a straight- 
forward transcription of the Euclidean space equation has no relation to the analytic continuation of 
the solution of the Euclidean space equation. This raises the question of how one may proceed between 
Euclidean and Minkowski space in nonperturbative studies. A possible answer has already been given 
above; i.e., one can define the [model] field theory in Euclidean space and solve for the Schwinger func- 
tions. Assuming then that the field theory "exists", in the sense that all of the necessary axioms are 
satisfied, the Wightman functions are obtained as the analytic continuation of the Schwinger functions, 
Eq. ( j2.82j ), from which one may obtain the Green functions and physical observables. It should be noted 
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that the observed invalidity of the straightforward transcription procedure in momentum space DSEs 
may be entirely due to the limitations of the various model systems that have been studied and that 
no such problem is manifest in the true tower of equations. At the present time there appears to be no 
definitive answer to this question. 

In the following we will adopt the point of view that nonperturbative studies of field theories are most 
easily defined in Euclidean space with physical observables being obtained from the Schwinger functions 
by analytic continuation. As we have described above, this approach is the same as that taken by the 
practitioners of lattice gauge theory and we emphasise that, at least in principle, all physical observables 
can be obtained in this way. 



3 Three-dimensional Quantum Electrodynamics 



It will be seen upon comparison of Sees. |2.1| and 2~2" that the DSEs in Abelian gauge theories are much 



simpler because of the decoupling and effective absence of ghost fields in all [linear] gauge fixing schemes. 
For the same reason the Ward-Identities relating n-point functions take a relatively simple form. This 
makes Abelian gauge theories an ideal framework within which to learn and gain experience in solving 
a field theory using the tower of DSEs. 

In this section we will focus on QED 3 which, in addition to being useful for pedagogical reasons, is an 
interesting theory in its own right. We will work almost exclusively in Euclidean space and use the 
conventions of Sec. |2.3| but we will not explicitly indicate Euclidean space quantities with the superscript 
E: this is to be understood where appropriate. As an illustration we note that in our conventions we 
then have, for example, 

S(p) = - , > , iIV(z/,p) = + *A M (j/,p) (3.1) 

ry ■ p + m + L(p) 

and D,M) = (V - q -f-^j D(q) + £ ^ , (3.2) 

where D{q) = WTWY] ' (3 ' 3) 

[We discuss our QED3 Dirac matrix conventions after Eq. ( |3.7| ).1 

For those concerned with QCD and strongly interacting theories QED 3 is interesting because, in the 
absence of fermion loop contributions to the photon polarisation tensor; i.e., in quenched approximation, 
the lattice version of the theory has been shown to be confining (Gopfert and Mack, 1982). This result 
can be seen heuristically by considering the classical potential: 

V(x) = jT dra / -0-e^ + ^e 2 D( q ) = J -0^ e ^e 2 D(q) . (3.4) 

Neglecting the vacuum polarisation; i.e., setting II(g 2 ) = 0, one obtains 

e 2 / x 

V(r) = —In (re 2 ) (3.5) 

which exhibits logarithmic confinement. 
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It is also of interest to those concerned with the dynamical generation of mass in grand unified theories. 
This is because the coupling in QED 3 , e 2 , has the dimensions of mass. In addition, the theory is 
super-renormalisable which means that there are no ultraviolet divergences whose regularisation would 
introduce a new scale that could complicate the relation between the scale of dynamical mass generation 
and the natural length scale provided by the coupling, e 2 . 

Another application is the study of field theory at finite temperature . It has been argued that the high 
temperature behaviour of a given four-dimensional gauge theory, which has a small coupling, at distances 
greater than some screening length is precisely described by the corresponding three-dimensional theory 
(Appelquist and Pisarski, 1981). It is worth noting that QED 3 at finite temperature, with fermions 
coupling to two Abelian gauge fields [i.e., the usual electromagnetic field plus a "statistical" gauge field] 
has been applied to a description of high T c superconductivity in the quasi-planar oxides La2Cu04 and 
YBa2Cu6 (Dorey and Mavromatos, 1992). 

In every one of these cases it is the dynamical generation of mass; i.e., the dynamical breaking of 
chiral symmetry, and its relation to, and effect on, the gauge boson propagator that is of interest. 
These quantities are most easily studied within the DSE approach. Indeed, an order parameter for 
DCSB is the fermion condensate, (if)i[>) , which is obtained from the fermion propagator, Sp(x), via: 
(ipip) = —tTSp(x = 0), when both the Lagrangian, m , and renormalised, m, bare masses are zero. 



£(p) = i 1 . v {^Z 1 -l)-{m-Z l mQ) + Z l e 2 J -^^Y^(p,q)D^{p-q]i)S F {q)-i v , (3.6) 



The DSE for the renormalised fermion self energy, which is discussed in Sec. |2.1| , is therefore a natural 
equation to study in order to address the issue of DCSB and in Euclidean space it can be written as: 

d 3 q 

where the Ward Identity Z\ = Z2 has been used, and in Euclidean space the general form of the solution 
can be written in the form 

S(j,) = -ry.po v {p)+a 8 (p) = - .7 • P p2A{p2)2 + B(p2)2 + p2A(p2)2 + B(p2)2 ■ (3-7) 

[Note that we are no longer using a tilde to denote these renormalised quantities.] In the absence of 
bare mass a non zero value of (ipip) signals DCSB and a plot relating (ipip) to relevant dimensionless 
parameters in a given model can be used to study the transition to the chiral symmetry breaking phase. 

In d = 3 dimensions there are two inequivalent 2x2 representations of the algebra {7^,7^} = 25^. 
Hence, to describe spinorial representations of the Lorentz group, two component spinors are sufficient. 
However, in this case any mass term, whether explicit or dynamically generated, has the undesirable 
property that it is odd under parity transformations. This can be avoided if one employs four component 
spinors and a 4 x 4 representation of the Euclidean Dirac algebra; for example, the set 71, 72, 74. It 
is clear that since both 75 = —71 72 74 and 73 anticommute with this set then the massless theory is 
invariant under two transformations: 

if; -> e iQri *lj and i/j -> e iai5 ijj (3.8) 

which are analogous to the chiral transformations in d = 4 dimensions. In this case there are two types 
of mass term 

miipip and m, 2 V4[73, 75]^ • (3.9) 



The first of these is invariant under parity transformations but not under Eqs. ( |3.8| ) while the second 
is invariant under the transformations of Eqs. (|3.8|) but not under parity. The mi term is clearly the 
d = 3 analogue of the mass term in four- dimensions and it is common practice to use only this type of 
mass term in QED3 (Pisarski, 1984). 
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3.1 Quenched Approximation 



In massless QED 3 the quenched approximation, which corresponds to neglecting fermion loop contri- 
butions to the vacuum polarisation; i.e., to setting 



I% 2 ) EE 



(3.10) 



in Eq. ( |3.3|) , leads to infrared divergences in ordinary perturbation theory based on an expansion in the 
coupling, e 2 . This divergence is evident in the lowest order vertex correction: 



A p (k,p) ~ e 2 J 



d 3 q 1 



5„ 



(2tt) 3 q 2 V 
In the limit p — > 0, q — ■> this is 



(1-0^ 7, 



1 



ij ■ (p — q)^ P i^i ■ (k — q) 



7* 



oc 



(«-»)/ 



d 3 q 1 
(2tt) 3 7 



(3.11) 



(3.12) 



which has a manifest infrared divergence. [In renormalising QED 3 this divergence is incorporated into 
the vertex renormalisation constant, Z 1; which is described in Sec. |2.1| . For d = 4 the gauge dependent 
factor is £ — 1 and hence Z\ — \'va. Landau gauge at 0(e 2 ) in Abelian theories.] 

One commonly used remedy for this problem is to soften the infrared behaviour of the gauge boson 
propagator by including the fermion loop contribution to the vacuum polarisation (Appelquist et ai, 
1986). For a fermion of mass m the lowest order [in e 2 ] contribution is given by 



d 3 q 



tr 



7/i 



■ ' 7i> ~ 

i'J ■ (q + \k) +m v i-f ■ (q — \k) + m 



(3.13) 



Evaluating this using a gauge invariant regularisation scheme, such as dimensional regularisation, one 
obtains [recall Eq. ( ggOfl ] 



S. 



fJtU 



kfj, k v 



U(k 2 



where the polarisation scalar is 

U(k) = 



k 2 — m 2 
2m H arcsin 



k 



y/k 2 + 4m 2 



(3.14) 



(3.15) 



In a theory with iV massless fermions the polarisation scalar is therefore 

TL(k) 



a 
k 



where a = Nfe 2 /8, in which case the gauge boson propagator is given by Eq. ( |3.2j ) with 

D{q) 



q 2 + aq 



(3.16) 



(3.17) 



which behaves as 1/q for q 2 0; i.e., the infrared divergence is softened significantly without altering 
the appealing ultraviolet properties. 

Allowing such a contribution from massless fermions introduces a dimensionless parameter, Nf, the 
number of fermions, into the fermion DSE. This breaks the scale invariance of this equation and makes 
possible the existence of a phase transition associated with DCSB; i.e., a situation where (ipip) is nonzero 
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for some values of Nf and zero otherwise. This is only possible when there is a dimensionless parameter 
in the equation. [There have been a number of studies of this issue and we will discuss them below.] 



However, the price of allowing such a contribution from massless fermions is to eliminate confinement. 
It is a simple matter to calculate the classical potential, Eq. (|3.4| ), in this case and one obtains (Burden 
and Roberts, 1991) 

V(r) = -j [H (ar) - N (ar)} , (3.18) 

where H (x) is a Struve Function and N (x) a Neumann Function (Gradshteyn and Ryzhik, 1980, Sees. 
8.4, 8.5). From Eq. (|3~T8|) one finds easily that 

e 2 1 

V(r) ~ ln(ar) at small r and V(r) ~ ^— at large r (3.19) 

2n ar 

and thus observes that the theory is no longer confining. This deconfinement is akin to that which 
occurs in lattice simulations with dynamical fermions: the fixed source charges dress themselves with a 
cloud of massless, charged fermions which strongly screen the source charge. 



3.2 Rainbow Approximation 



In studying DCSB using Eq. ( |3.6|) a commonly used approximation is to write 

rV(p,g) = 7„ (3.20) 

which is referred to as the "rainbow" or "ladder" approximation [it is the DSE analogue of the lad- 
der approximation in the BSE]. If this approximation is combined with the quenched approximation, 
Eq. ( |3.10p , then Eq. ( |3.6| ) decouples from the remaining DSEs and it can be studied in isolation. This 



is the merit of this approximation, however, there are significant problems associated with rainbow 
approximation; notably the loss of gauge covariance and a solution for the fermion propagator which 
has an unphysical singularity structure. 

A commonly used simplification of Eq. ( |3.6| ) is to set Z\ = 1. As remarked in connection with Eq. ( |3 . 1 2|) , 
this is true in QED4 at 0(e 2 ) but in QED3 it is just a convenient simplifying truncation whose validity 
can only be justified a posteriori. [The proponents of a 1/Nf expansion in QED 3 , where Nf is the 
number of light fermions [see Eq. (|3.13|) and the associated discussion] may argue that Z 1 = 1 and 



rainbow approximation form a consistent pair of approximations valid at leading order in 1/Nf .] 

In quenched, massless, rainbow approximation QED3 one finds easily that 

A(p 2 ) ee 1 (3.2i; 
in Landau gauge, £ = in Eq. ( |3.2|) , and obtains the following nonlinear integral equation for B(p 2 ): 

p + q 



9 -i- f r°° 



p-q 



(3.22) 



q 2 + B 2( q 2\ 



In other gauges A(p 2 ; £ 7^ 0) quickly approaches 1 as p 2 increases and it is a good approximation to use 
A(p 2 ) = 1 in most studies of quenched, massless, rainbow approximation QED 3 . 
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m = + - (3-24) 



An operator product expansion in QED3 suggests that when propagating in the presence of a condensate 
(ipip) 7^ the fermion propagator will receive a self mass contribution of the form (Burden and Roberts, 
1991): 

- 2tt 3 (^0) 5 3 (q) . (3.23) 

The factor 5 3 (q) is present because, by definition, the condensate does not exchange momentum with 
the fermion and the numerical factors are simply to ensure appropriate normalisation. Including this 
term as a perturbative contribution to the fermion propagator one finds 

ry • q 4 q 4 
and the OPE analysis then predicts that as p 2 — > 00 

This provides a means of checking, in the large-p 2 domain, the approximations and solution procedure 
used to analyse the QED 3 fermion DSE. 

The large-p 2 asymptotic behaviour of the solution of Eq. (|3.22|) can be obtained by analysing an ap- 
proximate differential equation that is valid in this domain. Approximating the kernel as follows: 

In p + q 
p-q 

which is a good approximation for p 2 <C q 2 or p 2 ^> q 2 (Roberts and McKellar, 1990), one obtains the 
following differential equation: 



-e(p-q) + -0(q- P ), (3.26) 
p q 



^-(p 3 ^-B(p)) + 2 -^B(p) = 0. (3.27) 
dp \ dp J ty z 



A differential equation valid over a greater p 2 -domain could be obtained using the methods discussed 
by Munczek and McKay (1990), however, this is not necessary for a comparison with Eq. ( |3.25p . 

The solution of Eq. ( |3.27|) that is consistent with the ultraviolet boundary condition B(p 2 ) — > as 
p 2 — > 00 is 




b( P ) = K-j 2 \r-^) - ( 3 - 28 ) 



where Jg(x) is a Bessel function of integer order and k is a constant that cannot be determined by the 
differential equation. It can, however, be determined by comparing Eq. (|3.28 ) with Eq. ( |3.25|) : 



7T 2 ,- 



K = -—(^) . (3.29) 



Burden and Roberts (1991) tested these predictions [Eq. ( |3.25| ) c.f. Eqs. ( p.28| ) and ( |3.29| )1 with the 
solution of Eq. (|3.6| ) in quenched, massless, rainbow approximation and with Zi = 1. In these studies 
A(p 2 ) was not neglected and the coupled, nonlinear integral equations for A(p 2 ) and B(p 2 ) were solved 
by iteration. In QED 3 , with a 4 x 4 representation of the Dirac Matrices, 



7T Z 



dp ■ 



p 2 B(p 2 



(3.30) 



p 2 A 2 (p 2 ) + B 2 (p 2 ) ' 

which is a convergent integral. As will be seen in Table |3]2| there is excellent agreement between the 
predictions and the numerical results. This indicates that at large spacelike p 2 the quenched, massless, 
rainbow approximation DSE with Z\ — 1 is a good approximation. 
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Table 3.1. Comparison of the asymptotic form of the fermion self mass 
with (ipip) to test the OPE prediction. [Adapted from Burden and 



Roberts 



(1991).] 



a 


P 

units of e 2 


units of 10 -3 e 4 


units of 10 -3 e 4 


-0.2 


511 


2.638 


2.638 




1000 


2.638 




0.5 


511 


1.775 


1.775 




1000 


1.775 




1.2 


511 


1.352 


1.352 




1000 


1.352 





3.3 Beyond rainbow approximation 



The loss of gauge covariance when using the rainbow approximation is a direct consequence of the 
violation of the Ward Identity: 

(p-q)^,(p,q) = S- 1 (p)-S- 1 (q) . (3.31) 

It is clear that if the solution of the fermion DSE, Eq. ( |3.6|) , is momentum dependent, as it will be if 
the kernel is, then this identity is not satisfied by the bare vertex. The correct form of the fermion- 
gauge-boson vertex is therefore crucial in restoring gauge covariance to the DSE. 

The structure of this vertex has been analysed to order a 2 in Abelian theories by Ball and Chiu (1980). 
One important result of this study is that the dressed vertex should be free of kinematic singularities; 
i.e., that T^p, q) should have a well defined limit as p 2 — > q 2 . With this in mind, quenched QED3 was 
studied by Burden and Roberts (1991) with the following "light-cone regular" Ansatz for the fermion- 
gauge-boson vertex: 



«IV(p,g) 



+ 



aA{p 2 ) + (1 - a)A{q 2 

(P + g)/x 

p 2 — q 2 



{t [A(p 2 ) - A(q 2 )] [(1 - a) 7 ■ p + a 7 • q] + [B(p 2 ) - B(q 2 )] } . (3.32) 



This vertex is a simple modification of that proposed by Ball and Chiu (1980). 

It is important to note that the constraint that the vertex be regular as p 2 — > q 2 entails that it cannot 
be purely longitudinal and hence that the solution of the DSE will be sensitive to the Ansatz even in 
Landau gauge; £ = in Eq. (|3.2[) . The parameter a in Eq. (|3.32|) was included in order to allow, in a 
simple way, for a variation of the transverse part of the vertex so that its effect on the solution could 
be studied. 



Burden and Roberts (1991) studied Eq. (|3.6| ) in the quenched, massless limit, with the vertex Ansatz 
of Eq. ( 3 . 3 2|) and with the additional approximation of setting Z\ = 1. In this case the fermion DSE 
reduces to a pair of coupled, one-dimensional integral equations: 



P' 



(A( P 2 



- 1 



q 2 dq 



4tt 2 Jo q 2 A(q 2 ) 2 + B(q 2 )- 



■x 
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and 



p2 _ g2 



{q 2 A(q 2 )-p 2 A(p 2 )}A(q 2 ) - {B(p 2 ) - B(q 2 )} B(q 2 



H In 

2pq 



p + q 



p-q 



p2 _ q2 



2 'i-t+tm 



{q 2 A(<?) + p 2 /t(p 2 )} A(q 2 ) - {B(p 2 ) - B( 9 2 )} B(q> 
p + q \ 



2pq 



p-q 



1 _ a y + aq 2 } { A(p2) _ A (q 2 )}A(q 2 ) + {B(p 2 ) - B(q 2 )} B(q 2 



B{p 2 



1 



q 2 dq 



An 2 Jo q 2 A(q 2 ) 2 + B(q 2 ) 2 



p2 _ g2 



H In 

2pq 



pq 



p + q 



p-q 



{A{p 2 )B{q 2 ) - A{q 2 )B{p 2 )} 



p + q 



p-q 



4(1 - a)A{q 2 )B{q 2 ) + (4a + ti)A{p 2 )B{q 2 ) + ^A{q 2 )B{p 2 



(3.33) 



(3.34) 



These equations are scale invariant. The mass scale is set by [i = e 2 and the solution for any value of 
H can be obtained from the fi = 1 solution by the scale transformation: 



A(p 2 ; f ,) = A(^;l 



and 



B(p 2 ]f i)=^B(^;l) . 



(3.35) 



It suffices therefore to solve the equations for fi = 1 and hence this choice has been made in Eqs. ( |3.33| ) 
and (|3.34|) . As remarked above, this feature of scale invariance means that there can be no critical 
coupling parameter in quenched, massless QED3 because there are no dimensionless parameters: once 
a chiral symmetry breaking solution exists for one value of e 2 it exists for all values of e 2 . 



The most commonly used procedure for solving this type of coupled, nonlinear integral equations is 
simple iteration: a guess is made for A(p 2 ) and B{p 2 ) and substituted; the result of the integration is 
then resubstituted and the procedure repeated until the functions in the integrand reproduce themselves. 
This is a practical and efficacious approach. There are analytic methods that can be used to establish 
the uniqueness of the solutions of certain nonlinear integral equations (McDaniel et al, 1972) but 
hitherto no attempt has been made to employ them in the DSE approach. The numerical stability of 
the solutions has been used as an implicit indication of uniqueness. 



This procedure can be used to solve Eqs. ( |3.33| ) and ( |3.34j ) with the result that for all values of the 



parameter a in Eq. ( 3.32 



A(p 2 ) ^ 1 even for f = . (3.36) 

This is a direct result of the transverse piece in Eq. ( |3.32| ); i.e., of the requirement that a vertex which 
satisfies the Ward identity be regular in the limit p 2 — > q 2 . [A study of a "light-cone singular" vertex 
that has been used by a number of authors confirms this (Burden and Roberts, 1991).] The vertex 
of Eq. ( |3.32|) also has the highly desirable feature that it leads to an [almost] gauge invariant fermion 
condensate. This is illustrated in Fig. 3.1. It will be observed that on the domain < £ < 1 there is a 
value of 

a w 0.53 (3.37) 

for which the explicit gauge dependence of the photon propagator is compensated by the implicit gauge 
dependence of the solution functions. This indicates that the transverse parts of the vertex are extremely 
important in restoring gauge independence to {ipip) . 
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Fig. 3.1. Plot of (ipip) , obtained from the solutions of Eqs. ( |3.33| ) and 
(|3.34|) as a function of the gauge parameter £ for various values of a, 
the parameter that varies the transverse part of the vertex. A value of 
a W 0.53 minimises the sensitivity of {ipip) to £ (Burden and Roberts, 
1991). 



3.4 Comparison with Lattice Simulations 



The approach described in Sec. |3]^ leads to an [almost] gauge invariant condensate which invites com- 
parison with (ipip) obtained in lattice simulations of QED3 . In the lattice formulation 

= {tlM y lU]) (3.38) 

where (■) represents a Boltzmann weighted average over the gauge field U, M~ 1 [U] is the inverse of the 
lattice Dirac operator in a given gauge field configuration and V is the number of lattice sites. Employing 
the formalism developed by Burden and Burkitt (1987) one obtains the following relationship between 
the lattice and continuum condensates: 



-<^> = rftW) (3-39) 

2 Or 

with Nf the number of fermion flavours which is one in the case we are considering here. 



For a proper comparison one must look to the lattice results obtained in a scaling window where they 
are thought to best represent continuum physics. One signal that a scaling window exists in the lattice 
theory is the observation that 



,2 — v lattice 



P*{>H>) = K. (a constant) (3.40) 



for large (3 = -r— . In the scaling window it should be that 



e 2 a 



1 ,-r,. continuum Nt „ ,„ 
— <^> = -fK.. (3.41) 
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0.015 



0.010 



0.005 




0.5 



1.5 



(3 = l/(e 2 a) 



Fig. 3.2. The chiral condensate from the quenched lattice studies of 
Dagotto et al. (1990) for 8 3 [diamonds] and 10 3 [squares] lattices. The 
solid line is (ipip) obtained in the DSE approach to quenched QED 3 . The 
vertical axis is —{ipip) /e 4 . This figure is adapted from Burden (1993). 

The simulations of Dagotto et al. (1990) suggest that in QED 3 a scaling window exists for (3 > 1. 



The calculations of Burden and Roberts (1991) correspond to the quenched approximation of the lattice 
formulation. This corresponds to the Nf = simulations of Dagotto et al. (1990). There exists now the 
possibility for confusion because of the Nf factor in (|3.41 

that the factor (detM)a ; present in Nf = 1 simulations, is not included in the Nf 



The quenched approximation actually means 

simulations. In 



this case the steps that led to ( p. 39 ) can be retraced and one finds that the proper comparison is 



continuum ^ l^ /=0 



(3.42) 



This comparison is presented in Fig. 3.2. It will be noted that, within the scaling window, the agree- 
ment with the 10 3 simulations of Dagotto et al. (1990) is very good. Indeed, there is little room for 
improvement which suggests that, in the quenched approximation, the Ball-Chiu vertex, Eq. ( p.32| ) with 
a = 0.5, is a very good approximation to the true vertex; or at least to that part of it which contributes 
to the fermion DSE. 



Computationally, the DSE calculations are much simpler and quicker than the lattice studies and Fig. 3.2 
illustrates that the DSE approach to solving Abelian field theories is becoming competitive with lattice 
studies. 



3.5 1/Nf Expansion and a Phase Transition in QED3 



As we mentioned in association with Eq. ( |3.17|) , including fermion loops in the photon polarisation 
tensor introduces a dimensionless parameter, Nf, the number of flavours of electrically active fermions, 
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into QED3 upon which (ipip) may depend. Indeed, there may be some value of Nf = N CT for which 
(ipip) — > 0. Such behaviour has been inferred from studies of non-compact QED 3 on 8 3 (Dagotto et 
al, 1989) and 10 3 lattices (Dagotto et al, 1990) with (0-0) ^ only for N f < N cv = 3.5 ± 0.5. The 
existence of a value of Nf = N CT above which there is no DCSB would limit the utility of QED 3 as a 
model for strongly interacting systems. 

The lattice studies referred to above followed a DSE study of QED 3 at leading order in an expansion in 
l/Nf (Appelquist et al, 1988). In this study, which used Landau gauge, it was argued that at lowest 
order in l/Nf, A(p 2 ) = 1; i.e., that corrections to A arise at higher order in the l/Nf expansion, and 
that B is given by the solution of 



B(p) 



kB{k) 
N f n 2 p Jo k 2 + B(k) 



dk ■ 



In 



k + p + a 



\k 



P\ 



a 



(3.43) 



We remark that Eq. fl3.43| ) incorporates the effect of massless fermion loops in the photon polarisation 
tensor; i.e., it goes beyond quenched approximation at leading order in l/Nf-. see Eq. ( |3.15| ) and the 
associated discussion. 



This equation could, of course, be solved directly by iteration, however, Appelquist et al. (1988) argued 
that the In term rapidly suppresses the integrand for p > at and that 



B(p) 



dk — ^ttttt (k + p-\k-p\) 



Nfir 2 pJo k 2 + B(k) 2 
is a good approximation. Equation (|3.44j) is equivalent to the differential equation 



d 

dp 



P 



dB{p) 
dp 



+ 



with the boundary conditions 

< B(0) < 00 and 



8 p 2 B(p) 
N f n 2 p 2 + B(p) 2 

dB(p) 







P- 



dp 



+ B(p) 



(3.44) 



(3.45) 



(3.46) 



If one assumes that in the domain < p < a there is a subdomain p\ < p < a on which B(p) <p then 
Eq. ( |3.45| ) can be linearised; i.e., the B 2 term in the denominator neglected, and the solution in this 
domain is 



B(p) oc p a where 




32 
Nfir 2 



(3.47) 



For Nf > 32/tt 2 one has — 1 < a < and it is not possible to satisfy the ultraviolet boundary condition 
at p = a in Eq. ( |3.46| ). However, for Nf < 32 /it 2 the solution has the form 



B(p) = — sin. 

VP \ 



32 
Nf^ 



In 



P 



B(0) 



+ 5 



(3.48) 



where 5 is some phase, which may depend on Nf, and the logarithm has been scaled by -B(O). Any 
dimensioned parameter would have done but all are proportional to B(0) so this is completely general. 
Imposing the ultraviolet boundary condition one finds that, for Nf —* 32 /ir 2 , 



B(0) 



~ 2+8 
ae exp 



-2nn 
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(3.49) 
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Since (ipip) vanishes if B(0) — > 0, the conclusion of this analysis is that in QED3 there is a critical 
number of flavours of massless fermions 

32 

iV cr = - (3.50) 

7T Z 

below which there is DCSB and above which chiral symmetry is restored. [This assumes that 5 is not 
singular as Nf — > iV cr .] Equation ( |3.43j ) was also studied numerically by Appelquist et al. (1988) and 



- 4M = 1 +4 37v^ 1 HIJ + ]v;' (3 - 51) 



the approximations associated with deriving the differential equation and its solution were found to be 
valid. The same result has been found by Maris (1993, Sees. 5.1 and 5.2) when A(p 2 ) = 1 is imposed. 

Higher order corrections in 1/Nf, 0(1/Nj), have been calculated Nash (1989). In this study it was 
argued that 

2 ~Ji] (E\ _£. 

where c is a constant, and an integral equation was derived for M(p) = B(p)/A(p) in which the kernel 
was simplified; the critical number of flavours was independent of the gauge parameter. The solution 
of the approximate equation for M(p) is of the form 

M(p) =p a (3.52) 

where a is the solution of 

32 / 7.81 \ 

The ultraviolet boundary condition of Eq. ( |3.46|) allows a non trivial solution for M(p) if a becomes 
complex. Neglecting the 0(1/ Nf) terms in Eq. ( |3.53| ) then this condition yields N CT = |32/7r 2 « 4.32. 
[The discrepancy between this result and that of Eq. ( |3.50|) arises because Nash (1989) took account of 
vertex corrections.] Solving the full equation leads to 

N cr « 3.28 (3.54) 

which is a 25% decrease from the 0(1/ Nf) result. [The solution Nf « 1.04 was discarded by Nash 
(1989) as being "unphysical".] This shift in N cv might be taken as evidence that the 1/Nf expansion 
converges reasonably well. 

The primary conclusion of this study is that the 0(1/ Nj) terms do not qualitatively change the nature 
of the solution: there is still a critical number of massless fermion flavours above which chiral symmetry 
is restored. 



Dynamical Chiral Symmetry Breaking for all Nf. The existence of a critical number of massless fermi- 
ons in QED 3 is not universally accepted. Pennington and Webb (1988) and Atkinson et al. (1988d) 
studied the pair of coupled, nonlinear integral equations for A(p) and B(p) obtained in QED 3 including 
the vacuum polarisation of Eq. ( |3.16| ) but working in rainbow approximation, Eq. ( |3.20j ). These studies 
suggest that chiral symmetry is dynamically broken for arbitrarily large Nf. 



The origin of this difference is the behaviour of A(p) for small p 2 . Indeed, it has been shown (Pennington 
and Webb, 1988) that in perturbation theory, when B(p 2 ) is replaced by its value at p = 0, m = B(0), 
in order to approximate its effect on A(p): 



A(p) = 1 - 



3N f ir 2 



a 2 



m 2 + p 2 



4 m 
3 p z 



m ( p 

— arctan — 
p \m 



1 



(3.55) 
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Considering Eq. ( |3.55|) in the limit p —>■ one finds 



m 



A ^) = l + ^T-2 U + la- . (3-56) 



c./. Eq. ( p.51| ) in Landau gauge, and with 



ln m 7^=j (3 ' 57) 



near the supposed critical point, see Eq. ( |3.49|) , then 



A (p) ; (3.58) 

N f JN CT - N f 



i.e., it is singular at the critical point. One observes then that for momenta such that m <C p < a, 
A(p) does indeed receive small corrections at "large Nf" . However, for small momenta, p <C m, which 
is the domain relevant to DCSB, 

A(p) ± 1 + 0(1/N f ) (3.59) 
and the 1/Nf expansion is therefore not a valid tool for analysing DCSB. 

There is further evidence that supports the assertion that DCSB in QED 3 cannot be properly analysed 
without considering A(p) in detail. A study of the Cornwall- Jackiw-Tomboulis (1974) [CJT] effective 
action in QED3 using the 1/Nf expansion reveals that this effective potential is always unstable against 
fluctuations away from A(p) = 1 and B(p) = (Matsuki, 1991). This result suggests the chiral 
symmetry is dynamically broken in QED 3 for all Nf. 

Going beyond perturbative arguments, the unquenched, rainbow approximation [Eqs. (|3.16|) and ( |3.20|) , 
respectively] DSE in QED3 can be solved numerically. This equation yields the following pair of coupled 
equations: 



N f 7T 2 p 3 Jo k 2 + B(k) 2 
\ ln ( jTr-Vrj ) +i(k + P~\k- p\) (1 + \p 2 - k 2 \) - kp 



\{k 2 -p 2 ) 2 lnf #^±^ l (3.60) 



(k+p)(\k-p\ + 1) 



\k — p\ + 



(3.61) 



where the momenta have been rescaled; e.g., p — > a p. Analysing the numerical solution (Maris, 1993, 
Sec. 5.3) one observes that the result of Eq. (|3.58|) does not survive but that instead A(p) ~ 1/^/iV/. 
This is, nevertheless, not of the form l+0(l/iV) and again suggests that the 1/iV/ expansion is not 
valid. Solving these coupled equations numerically yields the result that chiral symmetry is dynamically 
broken for all Nf. 



There have been a number of other studies of this problem which go beyond the rainbow approximation 
(Atkinson et at, 1990; Walsh, 1990). The approach in these studies is to make an Ansatz for the 
fermion-gauge-boson vertex, which satisfies certain reasonable constraints, and solve the DSE that 
results. Perhaps the most sophisticated of these studies is that of Curtis et al. (1992) who have solved 
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the coupled integral equations for A(p) and B(p) obtained using the following Ansatz for the fermion- 
gauge boson vertex which satisfies the WTI and ensures multiplicative renormalisability of the fermion 
DSE in QED 4 : 



where 



r M (M = rf(M + rf(M 

A(p 2 ) + A(k 2 " 



(P + k) t 



A(p 2 ) - A{k 2 ) 



(3.63) 



(3.64) 



[«7 • p + ry • k] 



p 2 — k 2 

is the vertex proposed by Ball and Chiu (1980) and 

lii{k 2 - P 2 ) - (k + p)n(h -k-ij -p) 



+ 



B{p 2 ) - B(k 2 



^ P (k,p) 



with d(k,p) 



(k 2 + p 2 ) 



2d(k,p) 
(k 2 -p 2 ) 2 + 



A{k 2 ) - A(p 2 ) 



B 2 {k 2 ) I 'rij, 2 ] 1 
A 2 (k 2 ) + A 2 (p 2 ) 



(3.65) 
(3.66) 



is the correction term introduced by Curtis and Pennington (1990) in order to restore multiplicative 
renormalisability in studies of the DSE in QED4 . 

This study encountered numerical problems for Nf > 7.2, simply due to the large momentum domain 
on which the solution must be known very accurately at large Nf. Nevertheless, the results obtained 
suggested that 

~ -i 

A(0) 



In ^M(O) = ) N f ; (3.67) 

i.e., that M(0) becomes exponentially small as Nf is increased but that DCSB persists for all Nf. 



We remark that this sort of behaviour of M(p) would be very difficult to extract from lattice simulations 
of QED 3 and, if real, would cast doubts upon the critical coupling inferred in those studies carried out 
to date (Dagotto et al, 1989, 1990). It should also be noted that the finite size of a lattice unavoidably 
leads to a suppression of M(p) for Nf > 1 (Kondo and Nakatani, 1992a) which may lead to an artificial 
phase transition. 

There is a criticism of the vertex Ansatz approach, however, and it is that it is not obvious how to 
systematically improve over a given Ansatz. The DSE for the vertex itself is a complicated equation 
and no attempt has as yet been made to solve it directly, which would be the natural, systematic step 
in the DSE approach. Failing this, however, even applying certain, very general constraints to restrict 
the form of the vertex, which we discuss in Sec. |3[7], the Ansatze remain ad hoc. As a consequence 
there are those who doubt the result that QED3 has dynamical chiral symmetry breaking for arbitrary 
numbers of massless fermions. Lattice simulations of QED3 with greater accuracy might therefore be 
valuable. 



3.6 Photon Polarisation Tensor and Mass Generation 

Another weakness of many of the studies of DCSB in QED3 undertaken to date is that they do not take 
into account the effect of the dynamically generated fermion mass on the photon polarisation tensor. 
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In Abelian gauge theories, once an Ansatz for the fermion-gauge-boson vertex is made then the DSEs 
for the fermion propagator and the photon polarisation tensor form a closed system of equations. 



The DSE for the photon polarisation tensor in Euclidean momentum-space is given in Eq. ( |3.13 ) and it 



is a simple matter to determine whether the dynamical fermion mass has any effect on the polarisation 
tensor. At first order one must simply evaluate the polarisation tensor using a solution of the fermion 
DSE evaluated with the same vertex Ansatz but with H(k) = 0. This is the first step in a self consistent 
solution of the coupled fermion and photon equations and has been carried out by Burden et al. (1992a) 
using the vertex of Eq. ( |3.32| ) with a = 0.5; i.e., the Ball-Chiu vertex. Such a calculation also allows for 



another simple test of gauge covariance in the DSE approach since the polarisation scalar, see Eq. ( |3.14|) , 
should be independent of the gauge parameter; i.e., 

-I n(fc; £) = . (3.68) 



In evaluating the polarisation tensor one must first note that with the vertex Ansatz of Eq. ( |3.32|) 



A(k 



oo 



1 



16Jfe 



and Bik — > oo) = — 



A 

£2 



(3.69) 



with A a gauge parameter dependent constant (Burden and Roberts, 1991), and hence Eq. ( |3.13|) must 
be regularised since a direct evaluation yields k^H^ v {k) as a difference of two logarithmically divergent 
integrals in all but Landau gauge. This problem can be avoided by contracting the polarisation tensor 
with 

*v = ( 3 - 7 °) 



k 2 



which is orthogonal to 5^ u in d dimensions (Walsh, 1990, Chap. 3). Since the divergent part of the 
integral is proportional to 5^ then, after contraction, only the finite piece of the polarisation scalar 
remains: 



U(k) 



d 3 q 

2k 2 J (2^)3 



tr J^S{q + \k)Y^q + \k,q- \k)S(q - \k) 
Z^Siq + \k) k v T u (q + \k,q- \k) S(q - \k) 



(3.71) 



At this point it is a simple matter to address the question of photon mass generation since there are a 
number of lattice studies of the gauge boson propagator in QED 4 (Coddington, et al, 1987; Mandula 
and Ogilvie, 1987; R. Gupta et al., 1987; Mandula and Ogilvie, 1988) which find a nonzero gauge boson 
mass after fixing a lattice Landau gauge. To study this one must consider 11(0) since a 1/p 2 pole in 
the polarisation scalar signals mass generation via the Schwinger mechanism (Schwinger, 1962). From 
Eq. (|3.16p it is clear that there is no such mass generation in perturbation theory. 



Writing 



k 2 U(k) = -N f -f(k) 



and using the differential form of the Ward identity [r At (g, q) = —id q S 1 {q)} one finds 



(3.72) 



f(0) = lim — / 



d 3 q 



tr 



7m 



7 • kk. 



diS(q) 



(3.73) 
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Fig. 3.3. The polarisation scalar, H(k), obtained from Eq. ( |3.71|) using 
the solution of Eqs. ( |3.33j ) and ( |3.34j ) with a = 0.5, is plotted for £ = 0: 
solid line; £ = 0.5; dashed line; £ = 1.0; dash-dot line. [The mass scale is 
fixed by setting e 2 = 1 .] The plot illustrates the deviation of the result 
from the perturbative one at small k. On the scale of this plot the fitting 
function of Eq. ( |3.75| ) lies exactly on top of the numerical results. The 
perturbative result of Eq. ( |3.16| ) is the short-dash line. [This figure is 
adapted from Burden et al. (1992a)] 



where the final limit must be taken with care after the Dirac trace is evaluated. Using Green's theorem: 

/(0) = lim — [ d 2 S ( 1 - 3 { 4^) 9A9 A } q) „ n ( 3 - 74 ) 
and since Jq d9 sin #[1 — 3 cos 2 6] = one has /(0) = 0. 

Hence, as long as the Ward identity is satisfied the photon remains massless independent of the gauge 
parameter and details of the transverse part of the vertex. This is a general result in Abelian theories, 
independent of the details of the interaction and, suitably modified, it also holds in QED4 . One 
expects therefore, that in the continuum limit lattice simulations should also yield this result. 

The form of U.(k) obtained by Burden et al. (1992a) for a range of values of £ is plotted in Fig. 3.3. 
The results can be summarised by the fit: 

n(Jb) = r + he~ ck i (3.75) 

8(P + e 4 a 2 )5 

with £ dependent fitting parameters a, b and c specified below, (£, a, b, c): 

(0.0, 0.2044, 8.760 x 10~ 2 , 7.767) , (3.76) 
(0.5, 0.2541, 5.055 x 10~ 2 , 11.70) , (3.77) 
(1.0, 0.2778, -9.349 x 10" 3 , 0.5285) . (3.78) 
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[Of course, in this calculation, a trivial multiplicative factor of Nf appears on the right hand side of 
Eq. fl3.75| ) if one uses Nf fermions. A complicated dependence on Nf can only arise when one carries 



out the simultaneous solution of Eqs. ( |3.33|) , ( |3.34|) and (|3.71|) . 



The first thing one deduces from Eq. ( |3.75|) is that the Ball-Chiu vertex, although yielding (ipip) essen- 



tially independent of the gauge parameter, £, it does entail Eq. ( |3.68D ; i.e., it does not lead to a gauge 



independent polarisation scalar. This is important and we will discuss it in some detail below in Sec. RT7 



One observes from Eq. ( 3.75|) and Fig. 3.3 that, at large k 2 , the nonperturbative result returns to the 



perturbative one of Eq. (|3.16| ). For small k 2 , H(k) differs from the perturbative result which was used 
in the analyses described in the preceding subsections and predominantly in Pennington and Walsh 
(1991). The difference is a quantitative softening of the infrared behaviour of IT(/c) and it has been 
reported (Pennington and Walsh, 1991) that this has little effect on DCSB when used simply as input 
to Eqs. ( |3.33| ) and ( |3.34| ). The difference is not unimportant, however. Simply assuming that II(/c) is 
bounded [in absolute value] and continuously differentiable on (0, oo) and that II(A;) ~ 1/k for k — ► oo 
is enough to ensure that, for large r, the potential associated with the dressed photon propagator is 



e 2 



V(r) = — —r — In e 2 r + const. + h(r) (3.79) 

1 + n(0) 2n 

where h(r) falls off at least as fast as 1/r (Burden et ai, 1992a). Hence, the large r behaviour of the 
potential is dominated by a confining logarithmic term. It will be recalled that this was true of the 
bare photon propagator but not of the perturbatively dressed photon propagator. The inclusion of a 
dressed fermion-photon vertex and dressed fermion propagator [obtained with this vertex] at this level 
has thus restored the confining nature of QED 3 . 

The first order response of the photon polarisation scalar to dynamical fermion mass generation is a 
modification on the small-/c 2 domain which is just that domain relevant to DCSB. This may have an 
impact on the question of the existence of a critical number of flavours in QED3 but a full analysis of 
the coupled fermion-photon DSE system with a physically reasonable vertex Ansatz would be the least 
necessary to address this. Preliminary results of such a project are described by Walsh (1990) but no 
firm conclusion may be drawn from this study. 



3.7 Fermion-Photon Vertex and Gauge Invariance 

The issue of gauge invariance in the DSE approach to the solution of gauge field theories is very 
important. Indeed, a manifest lack of gauge covariance in many works in this field has hindered its 
acceptance as an efficacious nonperturbative tool. [There are many contemporary works that still fail 
to address this issue.] We will discuss gauge covariance and invariance in the context of linear, covariant 
gauges which is the class of gauges in which most studies are undertaken. 

The coupled fermion-gauge-boson DSE system is closed if the fermion-gauge-boson vertex is known. 
Further, it has become clear that, in Abelian theories, the fermion-gauge-boson vertex is crucial to 
ensuring the correct gauge parameter dependence of the objects one calculates in the DSE approach. 
This indicates the importance of this vertex in the DSE approach. This presents a problem, however, 
because, in Abelian gauge theories, the equation of which the vertex is the solution involves the kernel 
of the fermion-antifermion BSE which cannot be expressed in a closed form; i.e., the skeleton expansion 
of this kernel involves infinitely many terms. [The problem is even greater in non Abelian theories.] 



39 



As we have remarked, there have been no attempts to solve the DSE for the vertex. Instead, recognising 
the importance of the vertex, a number of authors have invested some effort into enumerating constraints 
that this vertex must satisfy and proposed Ansatze which satisfy them. The DSE for the fermion 
propagator and, in some cases, the coupled fermion-photon DSEs have then been solved using a given 
Ansatz. 

There are six general constraints one can place on the fermion-gauge-boson vertex (Burden and Roberts, 
1993): 

(a) It must satisfy the WTI; 

(b) It must be free of any kinematic singularities [i.e., expressing T^(p, q) as a function of p and q and 

a functional of the fermion propagator, S(p), then should have a unique limit as p 2 — > q 2 }; 

(c) It must reduce to the bare vertex in the free field limit in the manner prescribed by perturbation 

theory; 

(d) It must have the same transformation properties as the bare vertex, 7 M , under charge conjugation, 

C, and Lorentz transformations [such as P and T, for example]; 

(e) It should ensure local gauge covariance of the propagators and vertices; and 

(f) It should ensure multiplicative renormalisability of the DSE in which it appears. 

Criterion (b) follows from Ball and Chiu (1980) and criterion (c) is related to this since together they are 
necessary to ensure that the vertex Ansatz has the correct perturbative limit. The charge conjugation 
element of criterion (d) is essential since it constrains the properties of T^(p, q) under p <-> q. 

Although condition (a) is a consequence of gauge invariance, it is only a statement about the longi- 
tudinal part of the vertex, and says nothing about the transverse part. By itself it is insufficient to 
ensure condition (e) (Burden and Roberts, 1991). A well defined set of transformation laws, which 
describe the response of the propagators and vertex in quantum electrodynamics to an arbitrary gauge 
transformation are given in papers by Landau and Khalatnikov (1956) and Fradkin (1956) [LKF] [see 
Appendix [A|. These laws leave the DSEs and WTI form-invariant and one can, in principle, ensure 
condition (e) by choosing an Ansatz for T which is covariant under the action of the LKF transfor- 
mations. Unfortunately, however, the transformation rule for the vertex is quite complicated, making 
this procedure difficult to implement. Alternatively, since the LKF transformation rule for the fermion 
propagator and polarisation scalar are straightforward, one can use the LKF transformations to check 
a posteriori whether solutions for propagators obtained using a particular vertex Ansatz transform 
appropriately and hence whether the Ansatz used can possibly have been consistent with criterion (e). 

The fermion-gauge-boson vertex in an Abelian gauge theory can be written in the form (Ball and Chiu, 
1980): 

r M (p,g) = rJ c (p, g ) + x; f\p\q 2 ,p-q)T;( P ,q) , (3.80) 

8=1 

where T® c (p, q) is given in Eq. ( |3.64j ) and T l are the eight transverse tensors of Eq. (3.4) in Ball and Chiu 
(1980) of which those with i — 1, 2, 3 are symmetric under p <->• q and the remainder are antisymmetric. 
The requirement (d) above implies that all of the f l are symmetric except for f 6 which is antisymmetric. 
All vertex Ansatze which satisfy the WTI can be written in this form. 
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Criteria (a)-(d) and (f) have been used (Curtis and Pennington, 1990) in an analysis of QED4 to argue 
that a minimal admissible form of the vertex has /* = V i 7^ 6 and f 6 chosen so as to yield Eq. ( |3.65| ). 



Gauge Technique and DCSB. This is an ideal point to remark that there is an approach to solving 
the DSEs that differs substantially in method to the studies described above: the "gauge technique" 
(GT) (Salam, 1963). This approach, based in Minkowski space, assumes that the elements of the DSE 
approach [i.e., the propagators and vertices] have spectral representations in terms of which the DSEs 
are reformulated and then solved for directly. For example, in the context of studies of DCSB, it is 
assumed that a spectral density, p^, exists such that the fermion propagator can be written: 

S{p)= T du P -^± (3.81) 

J-00 f — UJ 

and that the fermion-gauge-boson vertex has a similar form. In fact, a commonly used vertex Ansatz 
in this approach is (Delbourgo and West, 1977) 

/oo I I 
dup^u)- 7m - , (3.82) 
-00 p — 00 q — to 

which has been used in an analysis of QED 3 (Waites and Delbourgo, 1992). For comparison it is useful 
to note that Eq. (|3.82|) corresponds, in Euclidean space, to the vertex 



r« T (p,g) = i(a At 5- 1 ( 9 )-5- 1 (p)a A4 ), (3.83) 

with a M = [7 • pj^ + 7^7 • q]/[p 2 — q 2 }, or alternatively: 

GT p 2 A(p) - q 2 A{q) , Ajp) - A(q) . B(p) - B(q) 
*r£ {p,q) = 5 5 ^7 M + o 2 — z 7'.P7 j u7-? + - 2 n — (7-P7m + 7m7^)- (3.84) 

The vertex, T^ T , is easily seen to satisfy criteria (a) to (d) and must therefore be of the form Eq. ( |3.80|) . 

Inserting Eqs. ( 3.81| ) and (3.82) into the DSE for the fermion propagator one obtains a linear equation 
for the spectral density. This is a feature of the gauge technique: it can reduce the fermion DSE to a 
linear equation. In fact, in practice, Ansatze for the propagators and vertices are almost always chosen 
so that the DSEs under consideration are linear integral equations for the spectral density. 



In a manner similar to that in which the Ball-Chiu Ansatz can be improved using criterion (f) it is 
also possible to improve Eq. ( |3.82| ) by adding a correction term which ensures that the fermion DSE 
is multiplicatively renormalisable in this approach (King, 1983) and ensures gauge covariance of the 
fermion propagator for small- and large-p 2 [but not for intermediate momenta]. 

A property of this approach which may prove undesirable, however, is that the equations for the 
Dirac-vector and -scalar parts of the spectral density decouple in the limit of zero Lagrangian- and 
renormalised-bare-mass. In this case one always has the chiral symmetry preserving solution 

S w (p) = -t7'J*rv(p) (3-85) 

and, in addition, it is also probable that the equation admits a dynamical chiral symmetry breaking 
solution which would have the form 

S NG (p) = - ll -pa^(p) + ^ s G (p). (3.86) 
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Here the notation indicates the Weyl (W) and Nambu-Goldstone (NG) modes respectively. [This was 
the case, for example, in the phenomenological QCD studies of Haeri and Haeri (1991).] We remark 
that in Eq. (|3.85| ) and Eq. (|3.86|) the vector part of the propagator is necessarily the same. 

The fermion DSE is the stationary point equation for the effective action discussed in in Corwnall, Jackiw 
and Tomboulis (1974) which, evaluated at this stationary point is (Stam, 1985; see also Sec. |7. 2f ) : 



V[S] 



d d p 



i 



trln[l-£(p)S(p)] + -tr[E(p)S(p)] 



(3.87) 



One might measure the relative stability of these extremals by evaluating the difference ^[S* 1 ^ ] — V[S W ]. 
For an Abelian gauge theory with Nj flavours of fermion one finds [for d = 3 or 4 since we use 4 
component spinors] that 



V[S NG ] - V[S 



Wi 



2N f 



d d p 

(2^jd 



111 



1 + 



1 <rj(p) 
P 2 Oy(p) 



> 0, 



(3.88) 



since it is reasonable to assume that as and oy are real for real Euclidean p 2 . [Since the equation for 
as is homogeneous, this difference can, in fact, be made arbitrarily large: erg — ► \crs-] Hence, based 
on this effective action [which is the same as the auxiliary field effective action at the stationary point] 
one finds that the gauge technique cannot support dynamical chiral symmetry breaking (Burden and 
Roberts, 1993). This result is true of any approach in which the equations for the Dirac- vector and 
-scalar pieces of the fermion propagator decouple. [We note that an extension of the CJT effective 
action to include a functional dependence on the fermion-gauge-boson vertex may allow for DCSB even 
when the equations decouple (Haeri, 1993).] 



Gauge Covariance in the Quenched Approximation. Criterion (c) has an important consequence. 
The unrenormalised fermion DSE in QED 3 and QED 4 can be written 



(ry • p + m)S(p) + e z 



d d q 



D^(p - qh»S(q)T v (q,p)S(p). 



(3.89) 



If one works in the quenched approximation, Eq. ( |3.10j ), then since 



/ dQ d - (d - 3) p ■ q + 2 = 

J (p- q) 1 \ (p- qy J 

it follows that, in Landau gauge, Eq. ( |3.89| ) admits the free propagator solution 

S(p) = — 
ff ■ p 



(3.90) 



(3.91) 



for any vertex that satisfies criterion (c). It follows directly from this that if a given vertex Ansatz is to 
satisfy the gauge covariance criterion, (e), then, for arbitrary £, the associated fermion DSE must have 
the LKF transform of the free field propagator as its solution [Eq. (|A.2Q 1. 

We can elucidate this further by considering the massless fermion DSE in configuration space: 

8 i (x-y) = 'y-df B S{x-y)+ 
e 2 J d d zd d x'd d y'- ffl (D^x -z) + d z d z A{x - zj) S(x - x')T v {z; x', y')S{y' - y) (3.92) 
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where we have explicitly divided the gauge-boson propagator into the sum of a transverse, gauge inde- 
pendent piece, Dj u , and longitudinal, gauge dependent piece, A. Making use of the WTI 

d,T,(z; x', y') = S~\z - y')5 d (x' - z) - 6 d (z - y')S-\x' - z) (3.93) 

and the identity J x , 7 • d x S(x, x')S' _1 (x / , z) = 7 • d x 5 d (x — z) , one obtains the massless DSE in the fol- 
lowing form: 

5 d {x-y) = ^-d x S{x-y) 



e 2 



d d z[y ■ d x A{x - z)]5 d (x - z) - [7 • d x A{x - y)}\ S(x - y) 



+ e 2 d d zd A x'd d y' 1 ^D T (x-z)S{x-x')T v {z ] x',y')S{y'-y). (3.94) 



/IV 

Now it is clear by inspection that if 



J d d zd d x'd d y' lii Dl v {x - z)S(x - x')r u (z; x', y')S(y' - y) = ; (3.95) 

then Eq. (|A.2p , with S(x;£ = 0) given in Eq. ( |A.1U| ), is a solution of the massless DSE; i.e., it is a 
solution if the last term on the right hand side of Eq. ( |3.89|) is identically zero in Landau gauge. [We 



remark that Eq. ( |3.90D is only relevant in quenched approximation: beyond quenched approximation 



the integrand is modified and the integral is nonzero so that, even given (c), the free propagator is not 
a solution and gauge covariance plus Eq. ( [3.94D doesn't require Eq. (|3.95|) .] 



As we have seen, most studies of the DSEs are undertaken in momentum space and it is a simple matter 
to transcribe Eqs. (|3~94l) and ( ggg) . We see that the solution of the DSE is LKF covariant if 



/ -^ D Uv-q)l,S{q)Y u {q lP ) = ^ (3.96) 



where Dj u (k) = (<Lj, — k^k v /k 2 ) /k 2 in the quenched theory, in which case the propagator satisfies: 



d d q ry • (p — q) 
(2vr) d (p-qY 



1 = * 7 " pS(p) + ie 2 J f2 ^ d ' 7 fp _ q) 4 l S <P) ~ S ^ ( 3 ' 97 ) 
in the covariant gauge fixing procedure. 



Equation (|3.96|) provides a much needed additional constraint upon the vertex function which, while 
not a full implementation of criterion (e), is nevertheless a restriction on the form of the transverse part 
of the vertex: 

(e') In the absence of dynamical chiral symmetry breaking; i.e., for as = 0, the vertex must be such 
that Eq. ( |3.96|) is satisfied, 



where D^ u (k) is the transverse part of the quenched photon propagator. This is a necessary condition 
which must be satisfied by any vertex Ansatz if it is to confer gauge covariance on the quenched QED3 
and QED 4 DSEs. 



We can apply this check in the case of the Curtis-Pennington vertex, Eq. ( |3.63| ), for which the fermion 
DSE in QED 3 is, in the absence of DCSB, 



-e 2 ^ yoo 1 (p 2 A{p)-q 2 A(q) p 2 A{p) + q 2 A{q) 
A[T>) ' l -\~^k q W)[ 2^ ln 
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p + q 



p-q 



(3.98) 



in which the right hand side is clearly zero in Landau gauge. Hence this vertex, or at least that part 
of it which contributes to the DSE, has the form necessary to ensure gauge covariance of the chirally 
symmetric fermion propagator. [We remark that here an infrared regulator, Air, is used to define d(k,p) 
in Eq. ( p.66| ) when B(k) = and the limit Air — > taken before evaluating the integral because it is 
finite.] 

It is possible to solve this equation analytically (Burden and Roberts, 1993). The solution, for £ > 0, is 

—— — = 1 — arctan ( — — ^ | , (3.99) 

A(p) 8vrp V e W 

as it should be since this corresponds to the LKF transform of the massless, free fermion propagator in 
QED 3 [see Appendix [Aj. 

The fermion DSE for QED 4 using the Curtis-Pennington vertex is given in Curtis and Pennington 
(1991) and can be written formally as 



2 2x9 2 , of 2 2^P A (P) 



o( P 2 - q y-, + % 2 " P 



p 2 ' q 2 A(q) 



(3.100) 



In Curtis and Pennington (1991) this equation was solved by introducing an upper bound on the q 2 
integral with the result A R (p 2 ;£) = A R (/x 2 ; £)(p 2 /^i 2 ) 7 , 7 = a £/(47r)/[l + a £/(87r)], which illustrates 
that Eq. (|3.63|) ensures multiplicative renormalisability of the chirally symmetric fermion DSE for all p 2 . 
However, this solution only agrees with the LKF transform of the free propagator at 0(£). This defect 
arises because a hard cutoff was used, which violates Poincare invariance and generates a spurious term 
in the equation. A detailed discussion of this can be found in Dong, et al. (1994), where the actual 
solution is shown to be A R (p 2 ;£) = A R (/i 2 ; £)(p 2 //i 2 ) 1 ' with v = a £/(47r), which agrees with the LKF 
transform at all £. 

It is interesting to note that the Ball-Chiu vertex alone, Eq. ( |3.64| ), does not satisfy criterion (e') (Burden 
and Roberts, 1993). This is also a defect of the minimal gauge-technique vertex, Eq. ( |3.83| ), however, 
the addition suggested by King (1983) goes some way to eliminating this problem. 



4 Four-dimensional Quantum Electrodynamics 

Many of the remarks made and studies described in Sec. |3| are directly relevant to QED4 . In the 
absence of a Lagrangian bare mass for the electron, the action for QED4 has the same chiral symmetry 
as massless QCD. This being the case, one of the main goals of DSE studies in QED4 over the last 
twenty years has been to answer the question of whether there is a phase in which chiral symmetry is 
dynamically broken and, if so, what is the order of the phase transition. [As we will see, this is related 
to the question of whether QED 4 is trivial or not.] Naturally, the tool used to address this question is 
the DSE for the renormalised electron self energy in QED4 which, in Euclidean space, is 

f d^q 

E(p) = 17 ■p(Z 1 - 1) - (m - Z 1 m ) + Z x e 2 J j^-^T fl {p,q)D fa/ {p - q-,^)S F (q)y v ; (4.1) 

a minor modification of Eq. ( |3.6|) with significant consequences. 

The first of these we have mentioned above. In Landau gauge and at 0(e 2 ) one has 

Z 1 =l. (4.2) 



44 



This result has been used in many studies to simplify this equation. Further, it has been used as an 
approximation in many model DSE studies in QCD. 



4.1 Dynamical Chiral Symmetry Breaking 

The second consequence is of particular use in the rainbow, quenched approximation to QED 4 . In this 
approximation, it follows from Eq. Q3.90| ) that the angular integral for the Dirac-vector component of 
the electron self-energy vanishes in Landau gauge so that, subject to Eq. (|4.2|) , Eq. (|4.1| ) admits the 
solution A(p 2 ) = 1 with B(p 2 ) obtained as the solution of the following nonlinear integral equation: 

, „ ^ , f A w d A q 1 B(q 2 ) . . 

where Auv is a cutoff which is necessary if m ^ 0. 
In addition, one may make use of the result: 

where / d£l± = ^2 Jo d/3 sin 2 (3 Jq d6 sin 9 J 2v dcf>, to reduce Eq. fl4.3|) to the following differential 
equation: 

, B » W+ 2 B ' W + |-^M_^0, (4.5) 
where x = p 2 , with the boundary conditions 

lim ( — \x 2 B(x)] l =0 and lim ( — \xB(x)]\ = m . (4.6) 
x->o\dx L v n ) x^Al v \dx L y n ) K ' 

The problem of studying DCSB in the rainbow, quenched approximation to QED4 is thus reduced to a 
study of this differential equation. 



This is an opportune point to remark that Eq. ( |4.4|) has been used as the basis for an extensively used 
approximation in phenomenological DSE studies of QCD and QED 4 beyond quenched approximation. 
Many authors have assumed that it is reasonable to write [i.e., the angle approximation] 

dil A — — w 9{p -q) — — + 9{q -p) — — , (4.7) 

(p — q) A p z q z 

with the function a(p 2 ) assumed to be monotonically decreasing on [0, 00). This approximation has 
been critically analysed (Roberts and McKellar, 1990) and found to be reasonable only if a(0) is not 
large; i.e., of 0(1) or less. 



Miransky Scaling. Equation ( |4.5|) has been studied by Fukuda and Kugo (1976) and the results have 
a direct analogue in those of Sec. [37|. If one sets m = then Eq. (|4.5|) is obviously scale invariant. 
In this case the scale is set by choosing -B(O) = 1 and the solution then has the following power series 
expansion around x = 0: 

B(x)=l-—x-— --8 )x 2 + ... . (4.8) 



871 QAlT V 7T 
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The asymptotic behaviour of the solution for large-x is easily found to be 



B(x) ~ 



71 

Q£ 3 



COS 



-y3a/7r — 1 In x 



71 

«>3 



(4.9) 



If one demands that the ultraviolet boundary condition of Eq. Q4.6D be satisfied at some finite Auv then 
it is clear that a nonzero solution for B(p 2 ) is admitted only if 

a>a c = \ (« 1.047) (4.10) 

since the form of the solution valid for a < a c cannot satisfy this condition. In this case one has [see 



Sec. ET5 
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B(0) = A uv exp I - + 5 , (4.11) 

yja/a c -l J 

where 5 is some phase, and the conclusion is that chiral symmetry may be dynamically broken in 
rainbow-quenched QED 4 if the coupling exceeds a c . A study of the stability of the nonzero solution for 
B, the solution which corresponds to a Nambu-Goldstone mode realisation of chiral symmetry, relative 
to the solution B = 0, which corresponds to a realisation of chiral symmetry in the Wigner-Weyl mode, 
using the auxiliary field effective action (Roberts and Cahill, 1986) shows that the Nambu-Goldstone 
mode is indeed dynamically favoured for a > a c . 

As the critical coupling is approached, -B(O) and (ipip) are interchangeable as order parameters of DCSB. 
Equation ( f4.11| ) shows then that, in rainbow-quenched approximation, QED 4 has an infinite order phase 
transition. In models such as this one Eq. ( [4.1 1|) is often referred to as "Miransky Scaling" (Miransky, 
1985a). 

It has been pointed out by Roberts and Cahill (1986), however, that in the limit m = the integral 
equation admits a solution even when the cutoff is removed; i.e., when Auv — °°- in this case the 
ultraviolet boundary condition can be satisfied for any value of a and hence rainbow-quenched QED 4 
admits DCSB for all values of the coupling. [See also, Atkinson and Johnson (1987).] This problem 
has also been studied using the vertex Ansatz of Eq. ( |3.63| ) (Curtis and Pennington, 1992); i.e., beyond 
rainbow approximation. In this case it is convenient to solve the integral equations directly because 
of the complicated structure of the vertex but the same result is obtained; i.e., a dynamical mass is 
generated for all values of the coupling. 

In this light the existence of a critical coupling in massless QED 4 might be seen as the result of the 
imposition of a cutoff which can be viewed as an auxiliary condition, not intrinsic to the model, but 
relevant to the consideration of whether the model is a reasonable representation of QED 4 . For example, 
it might be argued that the appearance of the Landau ghost in perturbation theory suggests that the 
quenched approximation cannot possibly be valid at all momentum scales and hence that introducing 
a cutoff is a simple manner in which to model the effect of vacuum polarisation contributions. 

However, dynamical chiral symmetry breaking in theories with ultraviolet divergences is a subtle prob- 
lem (see, for example, Fomin et ai, 1983, pp. 9-10 and pp. 64-65; Miransky, 1985a and 1985b). In 
a recent reanalysis of this problem (Curtis and Pennington, 1993) it is argued that dynamical chiral 
symmetry breaking in quenched QED 4 cannot be studied simply by setting m = in Eq. ([4. 1| ) and 
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allowing the momentum integration to range over [0, oo) but rather that renormalisation plays an im- 
portant role. The vertex Ansatz of Eq. ( |3.63| ), which ensures multiplicative renormalisability of the 



fermion DSE in QED 4 , is therefore well suited to this study. These authors suggest that in order to 
ensure conservation of the axial vector current in the renormalised theory, which can be interpreted as 
a definition of masslessness in the theory, one must take the cutoff to infinity and the bare mass to zero 
in such a way that their product approaches zero. This ensures that, since the bare composite operator 
(ijjjsip) a uv depends on the cutoff, 

8,Jl = 2m (^) Auv = 2(Z 2 - 1 m ) $7bV0r Au ^°° , (4.12) 

where Z2 is the mass renormalisation constant, and can be implemented through a definition of the 
renormalised theory with = and with a cutoff Auv < 00. In this way a Nambu-Goldstone realisa- 
tion of chiral symmetry; i.e., DCSB, was found to be possible only for 

a > a c « 0.93 (4.13) 

and Miransky scaling was observed at the critical coupling. This result was only weakly dependent on 
the gauge parameter, a result that one can attribute to the efficacy of the vertex Ansatz [see Sec. |3[ 



Mean Field Scaling - Vacuum Polarisation Ansatz. A question that arises is whether the existence 
and order of the phase transition in QED 4 is sensitive to the nature of the approximations. The studies 
described above suggest that going beyond rainbow approximation does not qualitatively affect the 
results but what is the effect of the quenched approximation? Given the results discussed in Sec. ET5 



it is clear that going beyond quenched approximation can introduce another dimensionless parameter 
into the theory, the number of electrically active fermions, with respect to which there may exist critical 
behaviour. A number of authors have addressed these questions. 

The first step is to choose a form for the vacuum polarisation. Motivated by the large-fc 2 form of the 
vacuum polarisation in perturbation theory at one loop order, the usual choice is 



uj{k 2 ) = 1 + II(£; 2 ) 



, N f a 1 
1 H J — In 

3vr 



A 2 
k 2 + fi 2 



(4.14) 



where Nj is the number of active fermions and /j, 2 is an infrared regulator. In some studies /z 2 is set to 
(Kondo, 1990; Gusynin, 1990; Kondo, 1992; Kondo and Nakatani, 1992b; Maris, 1993, Sec. 4.3) while 
in others it is related to B(0), e.g., /z 2 = I?(0) 2 (Oliensis and Johnson, 1990), in order to incorporate 
some information from the fermion DSE in the expression for the vacuum polarisation. In any event it 
has little effect on the results. [Equation ( 4.14 ) can be obtained from Eq. ( 3.13 ) with d = 3 — * 4 using 
bare fermion propagators and a bare fermion-photon vertex.] Such a form has been used in almost all 
studies to date in concert with the approximation of Eq. ( |4.7[ ) in which case one obtains A(p 2 ) = 1 and 
the following differential equation for B(x = p 2 ): 



xB"(x)+B'(x) 



1 + 



to + to 
to 



1 + 



u' uj [J + 



UJ + UJ 1 [uj + Uj'\ 



n 2 



3a uj + uJ B(x) 
47r uj 2 x + B(x) 







(4.15) 



where u' and uj" represent derivatives of uj with respect to the variable t — In x, with boundary condi- 
tions 



x 2 B'(x) 



2 = 



and \(uj' + uj) B(x) + uj x B\x)} 



x=A 2 







(4.16) 
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This differential equation can be solved numerically (Oliensis and Johnson, 1990) and one finds that, 
with Nf = 1, it admits a nonzero solution for 

a > a c » 2.00 (4.17) 

and that the phase transition is of finite order. Further analysis (Kondo, 1990; Gusynin, 1990; Kondo, 
1992; Kondo and Nakatani, 1992b; Maris, 1993, Sec. 4.3) shows that in this case, in contrast to 
Eq. (TO , 

B(0) ~ A uv (a-a c )* ; (4.18) 

i.e., one has a second order phase transition with the classical mean field value for the critical exponent. 
Both (ipip) and B(0) behave in the same way near the critical coupling (Kondo and Nakatani, 1992b). 

It is obvious (Kondo and Nakatani, 1992b; Maris, 1993, Sec. 4.3) that if Eq. (|4.7| ) is not used in order 
to obtain an approximate differential equation but instead Eq. ( |4.14|) is used directly in the integral 
equation form of the fermion DSE in QED4 , then one encounters a problem with the Landau ghost: 
the kernel of the integral equation, which involves 

1 + U(k 2 +p 2 -2kp cos9) (4.19) 

where the argument of IT can attain a maximum value of 4 A^y, has a nonintegrable divergence unless 

3tt 
2~ln~2 

which restricts the number of flavours that one may consider. A first estimate based on Eq. ( |4.17| ) is 
that Eq. Q4.14 ) will not be a viable Ansatz for Nf > 3. 



a Nj < 771 - TT w 6.80 (4.20) 



Numerical studies of the coupled integral equations for A(p 2 ) and B(p 2 ) using Eq. ( |4.14 ) in rainbow 



approximation (Kondo and Nakatani, 1992b; Maris, 1993, Sec. 4.3) find nontrivial solutions only for 
Nf = 1,2 with critical couplings (Maris, 1993, Sec. 4.3) 

a c (N f = 1) « 2.07 and a c (N f = 2) « 2.82 . (4.21) 

Kondo and Nakatani (1992b) used an Ansatz for the vacuum polarisation that differs slightly from 
Eq. fl4.14j ), one which reflects the possible regularisation ambiguities, and in this case the value of the 
critical coupling depends on the choice of the additional parameters. Indeed, it is possible to choose 
the additional parameters and the number of active fermions, Nf ! , such that one does not have DCSB 
irrespective of the strength of the coupling. Avoiding such a complication, the numerical solution of the 
coupled integral equations leads to the same mean field scaling behaviour as the approximate differential 
equation. 



Mean Field Scaling - Solving for the Vacuum Polarisation. Given these last remarks, an obvious im- 
provement would be to solve the coupled fermion and photon DSEs in Landau gauge using rainbow 
approximation: 

A{p 2 ) = 1 (4.22) 



+ a 



dk 2 k 2 



A{k 2 ) 



Air k 2 A(k 2 ) 2 + B(k 2 ) 2 Jo it 



T d6 9 k 
2 sin 2 9 - 



3 (k — p) 2 cos 9 — 2 kp sin 2 9 



B(p 2 ) 
U(p 2 ) 



m +a 



--a N f 



dk 2 k 2 



B(k 2 



* d9 



4tt k 2 A(k 2 ) 2 + B(k 2 ) 2 Jo it 

2 .2 — in A I U2 



2 sin z 9 



p (k - p) A [1 + n((ife - p) 2 )] 
3 



(k-p) 2 [l + U((k-p) 2 )] 



A uv dk' 2 k 



Att 



* d9 , n A(k 2 + )A(k 2 _) (M8 cos 2 9-2]- 
, it \k\A(k\) 2 + B{kl) 2 } \klA{kl) 2 + B(k 2 _ 



(4.23) 



(4.24) 
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where ao is the bare coupling, mo is the bare electron mass and k 2 ± = k 2 + \p 2 ± kp cos 6. Such a study 
has been undertaken by Rakow (1991). 



The main body of the study concentrates on Nf = 1 and it is assumed that Z\ = 1, Eq. (|4.2|) , while 

Z 2 = -t^t and Z 3 = \— , (4.25) 

A(0) i + n(o) ' v ; 

which follows from = Z{Sf , D^ v = Z 3 D^ U and Eq. ( |2.34| ) in Sec. pTT] , with renormalisation at p 2 = 
rather than on the mass shell. Clearly, deviations from A(p 2 ) = 1 measure the violation of the Ward 
Identity in rainbow approximation, Eq. ( |2.32| ). The results are then analysed in terms of a renormalised 
mass and coupling: 

mR = 1(5) = 22 jB(0) and " R = A(o) 2 [i + n("" = 22 23 a ° ' (4 ' 26) 



respectively, with ckr measuring the effect of charge screening due to fermion loops. 
It is found that with m = and for Nf = 1 there is a critical coupling 

a c sa 2.25 (4.27) 

such that 

(Wauv = A w ("o - occ) 1 * ; (4.28) 

i.e., that one has mean field scaling behaviour in the model, but that the critical exponent approaches 
the mean field value of | very slowly. [See Creswick et al. (1992) for a pedagogical discussion of phase 
transitions and renormalisation group methods in physics.] It has been shown that 

(fij>) ~A UV ™ R (4.29) 

and hence, identifying a as a(Auv), Eq- ( fl-28[ ) leads to the following /5-function: 

/3(a) = A UV t^— In a(A uv ) = -2 [a(A uv ) - «J (4-30) 

CAuv 

which suggests that unquenched, rainbow approximation QED4 has an ultraviolet stable fixed point at 
«o = a c ; a result supported by lattice simulations (Kogut et al, 1988a, 1988b). 

In the limit ao — - ► a c from above it is found that — > logarithmically; i.e., recall that vanishes 

aS Q!o ~~ ^ a c; 

11 . . 

"roc-- cx — (4.31) 

In m.R In^ 

where ^ = Auv/ m R is the correlation length in QED4 . The continuum limit of the theory is defined 
as the limit £ — > 00. It follows from Eq. fl4.31|) that the continuum limit of unquenched, rainbow 
approximation QED 4 is trivial in the sense that it does not involve photons interacting with charged 
fermions, since it is not possible to obtain an infinite correlation length at nonzero coupling. 

It is clear from these observations that solving for the vacuum polarisation yields the same results as 
the Ansatz of Eq. ( [4.14|) with only minor quantitative differences. The problem of the "Landau-ghost" 
is not removed by a simultaneous solution of the fermion and photon DSEs in rainbow approximation. 
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Interacting Continuum Limit? The existence of a phase in which chiral symmetry is dynamically 
broken suggests that at strong coupling [a > a c ] in QED 4 there are Goldstone modes; i.e., massless, 
pseudoscalar e + -e~ bound states. In this connection, one notes that in addition to studying the DSEs 
for the fermion and photon propagators in rainbow approximation, Rakow (1991) solved the following 
approximate form of the integral equations for the pseudoscalar and scalar electron-positron scattering 
amplitudes: 

S 5 (pV) = a,/ o -2s.n ^ 2[1 + n(t2)] (4.32) 







4tt k 2 A(k 2 ) 2 + B(k 2 ) 2 Jo ti g 2 [l + II(g 2 )]' 



Sl(pV , )=ao /;^ 2sin ^_^_ (4.33) 



,^dk 2 k 2 on2 2 . k 2 A{k 2 ) 2 - B{k 2 ) 2 r d6 n . - 
+«o / —— S,(k 2 , r 2 ) - 2 A )> ■ I - 2 sm 2 (9 



o 



4tt ±v ' 7 [k 2 A(k 2 ) 2 + B(k 2 ) 2 ] 2 Jo ti g 2 [l + II(g 2 )] 



(4.34) 



where t 2 = p 2 + r 2 — 2 p r cos 9 and q 2 = k 2 + p 2 — 2kp cos 0. These equations have been averaged 
over p ■ r and describe [ip{p)', ip{— p)] - * [VK 7 *)) V'l - r )] scattering. 

In studying the solution of these equations it was observed that, at fixed values of «r, these amplitudes 
did not scale; i.e., there was a dependence on p/m^. Indeed, the pseudoscalar amplitude manifested a 
contribution that diverged in the limit mo — > obviously as a result of the formation of the Goldstone 
mode. This lack of scaling at constant «r suggests that another operator is relevant in the continuum 
limit. 

The fact that another operator is relevant in the continuum limit can also be inferred from Eq. (|4.9|) . 
In the absence of interactions the mass operator, ipip, diverges as A^y: 

dsstv[S(p)] A "oT°° A^y . (4.35) 

In the interacting case, using Eq. fl4.9|) , it diverges as A]^/^~ 3a ^^ from which one can infer that its 
interaction-induced anomalous dimension is 



3a , 



7T 

in quenched approximation and hence that the net dimension of the operator is 



In quenched, rainbow approximation [also called the planar approximation] the dimension of [ipijj] 2 is 
just 2 d^ and one is therefore lead to infer that, in the continuum limit [obtained, in this approximation, 
as a — > |] drp^ 2 ~~ * 4 and hence that it becomes a relevant operator. This inference is supported by an 
analysis of quenched QED 4 which goes beyond the rainbow approximation, using the dressed electron- 
photon vertex of Eq. ( |3.63| ) (Atkinson et al, 1993). 

Rakow (1991) studied QED4 with the additional chir ally-invariant, four-fermion interaction term 



,4 ^0 
a x — 
2 



ip(x)^(x)ip(x)ip(x) — "ip(x)j 5 ip ( x )' l P{ x )l5' l P ( x ) (4.38 
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which is just the interaction of the Nambu-Jona-Lasinio model [NJLM] (Nambu and Jona-Lasinio, 
1961). With this term, Eqs. Q4.23| ) and ( 4.32|) receive the additive corrections: 



Go f^v dk 2 k 2 B(k 2 ) dk 2 k 2 S 5 (k 2 ,r 2 ) 

and G + / — — tJijtv?. I U7T^3. > ( 4 - 39 ) 



n Jo An k 2 A(k 2 ) 2 + B(k 2 ) 2 u Jo An k 2 A(k 2 ) 2 + B(k 2 ) 2 

respectively, and the other equations remain unchanged. [This model was first considered by Leung 
et al. (1986).] With this additional interaction, Eq. ( |4.38| ), all of the functions scale with p/mR, the 
fermion and photon propagators and the pseudoscalar and scalar scattering amplitudes, showing that 
all of the interactions that are relevant in the continuum limit are now included. 

With «o = one simply has the NJLM in which chiral symmetry is dynamically broken for GoA^y > 
An 2 . In general the model has a two-dimensional phase diagram in the (GoA^y, ao)-plane. From Fig. 14 
of Rakow (1991) one infers that the model has DCSB on a set which, in the first quadrant, can be 
represented approximately as: 

Si = |(G A^ v ,a ) : GoA^ v > 0, a > 0, + ^ > 1 j . (4.40) 

This is only a part of the picture, however, and for a given value of «o £ [0, oo) there is always a value 
of G = Gq G (— oo, An 2 ] such that the model exhibits DCSB for G > Gq\ as a is increased Gq moves 
towards — oo. The domain in the (GoA^y, ao)-plane on which chiral symmetry is dynamically broken 
is similar to that found in quenched, rainbow approximation (Leung et al, 1986). 

One may now return to the question of triviality. 

One notes that, for a given value of «r, the correlation length increases as Go decreases and takes its 
maximum value at Go = — oo but that at any point not on the phase boundary this maximum value 
is finite. On the phase boundary the renormalised coupling is always zero and hence, even with the 
four-fermion interaction of Eq. ( |4.38| ), the theory doesn't have a continuum limit of interacting fermions. 

There is another possibility, however. There has been speculation for some time (Miransky, 1985a; Kogut 
et al., 1988a, 1988b) that QED 4 has an interesting continuum limit of interacting, strongly-bound e + -e~ 
states. One can convince oneself of this possibility simply by looking at the effective action for QED 4 
as obtained in Roberts and Cahill (1986), for example. Expanding this action, following the procedure 
of Roberts et al. (1988), for example, one finds interactions between the e + -e~ bound states whose 
strength is characterised by 

const. < oo , 
9y = j (4.41) 

where 



2tt 2 Jo (s + B(s) 2 ) 2 

It is clear then that as one takes the cutoff to infinity, in order to recover the continuum limit, the 
theory will be one of interacting bosons if / < oo and free bosons if it is not. 



It is quite clear from Eq. (4~S) that in quenched approximation / is indeed finite. In going beyond 
quenched approximation, however, one must use numerical solutions, or some approximation to them, 
in which case one finds that / = oo (Gusynin, 1990; Kondo, 1990; Maris 1993); a result which suggests 
that the continuum limit is one of noninteracting bosons. 
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4.2 Analytic Structure of the Electron Propagator 

The differential equation, Eq. Q4.5| ), has been studied by Atkinson and Blatt (1979). This study revealed 
a "pathology" of the ladder approximation, used in all of the studies discussed in the previous section 
[except Kondo (1992)], that has recently attracted a good deal of attention. Solving Eq. ( [4.5| ) one finds 
that instead of an electron propagator with a real branch cut on the timelike p 2 axis, as expected on 
physical grounds, the electron propagator in rainbow-quenched QED 4 has a pair of complex conjugate 
branch points in the complex p 2 plane whose position depends on a. The recent studies of rainbow 
approximation model DSEs for QCD (Stainsby and Cahill, 1990; Maris and Holties, 1992; Stainsby 
and Cahill, 1992; Stainsby, 1993; Maris, 1993) suggest that this pathology is an artifact of the rainbow 
approximation; a conclusion supported by the extended analysis of rainbow approximation QED4 by 
Maris (1993, Sec. 4) whose analysis also suggests that this "pathology" survives the inclusion of the 
one-loop vacuum polarisation, itself calculated in rainbow approximation. 



It may be that the resolution of this problem lies simply in using a dressed fermion-gauge boson vertex; 



the possibilities and constraints on which we discussed in Sec. |3.7| . However, there have been no studies 
to date that test this hypothesis directly in QED 4 . It should be noted, however, that, in a QCD-based 
model, Burden et al. (1992b) demonstrated that the structure of the vertex does significantly affect the 
singularity structure of the fermion propagator. 



In Sec. |3j we briefly discussed the gauge technique. This has also been applied to QED 4 by Delbourgo 
and West (1977) who advocate an iterative solution of the DSEs for the spectral densities. 



At lowest order, as it is defined in this approach, the only equation of interest is 



Z' 1 = (p l -m G )S {a \p)-ie 2 



(2k)' 



S {0) (p) If) (p, p-q) S<® (p - q) D""<°) (q) lu 



where is given by the expression in Eq. (ET5ID , by Eq. ( ggg) and 



D^\q) 



with £ the gauge parameter. 



qi 1 gz 



(4.43) 



(4.44) 



Making use of these definitions Eq. ( |4.43|) becomes 



Z" 1 



00 <L;^M 

-00 "p — OJ 



-mo + E(p,w)) 



with S(p, u) obtained in the lowest order of perturbation theory: 



e 2 (p 2 



Im E(p, a;) 
Making use of the definitions 

Z^ 1 = I dm p^{uj 



16 7rp 3 



£(p 2 + l j 2 )-(Z + 3)pu; 



and 



m 



duo to p^(uo) 



then Eq. (|4.45|) yields the following result in Landau gauge: 
S ( ~°\p) = - 1 



m 



m 



2\ 2r l 



r(i-77) 2 r(i + 2£) 



F 1 - 77, 1 - 77; 1; 



V 

m 2 



+ 



rj 



m 



F 1-77,2-7752; 



(4.45) 

(4.46) 
(4.47) 

(4.48) 

p 2 * 
m 2 
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where fi 2 is an infrared cutoff and r] = — 3a/(47r). The next iteration in the procedure uses this result 
coupled with Eq. ( |3.82| ) and Eq. (|4.44|) in the DSEs for the photon propagator and electron-photon 
vertex. 

The interesting feature of Eq. ( |4.48|) is that the fermion propagator has a branch cut at p 2 = m 2 ; i.e., 
it has the singularity structure expected on physical grounds. This is a feature that it has in common 
with the propagator obtained with the solution of a linearised form of Eq. ( ^4.5| ) (Maris, 1993, Sec. 4.2): 

x (x + m 2 ) B"(x) + 2 (x + m 2 ) B'(x) + B(x) = , (4.49) 

47T 

with m = B(0), which is a hypergeometric differential equation whose solution is 

„. s „ | 1 1 3a 1 1 3q x \ . M 

^)=^F^- + -^--,---^--;2;--^ . (4.50) 

This is, perhaps, not surprising since the operational procedure in the gauge-technique is to ensure that 
one has linear equations in the spectral density. The study by Maris (1993, Sec. 4.2.3) emphasises that 
the singularity structure of the solution of the nonlinear differential equation is quite different from that 
of the solution of the linearised equation: it has complex conjugate branch points. 

The application of the gauge technique to QED 4 has not been pursued further. It is not known 
whether the procedure, outlined briefly here, converges nor whether the desirable singularity struc- 
ture of Eq. ( |4.48| ) survives further iteration. Given the results of Stainsby and Cahill (1990), Maris and 
Holties (1992), Stainsby and Cahill (1992), Stainsby (1993, Sec. 3) and Maris (1993, Sec. 4) it would 
be interesting to answer these questions. In addition to the limitations mentioned in Sec. |3.7| , however, 
one must add that the procedure advocated is quite difficult to implement. 

At present, in the study of Abelian gauge theories the ingredients that are necessary and/or sufficient 
to ensure physical singularity structure in the propagators remain unknown. 



5 Gauge Boson Sector of QCD 



We have briefly discussed the DSEs for QCD in Sec. |^2| with the equation for the vacuum polarisation 
and gluon propagator represented diagrammatically in Fig. 2.5. In this section we will review the at- 
tempts that have been made to solve these equations, the approximations used and the results obtained. 
[See also, Hadicke (1991).] This is important because in the application of DSEs to hadronic physics 
a physically reasonable form for the quark propagator is needed and, in practice, this is obtained by 
solving the quark DSE with Ansatze for the gluon propagator and quark-gluon vertex, with the former 
often being motivated by studies of the DSE for the gluon propagator. 

Before proceeding, however, this is an opportune point to make a number of observations. We use our 



Euclidean space conventions throughout this section, see Sec. 2.3 



In an arbitrary covariant gauge, specified by £, the gluon propagator can be written as 



^(, HV _<^}_L_ + ^ (5 ,; 
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where, as usual, IT(g 2 ) is the vacuum polarisation. It will be observed that only the transverse piece of 
the gluon propagator is modified by interactions. This is a result of BRS invariance which yields, in 
the covariant gauge fixing scheme, the STI (Pascual and Tarrach, 1984, pp. 42-45) 

q,D, u (q)=^ . (5.2) 

This is an important constraint. 



The least ambiguity and/or uncertainty exists in the form of the gluon propagator at large spacelike-g 2 
since this is the domain on which perturbation theory and the renormalisation group can be easily 
applied. Indeed, the evaluation of the asymptotic form of IT(g 2 ) is a textbook exercise (for example, 
Pascual and Tarrach, 1984, pp. 70-76). Determining the form of the propagator at smaller spacelike-g 2 
is, however, a difficult problem. 

In this connection, one can employ the renormalisation group in QCD to obtain an approximate relation 
between the vacuum polarisation and the running coupling constant in QCD. In Landau gauge, the 
renormalisation group equation for 

^^TTW (5 - 3) 

is, with Za = (Zf h )~ 2 (Z^ h ) 2 where Zf h and Zf are renormalisation constants for the ghost-gluon 
vertex and ghost wave function, respectively, 



A R (g 2 ;a(/x)) = (5.4) 



where (3(a) is the QCD /3-function: 



£(«(/*)) = ^ ^_J1 + QUo^/vr) 2 ) . (5.5) 



7T 6 

If one assumes that 

Zf h = Zf , (5.6) 

which is sometimes called the "Abelian approximation", then the last term on the left-hand-side of 
Eq. ( |5.4p vanishes and A R satisfies the same renormalisation group equation as «r. Since at large 
spacelike momenta one has A(g 2 — > oo) = a then the equivalence of the renormalisation group equations 
leads to 

A( ' 2) ^(xt) ; (57) 

i.e., to an identification of Ar with the running coupling constant in QCD and to 

/2 n , <\ _ J r \ 47ra(g 2 ) 

[g D, v {q)) K « | Kv - -Ar j ~^T- + ^ ■ ( 5 - 8 ) 

[This argument is a minor modification of that presented in Bar-Gadda (1980); see also, Itzykson and 
Zuber (1980, Chap. 13).] Equation ( |5.8| ) provides a prescription for an extrapolation of the known form 
of the quark-quark scattering kernel at large spacelike-g 2 (Marciano and Pagels, 1978) to the entire 
g 2 -domain on which the running coupling constant is known. 
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At two-loop order 



a 2 {q 2 ) = «i(g 2 ) 



, 1 153- 19iV f aAq 2 ) , 
1 — m 



where 



2 33-2JV/ 



7T 



- In —7T 

2 A 2 



QCD , 



ai(g ) 



(5.9) 



(5.10) 



ln(g 2 /A^ CD ) 

with [cZjvf = 12/(33 — 2Nf)] [see, for example, Pascual and Tarrach, 1984, pp. 128-131]. The renormal- 
isation group invariant mass scale in QCD, Aqcd ; is determined by fitting data in a number of high 
energy experiments and, in the MS scheme and with four quark flavours, 



A QCD = 0.20 ± ° Q f 8Q GeV 
(Particle Data Group, 1990). It is worth noting that at <7 2 /Aq CD = 100; i.e., y/q 2 ~ 2 GeV, 

0(2 — Oil 



(5.11) 



-0.3 : 



tt'2 



(5.12) 



ct 2 ~ 0.25. This suggests that the three- loop contribution may, already at this q 2 , not be negligible and 
that the strength of the interaction extracted from this expression for the running coupling constant is 
unreliable for q 2 less than this. 



5.1 Infrared Behaviour of the Gluon Propagator 

The discussion above illustrates that the tools of perturbation theory cannot provide information about 
the structure of the gluon propagator at small spacelike-g 2 . This region is very important as there is an 
expectation that the structure of the gluon propagator at q 2 ~ has important implications for quark 
confinement: in an imprecise way one might say that the behaviour of the q — q interaction in this 
region determines the long range properties of the q — q potential and hence incorporates the physics 
of confinement. [The discussion of confinement in QED 3 , Sec. ^, illustrates this.] On this domain 
nonperturbative techniques are necessary and studying the DSE for the gluon propagator is one such 
approach. 



Singular in the Infrared? Axial gauge studies. Some often quoted studies of the DSE for the gluon 
propagator (Baker et al, 1980a, 1980b, 1983) adopt the axial-gauge formalism which implements the 
constraint n ■ A a = 0, with a the gluon-field's colour index and n M an arbitrary four-vector: n ■ n > 0, 
and in which ghost fields are absent, Sec. |2~2| . [With our Euclidean conventions, temporal gauge is 
specified by n ■ n < 0.] In this approach the gluon propagator has the form 

D»u(q,l) = Fi(q, 7 )M flu (q,n) + F 2 (q, 1 )N^(q,n) (5.13) 

with 

M, u (q, n) = 8, v - + q » n » + n 2 -Mi_ and N, v {q,n) = 8„ v (5.14) 

q ■ n [q • n) z n z 

and where 7 = [q ■ n] 2 /[q 2 n 2 } is the "gauge parameter". The free propagator in axial gauge has Fi = 
—q~ 2 and F 2 = 0. As remarked in Sec. |2~2|, the benefit of axial gauge is the absence of ghosts; the 
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drawbacks are the loss of Lorentz invariance and the gauge-dependent singularities in the propagator 
at q ■ n — 0. 

In axial gauges the DSE for the gluon vacuum polarisation, given diagrammatically in Fig. 2.5, does not 
have a ghost-loop term. Further, because of the structure of the bare 4-gluon vertex, the term involving 
the dressed 4-gluon vertex vanishes in the contraction n M IL^. Hence one can obtain a simplified integral 
equation for n^U^n^ which only receives contributions from the first two diagrams and the last diagram 
on the right-hand-side of Fig. 2.5. This is the approach commonly taken to simplify the analysis of the 
gluon DSE. Further, fermion loop contributions are often neglected leaving only the second and last 
diagrams and thus a closed equation for the gluon vacuum polarisation once the dressed 3-gluon vertex 
is known. 



The STI for the 3-gluon vertex, which takes the following simple form in axial gauges: 

Pa 9, r) D w {q) D va {r) = D pa (r) - D pa (q), 



(5.15) 



fixes its longitudinal part, however, the transverse part remains unknown. [The colour structure of the 
vertex is, of course, f abc .] One way to proceed is to neglect the transverse piece entirely and simply 
solve the STI which leads to the Ansatz (Baker et at, 1980a) 



^X/11/ 



(p,q,r) + E^ vX (q,r,p) + 



r,p,q) 



(5.16) 



where 



^X/11/ 



(p, q, r) 



5 



Xfj, 



Pv 



p 2 F 1 (p 2 ,>y p ) q 2 F 1 (q 2 ,'y q ) 



+ 



p 2 F 1 (p 2 ,-y p ) q 2 F 1 (q 2 ,>y q ) 



(5.17) 
q v -Pv 

p2 _ g2 



which is free of kinematic singularities. At this point the further assumption F 2 = has been made; 
i.e., that the dressed gluon propagator has the same tensor structure as the free one. Substituting this 
into the DSE and regularising to remove ultraviolet divergences one obtains an equation of the form 



1 



P 2 + J(&?K(p,q,n)F 1 (q, lq ) 



F x {p 2 , lp ) L(p, q, n) F{q, lq ) F{r, 7p ) 



(5.1* 



where K and L are complicated functions of their arguments. Baker et al. (1981b) constructed a 7 
independent, few-parameter form of F\ which satisfied this integral equation to an accuracy of a few 
percent. This form has the asymptotic behaviour: 



Fi{p\l P 



2 \ P 2 ^0 1 j n f __2 _ \ P 2 ~oo 



P S 



and Ft(p , 7 P ) 



p2 



21-H/16 



(5.19) 



The possibility that the gluon propagator has a double-pole at the origin is appealing to many because 
it implies area law behaviour of the Wilson loop (West, 1982) which is often regarded as a signal of 
confinement. In this connection, the "potential" obtained from a propagator with this form in the 
infrared, using the four- dimensional analogue of Eq. ( |3.4| ), has a linearly rising piece at large distances. 
[This Fourier transform can be defined as the solution of the differential equation obtained by operating 
on it with y 2 .] 
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The gauge technique can also be employed to solve the gluon DSE (Delbourgo, 1979a, 1979b, 1981; 

A/if 



Atkinson et al, 1983; Delbourgo and Zhang, 1984). Using an Ansatz for r| that is consistent with 



the STI, Atkinson et al. (1983) obtained the following equation; 

n 2 (l- 7 ) Fi(p 2 ) -1 -p 2 =/ dsp{s)R{p,s,n) (5.20) 
L J Jo 

where R(p, s,n) is a complicated function of its arguments and 

F ^ = r«°jMr- (5 ' 21) 



Equation ( |5.20|) is a nonlinear integral equation for the spectral density. The nonlinearity arises because 



of the demand that the Ansatz for the 3-gluon vertex satisfy the STI and distinguishes this analysis from 
the other gauge technique analyses referred to above. Analysing this equation, Atkinson et al. (1983) 
concluded that the axial-gauge gluon propagator is likely to have a q~ A singularity in the infrared but 
that F 1 is unlikely to be independent of 7, as claimed by Baker et al. (1981b). 

There is a caveat to be borne in mind when considering these axial-gauge studies. It has been argued 
(West, 1983), from a general study of the properties of the spectral representation of the gluon prop- 
agator in axial gauge, that the coefficient of the 5^ term in Eq. Q5.13 ) cannot be more singular than 



q . This suggests that neglecting the F 2 term, as was done in the studies described above, is a poor 
approximation because there are cancellations in the infrared between Fi and F 2 . 



Singular in the Infrared? Landau gauge studies. There have also been similar studies of the DSE 
for the gluon vacuum polarisation in Landau gauge (Mandelstam, 1979; Bar-Gadda, 1980; Brown and 
Pennington, 1988a, 1988b, 1989). In such approaches the caveat mentioned above is circumvented 
because of the presence of ghost fields, however, this approach has its own problems: 1) There is no 
way to eliminate the diagram involving the 4- gluon vertex as there was in axial gauge. In all of the 
studies to date this contribution has simply been neglected in the hope that the 3-gluon vertex alone will 
contain the essence of the infrared behaviour of QCD if not the details; and 2) The ghost contribution 
is present. Indeed, it is necessary to ensure that the vacuum polarisation is transverse. However, in a 
one-loop, perturbative calculation the ghost-loop only makes a numerically small contribution to the 
vacuum polarisation in Landau gauge. For this reason the diagram is neglected, again in the hope that 
this will not remove the dominant infrared contributions in QCD. 

Subject to these approximations the DSE for the gluon vacuum polarisation contains only the first two 
terms on the right- hand- side of Fig. 2.5. A final step of neglecting the fermion-loop contribution yields, 
as in the axial gauge studies, a closed equation for the vacuum polarisation once the dressed 3-gluon 
vertex is known. 

In a general covariant-gauge the STI for the 3-gluon vertex involves the ghost self-energy, b(q 2 ), and 
the proper ghost-gluon vertex function G^(fc;g,r) = r u G^ u (k;q,r), where k^ is the incoming gluon 
momentum. In Landau gauge the gluon propagator is transverse which entails that, in the limit of 
vanishing ghost momentum, 

G M (*; q, r) D- Xu \r) « D^(r) and G, x (k; q, r) D^(q) » D^(q) (5.22) 

(Marciano and Pagels, 1978, pp. 171-172). This result is commonly used to justify a simplification of the 
STI that is assumed, again, to capture the essence of the infrared behaviour of QCD. A final commonly 
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used assumption is simply to set b(q 2 ) = 0, which reduces the STI to the one given in Eq. (|5.15[) ; i.e., 
to the axial-gauge identity. This is solved by 



(5.23) 



with 

0\^(p,q,r) = 5 V 



Pv 



D u (p 2 ) D n (q 2 ) 



+ 



p2 _ qi 



D u (p 2 ) 



D 



{P^QX-P-Q {Pv - qv 



(5-24) 

where D n (p 2 ) = [1 + U(p 2 )] 1 . As usual, the STI has fixed the vertex up to the addition of a term which 
vanishes when contracted with any of the external momenta and also when any of the external momenta 
become zero. This last fact suggests that the 3-gluon vertex in Eq. (|5.23 ) may alone be sufficient to 
capture the essence of the infrared behaviour of QCD. 

Subject to the approximations discussed above, one now has a single integral equation for the gluon 
vacuum polarisation given by Fig. 2.5 with only the second diagram on the right-hand-side closed by 
using Eq. (|5.23|) for the 3-gluon vertex and this is the equation studied by Brown and Pennington 
(1989). This integral equation has the usual ultraviolet divergences but infrared divergences are also 
possible, especially if D u (p 2 ) ~ p~ 2 . The ultraviolet divergences were handled in the usual way and an 
integral equation, involving the QCD running coupling constant, obtained for the renormalised function 
DR(P 2 ). 

The possibility that has a term of the form p~ 2 was handled by using the "plus" definition in the 
theory of distributions: 



dx (- ) <f>{x,y) 
o \x; 



dx-(<f)(x,y)-6(y-x)<f)(Q,y)) . 



(5.25) 



This enabled Brown and Pennington (1989) to find a solution of the renormalised equation of the form 



p 



— * 

n=l 



P' 



where 



1 + 

d M Tr 



1 + 



P 

7 2 



P 2 +Po 



-28/53^ 



nb 



(5.26) 



(5.27) 



with fj, the renormalisation point, and where A, p , a n and b were determined by requiring self- 
consistency of this from under iteration in the integral equation; i.e., in Landau gauge, subject to 
the [severe?] approximations described above, the gluon propagator has a double pole at small-p 2 . This 
result is in accord with all of the other Landau-gauge studies. 

The effect of fermion loops on the form of the gluon propagator in the infrared has been studied (Brown 
and Pennington, 1988b) using a simpler Ansatz for the 3-gluon vertex, first proposed by Mandelstam 
(1979): 



F 3 

X/iu 



(p,q,r) 



,3(0) 
Xfiu 



(p, q, r) 



D n {p 2 



(5.28) 



where T 3 ^ is the bare 3-gluon vertex. Based on the fact that the results obtained with this Ansatz 
neglecting the fermion loop are qualitatively the same as those obtained using Eq. ( |5.23| ) (Brown and 
Pennington, 1989), this was judged to be a reasonable simplification. This study is not a self-consistent 
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solution of the full coupled DSEs for the fermion and gluon propagators and it has a number of flaws 
introduced by the approximations, however, it is a first step. The results suggest that fermion loops 
act to suppress the infrared singularity in the gluon propagator; i.e., that the double-pole found in the 
pure-gauge sector may be removed by nonperturbative fermion-loop corrections. 



Vanishing at q 2 = 0? In QCD there are seven superficially divergent proper vertices [quark-: T FF ; 
gluon-: Tyv] and ghost-: Tgg] self-energies; quark-gluon-: Tfvf] three-gluon-: Tyvv] four-gluon-: 
IVvvv; and gluon-ghost-: Tgvg] vertices]. The set of DSEs that couple these vertices is the starting 
point for an alternative approach to determining the form of the gluon [and quark] propagator (Stingl, 
1986; Habel et al, 1990a, 1990b; Stingl, 1992). The main difference between this approach and those 
discussed above is in the method employed to solve the DSEs but behind this there is a different 
philosophy. 

The DSEs for the vertices in QCD can be written in the general form 

r = <7 2 $[{r}] (5.29) 

where {T} = {T ff ,T V v,^gg,^fvf,^vvv,^vvvv^gvg}- Standard perturbation theory solves this 
system of equations by iteration, starting with r(°) peri and generating a power series: 

oo 

pipert pi (O)pert _|_ ^""^ ^2rij-ro (O)pert ^ gg^ 

71=1 

where the first iteration is obtained via 

g 2 $ [{ r (0)j»ert|] = g 2 T (l)pert + q^4) (53^ 

The philosophy of this alternative approach is not to abandon the organising principle of perturbation 
theory altogether; the solution of the DSEs is still obtained as a power series in the coupling, but to allow 
each of the vertices at each order to have an essentially nonanalytic dependence on the coupling. In 
employing this approach one assumes a form for the nonperturbative part of the zeroth order vertices, 
Y(o)nonp _ p(o) _ p(o)pert^ j n p rac tice a rational polynomial Ansatz is chosen, and requires that this 
reproduce itself under iteration, which places self-consistency constraints on the coefficients in the 
parametrisation via 

/$[r (0) ] = r {0)nonp + g 2 r (1) + o( g 4 ). (5.32) 

The outcome of this procedure, however, is simply a set of truncated/approximate DSEs which must 
be solved self-consistently. 

Using this approach Habel et al. (1990a) studied approximate coupled DSEs for the ghost self energy, 
b(k 2 ) and proper ghost-gluon vertex function G tMl/ (k;q,r), using rational polynomial Ansatze for these 
functions, the 3-gluon vertex and gluon propagator. A consistent solution was found to be 

b^°\k 2 ) = and G$(k; q, r) = 8^ ; (5.33) 

i.e., no zeroth-order nonperturbative parts is self-consistent at the one- loop level. The reason for this 
is that the loops in the integral equations are either convergent or purely perturbative and thus cannot 
generate a nonanalytic dependence on the coupling. 

This result was used by Habel et al. (1990b) in a study of the coupled DSEs for the gluon propagator 
and 3-gluon vertex. In this study the only diagrams that survived in the DSE for the gluon vacuum 
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polarisation, Fig. 2.5, were the second and fourth on the right-hand-side because all 4-gluon terms are 
of higher order in g 2 and fermions were neglected. The equation for the 3-gluon vertex was reduced to 
one involving itself, the gluon propagator and ghost loops, which are perturbative (Habel et ai, 1990a). 
With a rational polynomial Ansatz for the zeroth-order, nonperturbative elements: the gluon vacuum 
polarisation is assumed to be of the form 

q 2 T^\q 2 ) = c 2 + 2a 2 + -^— (5.34) 

q z + cr 

where, in this particular calculation the simplifying assumption a = = c was made, and the 3- 
gluon vertex, V\^(p,q,r), depends on 9 parameters, each of which multiplies a ratio of p 2 and/or q 2 
and/or r 2 , the coupled DSEs yield 11 equations for the 10 variables. The DSE for the 3-gluon vertex 
yields 9 algebraic, coupled, cubic equations and the DSE for the vacuum polarisation yields two linear 
equations. [The fact that the system is overdetermined is an artifact of the simplifying assumption 
a — = c] One solution is found which, although giving a 3-gluon vertex which is complex, is judged 
to be "reasonable" by the authors on the grounds that the fact that a solution exists at all, given the 
crudity of the assumptions, is very encouraging. This solution yields 

D n ^(q 2 ) 1_ 

q 2 ~ q 2 [l + rK°)(g 2 )] ~ q*T¥ ' 

i.e., a gluon propagator that vanishes at q 2 = 0. 

It has been argued (Stingl, 1986; Habel et al, 1990a, 1990b; Zwanziger, 1991; Stingl, 1992) that 
Eq. (|5.35| ) is a propagator which represents confined gluons because there are no poles on the timelike 
real axis in the complex-g 2 plane and it allows an interpretation of the gluon as an unstable excitation 
which fragments into hadrons before observation [in a time of the order of 1/6]. It is also argued (Habel 
et ai, 1990a, 1990b; Stingl, 1992) that in this framework such a gluon propagator should lead to a 
quark propagator with similar structure in the complex plane, and hence a similar interpretation, but 
this result has not been proven. 

It should be remarked that a possible flaw in these studies is that the Ansatz for the 3-gluon vertex, 
T 3 ( \pi,p 2 ,P3), has kinematic light-cone singularities of the form l/(p 2 ). Such singularities cannot 
arise in perturbation theory (Ball and Chiu, 1980) and hence such an Ansatz cannot reduce to the 
renormalisation group improved 3-gluon vertex in the deep spacelike region. 

It is interesting that the form in Eq. ( |5.35|) is suggested by a number of other studies. It has been 
argued (Zwanziger, 1991) that in order to completely eliminate Gribov copies (Gribov, 1979) and hence 
to fix Landau or Coulomb gauge uniquely in lattice studies, one must introduce new ghost fields into 
QCD in addition to those associated with the Faddeev-Popov determinant in the continuum. Analysing 
the lattice action thus obtained suggests that the gluon propagator vanishes as (g 2 ) 7 , with 7 > not 
determined. Subsequent analysis of a simplified model yields 7 = 1 and, in fact, a gluon propagator 
of the form in Eq. ( |5.35| ) with 6 a finite constant in Landau gauge. Similar considerations lead Gribov 
(1979) to the same result. 



Lattice Simulations. There have been a number of lattice simulations of the gluon propagator (Man- 
dula and Ogilvie, 1987a, 1987b, 1988; Bernard et ai, 1993). All of these studies have fixed a lattice 
Landau gauge, which is plagued by Gribov copies, but did not make use of the modifications suggested 
by Zwanziger (1991). Using 16 3 x 40 and 24 3 x 40 lattices at (3 = 6.0, Bernard et al. (1993) obtained a 
gluon propagator which allowed a fit of the form in Eq. ( |5.35|) at small q 2 but which could not rule out 
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a fit using a standard massive particle propagator. Other lattice sizes and values of (3 were also studied. 
The results at /3 = 6.3 on a lattice of dimension 24 4 were not inconsistent with these results but in this 
case the small physical size of the lattice was a problem. On a lattice of dimension 16 3 x 24 at (3 = 5.7 
it was found that the gluon propagator was best fit with a standard massive vector boson propagator 
with mass ~ 600 MeV. 

These studies represent an improvement in both technique and lattice sizes over earlier lattice studies of 
the gluon propagator (Mandula and Ogilvie, 1987a, 1987b, 1988) but the conclusions are not markedly 
different. Mandula and Ogilvie (1987a, 1987b) using (3 = 5.6,6.0 on a 4 3 x 8 lattice and (3 = 5.8 on a 
4 3 x 10 lattice, obtained results that were consistent with a free massive boson propagator with mass 
~ 600 MeV. 

Viewed as a whole, these studies appear to suggest that the Landau-gauge gluon propagator is finite 
and nonzero at q 2 = 0. This is consistent with an analysis of an approximate DSE for the gluon vacuum 
polarisation using the gauge technique (Cornwall, 1982). 



5.2 Summary 

A number of observations regarding the studies described above are in order. It is clear that the be- 
haviour of the gluon propagator in the infrared is poorly understood with the results obtained depending 
on the approximations/truncations used to obtain tractable DSEs and/or the Ansatze for the propaga- 
tors and vertex functions. From the point of view of a DSE based phenomenological approach to QCD, 
there remains a great deal of freedom in parametrising the infrared behaviour of the gluon propagator. 

If one chooses to believe that those studies which yield q~ A behaviour are the most reliable then a 
caveat must be borne in mind. This structure leads to an infrared divergence in the DSEs for both the 
gluon and quark. It is therefore only defined with respect to some regularisation procedure. One such 
procedure is given in Eq. (|5.25|) but this is not the only one possible; for example, another commonly 
used regularisation is to replace g~ 4 by S 4 (q), which is a distribution that is integrable on any domain 
containing the origin. It is not known whether such a prescription can reproduce itself under an iterative 
solution of the DSE for the gluon vacuum polarisation. The inference that q~ A behaviour corresponds to 
a linearly rising potential in configuration space is also subject to this caveat since the Fourier transform 
is not defined until a regularisation prescription is specified. 

The lattice studies of the gluon propagator are at an early stage; gauge fixing is a difficult problem 
in lattice simulations. However, these studies are of interest because they may be able to support or 
undermine one or another of the approaches to simplifying the gluon DSE. It seems reasonable to doubt, 
however, whether the finite lattice size and associated infrared cutoff will permit the identification of 
an infrared singular gluon propagator. 

Finally, as with QED3 and QED4 , the higher n-point functions have not really been studied at all. 
The studies described above show that they clearly have an important bearing on the structure of the 
quark-quark interaction in the infrared and, therefore, that they deserve further attention. 
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6 Fermion Sector of QCD 



In this section we work primarily in Minkowski space using Minkowski-space conventions unless explic- 
itly specified otherwise. Many elements of the detailed discussions of QED in the preceding sections 
are relevant here and we will draw on them as needed, introducing only those aspects of QCD which 
differentiate it from QED in order to proceed. 



6.1 Dynamical Chiral Symmetry Breaking 

The starting point for studies of DCSB in QCD is the DSE for the quark self-energy given in Eqs. ( |2.68| ) 
and ( p.69| ). We will discuss various models and approximation schemes, which attempt to incorporate, 
as far as possible, the symmetries and leading-log perturbative QCD results discussed in the previous 
section. Of course, it is of great interest to see how confinement can arise in such a scheme and we 
discuss this in Sec. |6~2. 



Proper Quark-Gluon Vertex. A common first approximation is to neglect the effect of ghosts in the 
quark sector and assume that this can be compensated for by fine-tuning the phenomenology in the 
gluon sector. The justification of this, such as it is, and its inherent flaws is described in Sec. |5.1| . It is 
undesirable and raises the question of how important the resulting violation of QCD gauge invariance 
is. Indeed, it can be argued that this approximation may compromise confinement. Nevertheless, it 
does allow good use to be made of what has already been learnt form QED studies. The neglect of 
ghosts in the quark sector implies a theory of quarks with some Abelian-like characteristics, with QCD 
gauge invariance reducing to the requirement that the colour-current be conserved at the quark-gluon 
vertex. For this purpose it is sufficient to require that Eq. ( |2.5S| ) [without ghosts] is satisfied for the 
renormalised quantities, which has then the same form as Eq. ( |2.30| ). [The discussion associated with 
Eq. ( |5.23| ) is relevant here.] Making use of the QED studies discussed previously [e.g., Ball and Chiu 
(1980), Curtis and Pennington (1990, 1991, 1992), and Burden and Roberts (1993)], we can write then 



r M (p' ,p) = r^ c (p' ,p) + Tqp(p',p) + transverse parts = 7 M + A M (p',p) 



(6.1) 



where we have indicated that only the longitudinal behaviour is specified by Eq. ( B.l|) . The two specified 
parts of the vertex are the so-called Ball-Chiu vertex r B c and the transverse addition suggested by 
Curtis and Pennington on the basis of requirements of multiplicative renormalisability. The transverse 
pieces are proportional to [g ua — [k v k u /k 2 )] where k v = (p' — p) v and so we have then k fJL T^ l {p\p) = 
kfj,T^ c {p' ,p). These are defined by 



with 



A(p 2 ) + A(p' 2 ) (p + pT 

r 7^ + 



Y(p'-p 2 )-(p' + p) 1 



A(p 2 ) - A(p' 2 ) 
A(p' 2 ) - A(p 2 ] 



4+ f) 



d(p',p) 



2d(p',p) 

(p' 2 - p 2 ) 2 + [M{p 12 ) 2 + M(p 2 ) 2 ] 2 
(p 12 + p 2 ) 



B(p 2 )-B(p' 2 )\j,(Q.2) 
(6.3) 
(6.4) 



where M(p 2 ) = B(p 2 )/A(p 2 ). The QED vertex proposed by Curtis and Pennington has the unspecified 
transverse parts in Eq. Q6.1]) set to zero; i.e., it is defined as 



rUp',p) = ^ BC (p',p) + r T CP ! (p', P ) 



(6.5) 
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It is straightforward to verify that (j>' ~V)^bc = S -1 ^') ~ S" 1 ^) and that (p' — p)/jTq P = as required 
by Eq. (fT30[ ). Since Eq. p37| ) implies that A(-fi 2 ) ~ 1, then in the limit that p' 2 ~ p 2 ~ -/j 2 the RHS 
of Eq. (|6.1| ) becomes 7 M plus terms proportional to dA/dp 2 and dB/dp 2 respectively, which for large 
fi 2 will be very small. Thus in this limit we have r CT ~ 7 CT which is the usual perturbative vertex. The 
only way of determining the nonperturbative transverse pieces and ghost contributions in Eq. (|6.1| ) is 
to attempt to simultaneously solve the DSE for the quark-gluon vertex. While Tcp does not include 
ghost effects it does give the correct perturbative limit, reduces to the free vertex in the absence of 
interactions, is free of unphysical kinematic singularities, transforms correctly under parity and charge 
conjugation, is Lorentz covariant, and preserves multiplicative renormalisability. 

It is unfortunate that there are difficulties in using ghost-free gauges, such as the axial gauge. Any 
approximation made in a noncovariant gauge tends to destroy the covariance of physical observables. 
In addition, in axial gauge additional singularities of the form n ■ p enter into propagators. The biggest 
difficulty, however, is that since one now has two four-vectors [p u and n u ] the most general Lorentz 
structure of the inverse quark propagator is S 1-1 = (Ap 1 — B) + [Cp 1 — D)fi, where the functions A, B, C, 
and D depend on each of the Lorentz scalars p 2 , p ■ n, and n 2 = ±1. The fact that the scalar functions 
A, B, C, and D can depend explicitly on n ■ p and n 2 significantly complicates the DSE analysis. For 
these reasons studies of DSEs where Lorentz covariance is important are most frequently carried out in 
a covariant gauge. 



Asymptotic Quark Self-Energy Behaviour. Recall that the renormalised DSEs relate quantities at 
the renormalisation point fi and that the appearance of the regularisation parameter A is through the 
renormalisation constants Z s (A,fi), Z D (A,fi), Z r (A,/z), Z m (A, fi), Z g (A,fi) and hence also through the 
bare quantities: the bare coupling constant go (A), gauge parameter £o(A), and quark mass m (A). So 
for the quark self-energy in Eqs. ( |2.68| ) and ( |2.69| ) we understand that 



-iTf(p) = lz r (A )A *) g 2 f '^(^(iSmiDnP-mm^)) ■ (6-6) 

In practice, numerical calculations of DSE solutions are carried out in Euclidean space where it is 
convenient to take the regularisation parameter A to be a momentum cut-off, Auv- In this case, it is 
further convenient to choose /z=Auv since then all of the Z-factors are unity; i.e., Zs(Auv, Auv) = 1, etc. 
Thus for Auv an d fi sufficiently large we can use the known asymptotic behaviour to obtain information 
about quantities for fi ^ Auv- It is of course then a standard check on numerical solutions that physical 
results so obtained at fixed renormalisation point fi are independent of Auv f° r Auv — * 00 • Discussions 
of these issues can be found in Fomin et al. (1983) and references therein. An added advantage of 
choosing /i = Auv with Auv ^> Aq C d is that S'(p) ~ at p 2 = —fi 2 and so this can be essentially 
neglected, which gives S(p) — S'(j»). Using the above arguments we have arrived at the commonly used 
expression for the quark self-energy DSE [with fi = Auv] 

4 r A vv d 4 k n 2 
E(p) = • - An J — 7ct S(k)<L- D™(k - p)T u (k,p) . (6.7) 



The asymptotic behaviour of the quark-quark scattering kernel (K) is known from renormalisation 
group analysis and to leading order [in Landau gauge] is [Q 2 = — (q' — q) 2 } 

Knm;n>m>(q', q; P) ^ -* 9 2 (Q 2 hn m D 0au(q - j)^, (6.8) 

where to this order K is independent of the center-of-mass momentum, P, and the renormalisation point, 
fi, and where n,m,n', and m' are spinor indices, [see, e.g., Nakanishi (1969)]. D is the perturbative 
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gluon propagator; i.e., D with II = 0. In the asymptotic region K, as expected, is dominated by one- 
gluon exchange (K ~ — i g 2 TDT) which leads to Eq. Q5.8| ). Since in this region T ~ 7 up to logarithmic 
corrections [which are unimportant for generating the leading behaviour of the self-energy] , we can use 
Eq. QS.8| ) for g 2/ ~fDT in Eq. ( |2.68| ). It is relatively straightforward to verify that this produces the correct 



asymptotic behaviour for the running mass [for a detailed discussion see, e.g., Fomin et al. (1983)]. 

As we showed in Sec. |5], if one makes the assumption Zf h = Zf then one may write a s (Q 2 )D^ vo (q) 
in place of the transverse part of (g 2 /47r)D 1/ ' T (q) in Eq. (|6.9| ) and so, in this approximation, the non- 
perturbative part of the gluon propagator is absorbed into the nonperturbative structure of a s {Q 2 ). In 
addition, it can argued that [at least part] of the ghost effects omitted from r CP and/or r B c can be ab- 
sorbed into the nonperturbative behaviour of a s (Q 2 ). This then gives rise to the semi-phenomenological 
expression for the quark self-energy [in Landau gauge] 

4 Muv d A k 

E(p) = i - 4vr J — - la S(k)a s (-(k - p) 2 ) D™(h - p)T v (k, p) , (6.9) 

where T u is then to be chosen consistent with the WTI and where nonperturbative gluon and ghost effects 
have been [as far as possible] absorbed into a s . Hence, the problem has been reduced to motivating 
particular nonperturbative gluon and ghost modifications and exploring the resulting forms for the 
nonperturbative behaviour of the function a s . Equation ( |6.9| ) is a form of the quark self-energy DSE 
in Landau gauge that we will consider in some detail here. The advantages of Landau gauge are that 
the gluon propagator is purely transverse and that the gauge parameter does not run in Landau gauge 
£,(Q 2 ) = 0, [£ = is an ultraviolet fixed point, albeit not a stable one in general. See Marciano and 
Pagels (1978) and Pascual and Tarrach (1984)]. 



Quark Dyson- Schwinger Equation. Using Eq. 
Euclidean space in Landau gauge [£ = 0]: 

16tt f A uv d 4 k a s ((p-k) 2 ) 



A(p 2 



1 + 



A(k 2 



in Eq. 
1 



and using Y h 



we obtain in 



(2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 



A(k 2 ) + A(p 2 ) 1 
' 2 p~ 2 



[3p • k — h(p, k)] 



-A(k 2 )AA(k 2 ,p 2 " 

WAuv + 

W{k 2 ) 



B{p 2 ) 



k - 



(k ■ p) 2 k ■ p 



p^ 

16tt f A uv d 4 k a s ((p-k) 2 ) 



-h(p,k) 



B(k 2 )AB(k 2 ,p 2 )^^- 



(6.10) 



1 



(2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 



A(k 2 ) + A(p 2 ) 



+ 



B(k 2 )AA(k 2 ,p 2 ) - A(k 2 )AB(k 2 ,p 2 ) h(p, k) 



(6.11) 



where h(p, k) = 2 [k 2 p 2 — (k ■ p) 2 } /(k—p) 2 and AF(k,p) = \F{k 2 ) — F(p 2 )}/[k 2 — p 2 ] and where mA uv is 
the running quark mass when \i = Auv- Calculations to be reported later typically have no ECSB [i.e., 
^a uv = 0]. For a calculation that explicitly includes ECSB quark masses see Williams et al. (1991). 
Including the additional Curtis- Pennington term in the vertex, [i.e., using T u = T^p] means that these 
equations are modified as follows: 

A{p 2 ) = RHS of Q610D 

16tt /- a uv d 4 k a s {{p-k) 2 ) A{k 2 )AA{k 2 ,p 2 ) (k 2 - p 2 ) 3(k 2 - p 2 )k ■ p 



+ 



(2tt) 4 (p-k) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 2d{k,p) 



P 



(6.12) 
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b(p 2 ) = rhs of (fnn 



167T Muv tfk a 8 ((p-fc) 2 ) £(fc 2 )AA(fc 2 ,p 2 ) (fc 2 -p 2 ) 
3 J (2tt) 4 (p-A;) 2 A 2 (k 2 )k 2 + B 2 (k 2 ) 2d(k,p) 1 PJ ' l ' J 

These Euclidean space equations are obtained by a simple transcription from the corresponding Minkow- 
ski space equations as discussed in Sec. In view of those discussions of rotations to Euclidean space, 
it is understood that the semi-phenomenological nonperturbative part of a s is a Euclidean space Ansatz. 
The results obtained are independent of the ultraviolet cut-off [Auv] provided it is chosen sufficiently 
large since the integrals in Eqs. ( |6.10 ) and ( 6.11|) are convergent. Note that this would not be the case 



were we to take a s (Q 2 ) constant. If we were to make the the approximation a s ((p — k) 2 ) ~ ck s (p>) where 
p 2 = max(p 2 , k 2 ), then we would find, after an angle integration, that A(p 2 ) = 1. This unnecessary but 
simplifying approximation [already discussed in Eq. fl4.7p ] has sometimes been used in the literature 
and we refer to it as the angle approximation, [see, e.g., Fomin et al. (1983) and references therein]. 



Chiral Symmetry Considerations. When there is no ECSB renormalised quark mass [m = 0], we 
have exact chiral symmetry and conservation of the axial- vector current leads to [in Landau gauge] 



M(-Q 2 



Q 1 



HQ 2 /k 2 QCD ) 



d M -l 



(6.14) 



where c is some constant independent of fi. This is in contrast with the form of the running mass given 
in Eq. ( |2.62f) . Then the asymptotic form for the quark mass can be written in convenient shorthand as 



M(-Q 2 



c 

~Q 2 



ln(Q 2 /A| CD 



+ m 



HQV^qcd) 



(6.15) 



where again the /i-dependence is not explicitly indicated. Note that Eq. (|6.15| ) is to be understood in 
the sense that for exact chiral symmetry the second term on the RHS is zero and the first term is the 
dominant one, while in the presence of ECSB [i.e., m/0] the second term is the dominant asymptotic 
behaviour. With ECSB there will be many terms which are suppressed by powers of ln(Q 2 ) with respect 
to the second term but which dominate the first as Q 2 — > oo. The [renormalised] quark condensate (qq) 
is a measure of the degree of DCSB and can be defined as (qq) = (vac| : <?(0)g(0) : |vac), where |vac) 
refers to the nonperturbative vacuum and the normal ordering of the [renormalised] operators is with 
respect to the perturbative vacuum. Denoting S^ 1 (p) = Z^ 1 (p 2 )[p' — M(p 2 )] as the nonperturbative 
quark propagator we have, in the absence of ECSB, 



<gg>=-lim tr {S(x, 0)} = -12i 

x— >0+ 



A uv d A p Z{p 2 )M{p 2 

(2vr)V -M 2 (p 2 



-Yli 



A uv d A p 



B{p 2 



(2n) 4 A 2 (p 2 )p 2 - B 2 (p 2 ) ' 
(6.16) 

where S(x, y) is the coordinate space quark propagator and where the trace over spinor [4] and colour 
[3] indices gives the factor 12. In Appendix |B| we show that c, as defined in Eq. (|6.14j) , satisfies, in the 
limit of exact chiral symmetry [m = 0] in Landau gauge, 



A-n 2 d 



M 



<qq> 



(6.17) 



which since c is a constant independent of the renormalisation point fi also implies that < qq >~ 
ln(/i 2 /Ag CD ) M . Eq. ( |6.14| ) is a result also obtained in discussions of the operator product expansion 
OPE] and QCD sum rules [see, e.g., Politzer (1976,1982), Gasser and Leutwyler (1982), and Reinders 
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Fig. 6.1. Shown are (a) the proper quark-pseudoscalar vertex coupling 
to a pion, (b) the proper quark-axial-vector vertex, and (c) the pion 
decay component of the quark axial-vector vertex which dominates in 
the soft-pion, chiral limit. 



et al. (1985)]. In the OPE the quark condensate is taken to be a coefficient with the above /^-dependence 
irrespective of whether or not ECSB is present. From the discussions of Appendix [B] it is clear that the 
OPE quark condensate only corresponds to the "natural" definition of Eq. ( 6.16Q in the limit of exact 
chiral symmetry. Thus in the limit of exact chiral symmetry there are two ways of obtaining the quark 
condensate: 1) directly from the integral of Eq. ( 6.16j ) and 2) from the asymptotic behaviour of the 
quark self-mass as follows from Eq. ( |6.15|) and Eq. (|6.17|) . 



In the limit of exact chiral symmetry we have the WTI of Eq. ( 2.59| ) which gives 



W(p',p) 



k-*0 



r l5 B(p 2 ) . 



r r 



A(p' 2 ) p" - A(p 2 ) p 1 - (B(p' 2 ) + B(p 



75 



(6.1* 



Note that the usual perturbative axial- vector vertex is (r/2) 7^75. 



In the chiral limit as k — > the axial- vector quark proper vertex [see Fig. 6.1(b)] becomes completely 
dominated by the pseudoscalar coupling of a massless pion to the quark [r 5 in Fig. 6.1(a)] and the 
subsequent weak decay of the pion into an axial- vector current, which is depicted in Fig. 6.1(c). The 
pion decay constant f n ~ 93MeV is defined by the axial- vector transition amplitude for an on-shell pion 
[k 2 = m 2 and m, n are isospin indices] 



(0|i4r(0)|7r»(A;)> = 
From Fig. 6.1(c) and Eq. (|6.19|) we find then that as k — > 



ifl{p\p) —> [XT(j/,p)] 



z(0\A°\n m (-k)) 



r 5 (p',p)fAk a /k 2 



(6.19) 



(6.20) 



Acting on both sides of Eq. ( |6.2(J| ) with k a , we find the that k a T^(p',p) if n T 5 (p',p) as k — > 0, which 
when combined with Eq. ( |6.18| ) gives the Goldberger-Treiman relation for the quark-pseudoscalar vertex: 



f 5 (p,p) = irj 5 



B{p 2 



ZT75- 



Z-\p 2 )M(p 2 



(6.21) 



66 





i(xl2)f% 



Fig. 6.2. The integral equation for the pion decay constant (/„-). 



The simplest generalisation of Eq. ( |6.21| ) to k ^ is perhaps 



f 5 (p',p) = i r l5 -[B(p' 2 ) + B(p 2 )]/f 7T 



(6.22) 



However, to properly determine the pseudoscalar vertex away from the soft-pion limit obviously requires 
solving the pseudoscalar vertex DSE. This is equivalent to the solution of the pion BSE, where the pion 
Bethe-Salpeter amplitude is defined as 



X ^(P,P) = I dSe^-™]-(0|T^(0)^(z)|7r"(P)) 



(6.23) 



where \ir n (P)) denotes a pion state with momentum P, n is the isospin label, and x is related to r§ by 



x(p, P) 



P P 

* F 5 [P+^,P-^ 



(6.24) 



In Eq. ( p.23| ) all colour indices are suppressed. In the amplitude x the incoming pion has 4-momentum 
P, the quark has 4-momentum p + (P/2), and the antiquark has 4-momentum p — {P/2). The parity 
and time-reversal invariance properties of a pseudoscalar imply that we can write 



Xip, p ) = 15Xp(p, P) + PJ5Xa(p, P) + ThsXA'Vip, P) + \jf, P]J5Xt(p, P) ■ 



(6.25) 



Here Xi{Pi P) f° r * — Pi A, A', and T are scalar functions of p 2 , (p- P), and P 2 . The subscripts P, A, and 
T denote the spinor matrix structure of the particular component of X- The three functions XPi Xa, 
and xt are even functions of p ■ P, whereas xa> must be an odd function of p ■ P. A similar expansion 
can also be made for the proper quark-pseudoscalar vertex IV In the chiral limit where the pion is 
massless we obviously have for an on-shell pion P 2 = m 2 = 0. 

Using the above it is possible to obtain an expression for f 2 from the integral equation for pion decay 
illustrated in Fig. |6.2| . Our discussion here is a relatively straightforward extension of the arguments 
presented by Pagels and Stokar (1979,1980) and Cornwall (1980) to the case where A(p 2 ) ^ 1. From 
Fig. IO we find that 



i{o\Ar(o)\K n (k)) 



i(ik a f n S mn ) = -k°fj 
d A p 



(2vr) 4 
d A p 



trj [iS(p + k)} \iTf(p + k, p)] [iS(p)) [*(t72) 7 ct 75] 
tr\\x?(k,p))[i(T n /2h*7 5 ) 



(6.26) 
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where, as usual, the factor [-1] on the RHS arises from the fermion [i.e., quark] loop. In order to avoid 
double counting we need to take the axial-vector vertex as a perturbative vertex for the same reason 
as this is done in the photon self-energy diagram in Fig. |2.1| . [The argument for f n in Williams et 
al. (1991) is erroneous for this reason, although the error actually has very little effect on the numerical 
results reported therein.] If, in the limit k — > in Eq. ( |6.26|) , we equate coefficients of k a , using the 
approximation: T 5 (p',p) « T 5 (p,p) + 0{k)^^ i.e., assuming the Dirac structure of the O(k) correction 
is purely 75, and use the Goldberger-Treiman relation Eq. ( |6.21| ) for T^, we obtain, in the limit of exact 
chiral symmetry, [see Appendix |C| for details] 



d 4 p Z(p 2 )M(p 2 



(2vr) 4 [p 2 - M 2 (p 2 )] 2 



M(p 2 ) - 



p 2 dM 
~2djP 



(6.27) 



which as a physical observable must be independent of the choice of renormalisation point, as can be 
verified numerically for solutions. It should be emphasised that, while this derivation is quite general, 
the assumed form of T§(p',p) for p' 7^ p influences the result, which is therefore model dependent. 
[Another form appears in Eq. ( |7.57|) . Both forms reduce to the Pagels and Stokar (1979) result when 
Z(p 2 ) = A~ l (p 2 ) = 1 which indicates that the O(k) corrections to T 5 (p',p) measure deviations from 
A(p 2 ) = 1, consistent with the result of Delbourgo and Scadron (1979).] It is straightforward to write 
down the Euclidean space versions of Eq. ( |6.16| ) and Eq. fl6.27p by inspection and for completeness these 
have been given in Appendices [B| and Q respectively. 



Numerical Results. As a specific example of the above arguments for studying DCSB in QCD we 
describe numerical results from one model DSE for the quark propagator (Hawes et al. , 1993). [Other 
models and results will be compared and summarised below.] This study uses the Ball-Chiu and Curtis- 
Pennington vertices together with an Ansatz for the gluon propagator in Landau gauge of the form 

(g 2 /4n)D^(q) -> a s {Q 2 )D^{q) (6.28) 

with 

« s (g 2 )=a s (r;Q 2 )-^- and a s (r-Q 2 ) = — ^ , (6.29) 

{Q ) +0 In [r + (Q /AqcdJ 

where lnr > [i.e., r > 1] and b are real constants to be determined. This is to be understood in terms 
of the arguments which lead to Eq. (|6.9|); i.e., it is understood that the ghost and nonperturbative 
gluon propagator behaviour have been absorbed into the nonperturbative structure of a s {Q 2 ). It is 
clear that this Ansatz has the asymptotic form of Eq. (|2.56|) as Q 2 — > 00. In the infrared, we see that 
the effective gluon propagator vanishes as Q 2 — *> 0; i.e., D(q) ~ q 2 /[{q 2 ) 2 + b ]. [The origins of this form 



of gluon propagator are discussed in detail in Sec. pTT] .] A closely related study, where the infrared gluon 



propagator has a Gaussian form and where the Ball-Chiu vertex was used, is described in Hawes and 
Williams (1991). 

We see that the parameter r > 1 plays the dual role of regulating the infrared behaviour of the 
logarithmic term and of determining the strength of the nonperturbative infrared interaction. Using 
the quark condensate as an order parameter, we find that there is a critical value of b = b c such that 
the model does not support dynamical chiral symmetry breaking for b > b c . 



In Sec. |6.2| we discuss confinement and illustrate a confinement test by applying it to this model. The 
results suggest that, for all values of b, the quark propagator in the model is not confining. Two 
properties of a particle's propagator which together are at least sufficient to ensure confinement of the 
particle are: 1) the absence of a Kallen-Lehmann representation and 2) no singularity on the timelike 



68 



q 2 axis. This will be discussed in more detail in Sec. |6.2| . The gluon propagator implied by Eq. ( |6.29| ) 



satisfies both of these requirements. Although it can be argued that it describes a confined gluon, it may 
seem that quark confinement is counterintuitive, since it would appear to provide a weak interaction 
between quarks at small p 2 , corresponding to large distances. Nevertheless, as discussed in Sec. pTT , 
quark confinement is not implausible in this model. 

We note too that there has been an attempt to employ a gluon propagator similar to Eq. ( |6.29| ) in a 
study of quarkonium spectra by Becker et al. (1991). A Blankenbecler-Sugar reduction of a ladder-like 
approximation to the BSE is used and it is argued that, in the bound state equation that results, one 
can approximate the effect of Eq. ( |6.29| ) by a Coulomb potential for all r. 

We can now study the implications of Eq. ( |6. 9[ ) and Eq. (|6.29| ) for the structure of the quark propagator. 
A similar study has been carried out by Alkofer and Bender (1993). Using Landau gauge [£ = 0] since, 
as discussed earlier, the form of gluon propagator we consider here is most realistic in this gauge, we 
find that in both the Ball-Chiu and Curtis-Pennington cases there are regions of DCSB and unbroken 
chiral symmetry characterised by a two-dimensional phase diagram in (6 2 ,lnr) space, where h and r 
appear in Eq. ( |6.29|) . The phase transition is second order. 



For the most part in this study r^ p (/c,p) was neglected; i.e., a minimal Ball-Chiu Ansatz was used. 
This has the virtue of simplifying the integral equations. This term was included for a single value 
of r and a number of values of b and found to generate a small quantitative change in some of the 
characteristic quantities calculated in the model but not to alter its qualitative features. The integral 
equations were solved numerically by iteration on a logarithmic grid of x = p 2 /Aq CD and y = A; 2 /Aq CD 
points. This ensured that the results were independent of the seed-solution and grid choice. The results 
were independent of the UV cutoff, which was Ayy = 5 x 10 8 Aq CD , and this value was also sufficient 
to ensure that the leading- log behaviour of the mass function had become evident. 

We are interested in determining whether the model gluon propagator specified by Eq. (|6.29|) can 
support DCSB - a crucial feature of QCD. The quark condensate, which is gauge- invariant, was studied 
earlier and is an order parameter for DCSB. In cases for which our iterative solution procedure for the 
DSE converged quickly, with relative errors of less than 1 x 1CT 6 , the condensate could be obtained 
easily. However, for values of b 2 near a phase transition the convergence could be extremely slow. In 
those cases, the numerical solution was examined at constant intervals through the run [say, every 50th 
cycle], and the condensate evaluated in each case. Aitken extrapolation (Wimp, 1981) was then used 
to find the "infinite-cycles" limit. In several cases the program was subsequently run until the solutions 
had converged to within 1 x 10~ 6 and the extrapolated result always matched the actual result to within 
a few parts in 10~ 6 . 

We solved Eqs. ( |6.10| ) and ( |6.11| ) for values of lnr in the domain [0.0, 0.7] and b 2 in [0.1, 1.0] using the 
minimal Ball-Chiu vertex and we plot the condensate obtained from our solutions in Fig. 6.3. This figure 
shows regions of unbroken and dynamically broken chiral symmetry. Our numerical results suggest that 
the condensate rises continuously from the transition boundary and hence that the transition is second 
order. As a consequence we assumed that the order parameter, (qq)^, behaves as 

(qq)^z)^ C (l-^y (6.30) 

for z — * z~ [for z equal to either lnr or b 2 ] and extracted the critical points, z c , and critical exponents, 
/3, using ratio-of-logs methods. We list these quantities in Table [□] and, in Fig. 6.4, plot the critical 
curve in the (6 2 ,lnr) plane. 
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Criticality plot, -<qq> /J vs. b , ln(x) 
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Fig. 6.3. Criticality plot for (—(qq)^)^ as a function of lnr and b 2 . The 



condensate, (— {qq)^) 3 , is in units of MeV, scaled to /i 2 = 
is in units Aq CD ; the gluon regulator r is dimensionless. 



lGeV, and b 2 



We also solved Eqs. ( |6.10|) and ( p.ll| ) with the Curtis- Pennington additions, Eqs. ( |6.12| ) and (|6.13|) , 
using lnr = 0.6. The critical curve [in b 2 ] in this case is shown in Fig. 6.5 along with the minimal 
Ball-Chiu results for the same value of lnr. The effect of the Curtis- Pennington addition is to lower 
the critical value of b 2 but, as we show below, the critical exponent is unchanged. This curve illustrates 
the point that the qualitative features of the model are not affected by this modification of the model 
quark-gluon vertex. 



From Table 6.1 we find: 



(3 BC = 0.575 and a p = 0.024 . 



(6.31; 



We note that the critical exponent obtained with lnr = is quite different from the others. This is a 
special case since for this value the propagator does not vanish in the infrared: 



f D (Q 2 = 0) = d M rr^ 
Air o 4 



If we neglect this point in our analysis then we find 

(3 BC = 0.572 and a p = 0.020 . 



(6.32) 



(6.33) 



The results in Eqs. (|6.31| ) and (|6.33| ) are in agreement with those of Roberts and McKellar (1990) where 
it is argued that f3 = 0.589 ± 0.031. That study used b 2 = and found a critical value of lnr = 1.69 
which complements the results reported here. 
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Table 6.1. The critical points and exponents extracted for various values 
of In r; the cumulative result is /?bc = 0.575, with erg = 0.024; excluding 
the point with lnr = 0.0, f3 BC = 0.572 with g p = 0.020. 



In t 


6q - Critical b 2 value 


j3 - Critical Exponent 


<jp - standard deviation in (3 


0.00 


0.6439 


0.609 


0.030 


0.10 


0.5448 


0.579 


0.021 


0.20 


0.4642 


0.570 


0.021 


0.30 


0.3932 


0.573 


0.021 


0.40 


0.3289 


0.567 


0.021 


0.50 


0.2706 


0.567 


0.021 


0.60 


0.2180 


0.561 


0.021 


0.70 


0.1710 


0.579 


0.021 



We also calculated the critical exponent using our numerical DSE solutions obtained with the Curtis- 
Pennington addition to the vertex at lnr = 0.6: 

p CP = 0.579 , ap = 0.015 . (6.34) 

This suggests that the vertex modification does not alter the critical exponent of the transition; a 
conclusion that is also supported by the observation that the vertex used in Roberts and McKellar (1990) 
was not of either of the above forms but was, effectively, a simple modified rainbow approximation: 
T ll (k,p) = A(k 2 h ll . 



Other Studies and Results. There have been numerous studies of DCSB in QCD using the quark 
DSE and it is unfortunately not possible to give an exhaustive summary of all of them herein. However, 
all of the studies are well represented by the above exemplary study: each necessarily assumes or argues 
for one or another form of effective gluon propagator and makes some assumption about the quark-gluon 
vertex. In addition, many of the earlier studies used the approximation of Eq. (|4.7|) in order to obtain 
an approximate differential equation for the quark propagator. Many also used rainbow approximation, 
Eq. ( |3.2Cp , in Landau gauge. 



There are some interesting early studies by Higashijima (1983, 1984) in Landau gauge and using rainbow 
approximation with the following Ansatz for the running coupling: 

a s (Q 2 <Ql)= , dM7T r and a s (Q 2 > Ql) = . ^ v . (6.35) 

° ln(Q^ CD ) ° ln(Q 2 /A^ CD ) 

Employing Eq. ( |4.7| ) in order to obtain a differential equation it was shown (Higashijima, 1984) that for 
lnfQo/AqcD] < 0.88, corresponding to a s (0) ~ 1, the model admits DCSB. Importantly, it was shown 
that in order to obtain a differential equation for the mass function whose solution agreed with the 
one-loop renormalisation group result, Eq. (|6.14 ), one must not neglect derivatives of a s (Q 2 ). 



Subsequently there were a series of studies by Atkinson and Johnson, (1987, 1988a, 1988b). It was 
found that in Landau gauge, with 

a s (Q 2 ) = dM7T - N , (6.36) 
In ( r + jf- 
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Critical curve for propagator with Ball-Chiu vertex 
2.0 I . 1 . 1 . . . 



*Endpoint: b 2 =0, ln(T)=1.69 
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Fig. 6.4. Critical curve for the phase transition in the (lnr, b 2 ) plane. 
The asterisk is the result extracted from Roberts and McKellar (1990). 

with a vertex Ansatz which satisfied the WTI, but which involved kinematic singularities, and using 
Eq. ( f4.7| ), one obtains a differential equation for the mass quark mass function whose solution is such 
that the model manifests DCSB for a s (0) > a£(0) ~ 0.9 (Atkinson and Johnson, 1988b). 

This model was reanalysed by Roberts and McKellar (1990) who solved the integral equations for A(p 2 ) 
and B(p 2 ) directly and analysed the reliability of Eq. (|4.7|) . It was found that using Eq. (17) the model 



manifests DCSB for a s (0) > a^(0) ~ 0.78. The difference between this result and that of Atkinson and 
Johnson (1988b) being due to a more detailed analysis by Roberts and McKellar (1990) of the behaviour 
of the mass function near the critical coupling. Upon solving the integral directly it was found that 



the critical coupling increased to a c s (0) ~ 0.89: a 12% rise. The approximation of Eq. ( f4.7|) is thus 
seen to be qualitatively reliable in this case. This is the best possible application of the approximation, 
however. It was shown that Eq. (|4.7| ) becomes unreliable as a s (0) is increased making it suspect in 
any study in which the model gluon propagator is large or singular at Q 2 = 0. An important result of 
this study was that, even in Landau gauge, A{p 2 ) ^ 1 when the integral equations are solved directly 
without further approximation. [The same it true in the QED 3 studies described in Sec. |].] 

Equation (|4.7|) was used by Krein et al. (1988, 1990) in a study of a 1/C} 4 infrared form for the effective 
gluon propagator, which lead to a confined quark. However, these results were later seen to be completely 
dependent on the use of the angle approximation to regulate the divergent infrared gluon behaviour. 
This emphasises the limitations of this approximation discussed by Roberts and McKellar (1990). 

The difficulties associated with using a model gluon propagator with a nonintegrable singularity in the 
infrared in a covariant quark-DSE were are also illustrated in a study by von Smekal et al. (1991). In 
this study 

«s(Q 2 ) = 7 ^ x + d M 7if(Q 2 ) , (6.37) 
In + 

\ QCD 
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Quark Condensate Results for Two different Proper Vertices 




1 



Fig. 6.5. Comparison of the (—(qq)^)^ condensate curves for the min- 
imal Ball-Chiu and Curtis-Pennington Ansiitze for the proper vertex. 
Both curves have lnr = 0.6. Diamonds, o, connected with solid lines are 
the results from the B-C vertex; plus-signs, +, connected with dashed 
lines are results from the C-P vertex. 

with f(Q 2 ) regular as Q 2 —>■ 0, and an Ansatz for the quark-gluon vertex was chosen such that it 
satisfied the Ward identity. This was subsequently approximated in order to simplify the analysis. 
As in all cases with a Q~ A singularity in the model gluon propagator a regularisation procedure was 
necessary and these authors chose to make the replacement Q 2 — > [Q 2 + fi 2 } in the singular part of the 
running coupling. [A value of /i = 10~ 4 Aq CD was found to reduce the sensitivity of the results to this 
parameter to less than 1%. This might be compared with the prescription employed by Brown and 
Pennington (1989), Eq. ( |5.25| ).1 The coupled, nonlinear integral equations for A(p 2 ) and B(p 2 ) obtained 



with / = were solved and the model shown to manifest DCSB. However, this version yielded a rather 
poor phenomenology in the sense that the quark condensates and pion decay constant, f ni calculated 
from the solution for the quark propagator, did not agree well with the accepted values. It was found 
that this could be remedied by using f(Q 2 ) = 7 2 /[Q 2 + m 2 } with 7 = 2.5Aq C d, m g = 0.2Aq C d and 
Aqcd — 0.5 GeV. [The value of the effective gluon mass implied by this parameter choice is smaller 



than the value of ~ 0.5 GeV suggested by other studies; for example, those described in Sec. |0 



In many respects this study is qualitatively and quantitatively similar to that of Williams et al. (1991' 
in which the Q~ A singularity was regularised by replacing it with an integrable distribution: 



a s (Q 2 ) = C A^ CD Q 2 5\Q) + ^ . (6.38) 

ln(r + # ' 



QCD 



[The first term alone was studied in rainbow approximation by Munczek and Nemirovsky (1983) and 
Cahill and Roberts (1985).] Williams et al. (1991) used Landau gauge, a simplified vertex: T u = 
(l/2)[A(p' 2 ) + A(p 2 )]j u , and solved the coupled, nonlinear integral equations for A(p 2 ) and B{p 2 ). It 
was found that once f n was fitted by adjusting the infrared coupling strength, C, then the quark 
condensate assumed a value in the typical range, {-{qq)) 1/3 ~ 225 ± 25MeV at y, = lGeV. A similar 
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study with the Ball-Chiu vertex and a Gaussian form for the gluon propagator in the infrared reached 
similar conclusions (Hawes and Williams, 1992). 



An effective long-distance quark-quark interaction was derived from dual QCD, which is a long-distance 
treatment of QCD where the vacuum has the properties of a dual superconductor. The roles of the 
electric and magnetic fields are interchanged in this picture [for a detailed review see Baker et al. (1991)]. 
This effective quark interaction was then used as input into studies of DCSB [see Baker et al. (1988) and 
Krein and Williams (1991)] with some success. It is unfortunate, that a sign error was recently detected 
in the original derivation by S. Kamizawa [Baker et al. (1988, erratum 1991)], which leaves open the 
question of whether or not DCSB occurs naturally in dual QCD. Studies of DCSB have also been 
generalised to a chromodielectric model of QCD where a mean field treatment of the colour-dielectric 
constant and DCSB arguments leads to a chiral, nontopological soliton model (Krein et al., 1988, 1991). 

A model gluon propagator that is not singular in the infrared has also been studied (Haeri and Haeri, 
1991; Papavassiliou and Cornwall, 1991). Haeri and Haeri (1991) use 

a s (Q 2 ) 1 1 rfi , q v 

Q2 ~ Q2 +m 2 f Q 2 +4m 2 ^ 



In 

V "qcd 

based on the analysis by Cornwall (1982), where m is an effective gluon mass, and a vertex Ansatz 
constructed so as to eliminate overlapping divergences in the quark DSE. As one would expect in 
such a model, there is a critical value of the ratio m/Aq CD = R c m such that only for m/Aq CD < R° m 
does the model manifest DCSB. In this particular model, however, such a result is phenomenologically 
inconsistent because it means that DCSB requires an effective gluon mass which is significantly less 



than the value of ~ 0.5 GeV suggested by other studies [see Sec. |5.1|| . A similar internal inconsistency 
was found by Papavassiliou and Cornwall (1991) in their study of approximate coupled DSEs for the 
3-gluon and quark-gluon vertices and subsequent analysis of the DSE for the quark propagator. The 
conclusion in both of these studies is that confinement is not incorporated in the model specified by 
Eq. ( |6.39| ) whose infrared behaviour therefore requires modification if a reasonable phenomenological 
model of DCSB is to be obtained. 

There have also been studies of the non-covariant, Coulomb-gauge, quark DSE (Finger and Mandula, 
1982; Govaerts et al, 1983, 1984; Alkofer and Amundsen 1988; Langfeld et al., 1989). As in the covariant 
studies, it is found that one has DCSB if the coupling at zero momentum transfer exceeds a critical value. 
We remark that Alkofer and Amundsen (1988) used a model interaction with ct(|(5| 2 ) ~ \Q\~ 2 which was 
regularised using a cutoff mass: \Q\ 2 — > [\Q\ 2 + (J?]. In all cases a not unreasonable phenomenological 
description of DCSB was obtained. 



Summary. The body of work exemplified and summarised in this section clearly demonstrates that 
a very good phenomenological description and understanding of DCSB in QCD can be obtained using 
the DSE for the quark self energy. Given a gluon propagator whose behaviour at large spacelike-g 2 ; 
i.e., in the ultraviolet, is given by the renormalisation group in QCD and whose behaviour at small 
spacelike-g 2 ; i.e., in the infrared, yields an effective coupling strength that is > 1, then one will have 
(qq) 7^ 0. With just one or two parameters to describe the transition from the known behaviour of the 
gluon propagator in the ultraviolet to the unknown behaviour in the infrared [see Sec. |5.1| 



one can 



obtain a good phenomenology of DCSB; i.e., one can obtain values for quantities such as the quark 
condensates and the pion decay constant that agree with those extracted from experiment. This will 
be discussed in more detail in Sec. [7]. 
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6.2 Quark Confinement 



General Observations and Remarks. It is now generally accepted that QCD is the appropriate theory 
for describing the strong interaction. The lack of observation of free quarks and gluons [and ghosts] 
leads then to the assumption of colour confinement; i.e., that only quantities which transform as colour- 
singlets can be physical observables. Theoretical studies also lend credence to this expectation, where, 
for example, lattice gauge theory studies are consistent with a string tension or linearly rising potential 
between coloured objects. The fact that QCD is a gauge theory also implies that physical observables 
must be gauge invariant. The BRS invariance properties of QCD can be used to show that the S- 
matrix is gauge invariant, [see Itzykson and Zuber (1980, pp. 604-605)]. At the observational level, the 
confinement of colour is the statement that the cross-section for the production of any combination of 
asymptotic coloured states from a colour- singlet initial state must be zero. A more extensive treatment 
of some of the following discussion can be found in Roberts et al. (1992). 

The failure of the cluster decomposition property [CDP] of vacuum expectation values of the coloured 
fields in QCD [the QCD Wightman functions] is certainly sufficient to ensure confinement [Marciano 
and Pagels (1978), Greenberg (1978), Strocchi (1976,1978); also see Streater and Wightman (1980) for 
a discussion of the CDP]. The CDP can be associated with charge screening and so, in simple terms, the 
failure of the CDP means that, irrespective of how large one makes the spacelike separation between two 
coloured objects, the interaction between them never becomes negligible. This is the common intuitive 
basis for understanding confinement that underlies the idea of a string tension in QCD and provides a 
basis for the linear and harmonic-oscillator potential models. In QED the locality property of the field 
theory is sufficient to ensure that the CDP is satisfied (Strocchi, 1976, 1978), however, the proof used 
in QED to demonstrate this property fails in QCD. Crucial to the failure of this proof is the fact that 
an indefinite metric is unavoidable in local gauge theories and further that QCD is asymptotically free. 
Of course, the failure to prove the CDP is not equivalent to proving its failure. It is interesting that in 
QED2, which has a linear potential, the failure of the CDP can be explicitly demonstrated (Strocchi, 
1976). 

The n-point Green's functions of a quantum field theory are the time-ordered vacuum expectation values 
) = (vac|T(0(xi) . . . 0(x n ))|vac), where here the field operators generically represent 
gluon, quark, and ghost field operators in QCD. The Wightman functions are simply Green's functions 
without time-ordering; i.e., W^ 1 ' = (vac|0(xi) . . . 0(x n )|vac). The CDP is that as A — ► +00 

W {n \x u . . . , x h x j+1 + Xa,...,x n + Xa) — -> W u \x u . . . , Xj ) W {n ^\x 3+U . . . , x n ) , (6.40) 

where a is an arbitrary spacelike vector. The actual CDP requirement is more than this. In theories 
with a mass gap [no zero mass states] the difference between the LHS and RHS of Eq. ( |6.4U ) must 



decrease exponentially, whereas for theories without a mass gap [e.g., QED 4 ] the decrease must be at 
least as fast as 1/A 2 , (Strocchi, 1976, 1978). No explicit proof that the CDP fails for coloured Wightman 
functions in QCD is currently known, although the fact that the usual QED proof of this property does 
not apply in QCD is certainly suggestive. Clearly, if all the W in Eq. (|6.40|) are colour singlets then 
the CDP is expected to hold. The failure of the CDP property does not necessarily entail the absence 
of free coloured states in the spectrum of the theory, e.g., in Greenberg (1978) an effort was made to 
construct an explicit confining, nonrelativistic, harmonic-oscillator model in which the CDP fails. 

A different perspective is discussed in Nishijima (1986) where the non-Abelian nature of QCD and the 
necessity of an indefinite metric give rise to confinement as a result of the complete cancellation of colour 
states from the S-Matrix, in a manner analogous to that in which the longitudinal and scalar photons 
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decouple in QED. The mechanism is associated with the extended gauge transformations manifest in 
the BRS invariance of QCD. In simple terms, this picture of confinement corresponds to the idea that 
after averaging over gauge field configurations only the colour singlet pieces of operator expectation 
values remain. This can also be related to the breakdown of the CDP for coloured states. 

Another sufficient condition for confinement is if QCD has only colour-singlet free [i.e., on-shell] states. 
Hence confinement might be equivalent to the statement that the propagator of a coloured state should 
have no singularities on the real, positive p 2 -axis, [see, e.g., Cornwall (1980), Krein et al. (1988,1990), 
Gogokhia and Magradze (1989)]. Note that any singularities on the negative [i.e., spacelike] p 2 -axis 
would imply tachyonic behaviour and hence are not allowed by causality considerations. To discuss this 
possible form of confinement we will concentrate on the quark propagator, the most general form of 
which [in an arbitrary covariant gauge] is given in Eq. ( |2.68|) . This statement means that the quarks 
have no mass-pole M 2 (p 2 ) ^ p 2 for any p 2 > 0, and so cannot go on mass shell. It might be argued 
that if the propagator has a timelike singularity, then it should have a mass-shell and hence there 
appears to be nothing to prevent us from defining the corresponding free state. [It is always possible to 
imagine M(p 2 ) having some very unusual behaviour at the singular point of the propagator, which would 
complicate this argument.] This suggests that confinement might be synonymous with the requirement 
that coloured particles and bound states have no poles in their propagators. [As a counter-example 
one might consider QCD 2 in which, at leading order in a 1/N C expansion, confinement is ensured in a 
conspiratorial fashion: the quark propagator has a mass pole but the vertex in the BSE has zeros at 
the momenta which correspond to one of the quarks going on-shell and hence there is no scattering out 
of the bound state (Einhorn, 1976).] 

We remark here that we are discussing the gauge-fixed quark propagator which is a well defined and 
nontrivial element of field theory and is an important input to any covariant calculation of scattering 
or bound state amplitudes. This propagator is gauge-dependent but does contain gauge invariant 
information; e.g., the quark condensate. Of course, a nonzero value of (qq) does not imply confinement. 
Our remarks here are meant simply to emphasise that there is gauge invariant content in propagators. 
A separation between the scales for confinement and dynamical chiral symmetry breaking is, in general, 
to be expected and DCSB is certainly possible without confinement. 

It might be argued that the analytic structure of the gauge-fixed quark propagator would be a genuine 
signal of confinement in QCD only if the presence or absence of a pole [or branch cut] was also a gauge 
invariant property. In this connection the Nielsen Identities (Nielsen, 1975) of QCD are important. 
These identities, derived using the BRS invariance of QCD, ensure that, for example, the quark mass 
value obtained from the QCD effective action is gauge invariant. This gauge invariance is a result of the 
cancellation between the gauge dependence of the effective action and that of the vacuum configuration. 
Hence, one can infer from this that the absence of a mass pole [or branch point] is also a gauge invariant 
property of the theory. Practically, however, there is a problem with the effective action approach: a 
non-perturbative calculation of the QCD effective action is prohibitively difficult. 

Given the possible implications for confinement it is of interest to examine the consequences of assuming 
that the quark propagator has no singularities on the real p 2 -axis. We see that the quark propagator 
will have a mass pole when M 2 (p 2 ) = p 2 provided that Z(p 2 ) is not zero in the neighborhood of this 
value of p 2 . Alternatively, we might have that Z goes to zero in just the right way to cancel a singularity 
of the type l/[p 2 — M 2 (p 2 )] or that Z should be zero over a finite region containing just this value of p 2 . 
In the absence of compensating zeros in Z then we would be left with the requirement that M 2 (p 2 ) ^ p 2 
for any real p 2 . On perturbative grounds we expect the same behaviour in the asymptotic timelike and 
spacelike regimes. The only way that M(p 2 ) can be finite in the deep timelike and spacelike regions of 
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p 2 and yet never have a pole at M 2 (p 2 ) = p 2 , for any real p 2 , is if M(p 2 ) is discontinuous across the p 2 
line or if Z(p 2 ) has a compensating zero. The possible implications of complex poles in the propagators 
are discussed in Sec. f| 

As a further observation regarding the absence of mass poles for coloured states we remark that this 
would imply, for example, that the quark and gluon propagators have no singularities on the real axis. 
Consider then Eq. ( |6.9| ) and assume that T u has no poles [e.g., Y v — > 7^ as for the rainbow or ladder 
approximation] and that the integral exists. We see that the RHS of this equation will be pure imaginary 
if X is initially pure real, and hence in general this DSE will give rise to a complex self-energy. This 
represents an additional complication if this picture of confinement is pursued. We note, however, that 
if we assume that it is reasonable to begin with the Euclidean transcription of Eq. ( |6.9j ) then the quark 
self-energy will be real. Apparently then, in this confinement scenario the naive Wick rotation is not 
valid, which should not be surprising. 

There is one further point that it is interesting to consider here. Marciano and Pagels (1978) have 
discussed the possibility of confinement occurring through the power counting of divergences, which 
were presumed to arise in a full solution for the various propagators and proper vertices of QCD. If 
confinement results from the non-Abelian nature of QCD then these divergences would also be expected 
to have non-Abelian origins. At the simplest level one could imagine, for example, associating with each 
coloured leg of an n-point Green's function a factor y/e, while for colour- singlet bound-state legs the 
factor would be unity. Similarly, for proper [i.e., one-particle irreducible] vertices coloured legs would 
carry a factor l/y/e and colour- singlet legs again would have the factor 1. It follows that S and D ~ e, 

T v ~ 1/ e^ 3 / 2 ) for the quark-gluon vertex, T§ ~ 1/e for the quark-pion vertex, etc. It is straightforward 
to see that any process with only colour- singlet external legs is of 0(1), while any process with m 
coloured external legs will be of 0{e <jri / 2 ' > ) and so will be completely suppressed as e — ► 0. In any 
perturbative QCD calculations the selection of only those soft final state interactions which give rise 
to colour singlet asymptotic states is thus ensured. [Many other power counting schemes are possible 
and may incorporate ghost fields etc. but this illustration exemplifies the general idea behind them.] 
This is an example of the failure of the CDP for non-colour-singlet asymptotic states, whereas it would 
not fail for allowed colour-singlet processes. It can simultaneously be considered as an example of the 
absence of coloured states since free quarks cannot propagate as e — > 0. We do not argue that this is a 
realistic picture of confinement, rather we take it as a lesson that confinement might arise from subtly 
cancelling and possibly momentum-dependent divergences in the Green's functions and proper vertices 
of QCD. 

To summarise : The failure of the cluster decomposition property is sufficient to ensure confinement and 
forms the basis for many intuitive arguments and simple nonrelativistic models of confinement. One 
can think of the CDP as a measure of charge screening: having the CDP in a theory means that the 
interaction is screened and falls off with separation. Failure of Eq. (16.40 ), when W^' and W^ n ~^ are 



not colour singlets, is consistent with the notion that colour-charge is anti-screened. There is as yet, 
however, no proof that the CDP fails in QCD. In simple models the failure of the CDP appears linked 
to this infrared behaviour. The absence of a pole at timelike momenta in a particle's propagator is also 
sufficient to ensure confinement [although not necessary, as the QCD2 example shows (Einhorn, 1976)]. 
It is certainly conceivable that the failure of the CDP and the absence of a pole in the propagator of 
confined particles are related and that one is not possible without the other. This interesting speculation 
receives some support in the simple model described in the next section. Explicit demonstrations of 
confinement in more realistic models are still wanting, however. 
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Examples of Confining Models. Let us consider two simple models of the quark DSE which indeed 
lead to a confined quark and gluon. More details of this work can be found in Burden et al. (1992b). 
We work in Euclidean space and use the notation of Burden et al. (1992b), where we write for the quark 
self-energy E(/c) = S~ 1 (k) — Sq (k) with the renormalised quark propagator 

S(k) = 1 - Z ^ (6 41) 

y ) 27 • kA(k 2 ) + B{k 2 ) i^-k + Mik 2 ) V ' } 

= -i<y ■ ka v (k 2 ) + a s (k 2 ) (6.42) 
and the perturbative quark propagator So (A;) = 1/(27 ' P + m )- 

The first model of interest is obtained with the choice (Munczek and Nemirovsky, 1983) 

g 2 D^(k) = D%{k) Stt^D k 2 5\k) (6.43) 
where D° v (k) = (5^ - k^/k^/k 2 and 

T,(k,p) = ^. (6.44) 

This is an infrared dominant model which underestimates the strength of the interaction in QCD at 
\aige-k 2 . In a study whose focus is confinement this is not unreasonable since confinement is expected to 
be a predominantly infrared effect. The quark-quark interaction resulting from Eq. ( |6.43| ) is a constant 
in configuration-space and so the usual proof of the CDP will fail. Using Eq. ( |6.43| ) and ( |6.44| ), we find 
that Eq. ( |6.7|) reduces to the following pair of coupled, algebraic equations for Z(p 2 ) and M(p 2 ) 0: 

h = 2M 4 -3mM 3 + M 2 (2p 2 -D + m 2 )-3mp 2 M + p 2 m 2 . (6.45) 



2M-m 

For m = there is a dynamical chiral symmetry breaking solution of these equations (s = p 2 ) 




M(s) = { V 2 D s ' s < \D ( 6 .4 6 ) 

s<\D 

Z{s) = { s ( ^ , 4D_ 1 > \ ^ s>1 _ D (6.47) 

The CJT effective action [Cornwall et al. (1974)] evaluated at its extremals is (Stam, 1985; Sec. |7.2j ) 
when evaluated in Euclidean space: 

V[S] = li0f {^HS,\q)S{q)] + fa\Z(k)S(k)]} (6.48) 

and is the same as the auxiliary field effective action in this case (McKay and Munczek, 1989). To 
determine whether chiral symmetry is dynamically broken for m=0 one must evaluate the difference 



(3 = V[Z = (s/2D) ( yjl + AD/s - 1 ) , M = 0] - V[Z = Eq. ( $Af} ), M = Eq. (|6\46|) 1 
D 2 ( . I I 



NcNf 32^{ Aln2 ~T> >0 (6 ' 49) 



"In deriving these equations we have made use of the fact that 



d\f{{ q + p) 2 ) [ ^5%) = \ftf) 
q z p z 4 



We remark that Eq. (3.43) represents a Landau gauge propagator in contrast to the work of Munczek and Nemirovsky 
(1983) where a Feynman gauge propagator was used. To reproduce these earlier results simply write D = rf /2. 
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which reveals that the ground state with a nonzero fermion mass function is dynamically favoured and 
hence chiral symmetry is dynamically broken in this approximation. The dynamical chiral symmetry 
breaking solution is also a confining solution since S(p) constructed from these functions has no pole 
on the real s axis. This feature survives for m > with M and Z smooth on the real s axis: in the 
spacelike UV limitB 

M(s) =m(l + —) and Z(s) = 1 - — ; (6.50) 



s J 



while in the timelike UV limit 



M(s) = M (s) + m^j^ and Z{s) = Z (s) + m~ , (6.51) 



where M and Z are the m = solutions of Eqs. ( |6.46| ) and (|6.47 ). In this case the integrable 5 A (k) 



singularity in the model gluon propagator of Eq. Q6.43) lies behind the absence of singularities in the 



quark propagator and has been used as the starting point for confining, phenomenological models of 
nonperturbative QCD [see, e.g., Praschifka et al. (1989) and Frank et al. (1991)]. 



The second model has the same form for the gluon propagator as Eq. ( |6.43| ), but uses the Curtis- 



Pennington vertex of Eq. fl6.5|) . This represents a test of whether the absence of a pole in the quark 
propagator is sensitive to the details of the quark-gluon vertex. Substituting Eqs. (|6.5|) and ( |6.43| ) 
into fl6.71 ) then, instead of Eqs. (|6.45|) , we obtain the following pair of coupled, nonlinear, differential 
equations: 

^(x) = ^+2 X + f (x) (C(x)m- M (x)), (6.52) 
X xQ(x) 

C(x) = ^ X) ^f- I ax) + 2(1-C(x)) (6.53) 
x + n A [x) 

where x = s/(2D), m = m/\/2D, M(s) = \ / 2Dfi(x) and Z(s) = ((x). [Note that the transverse addi- 
tion Eq. (|6.3|) does not survive when Eq. ( |6.43|) is used and so the Ball-Chiu vertex of Eq. ( |6.2j ) gives 
the same results.] The boundary conditions for these equations are, of course: 

fi(x — > +oo) = m and ((x — * +oo) = 1 (6.54) 

consistent with Eq. ( |6.50| ). The pragmatic benefit of Eq. ( |6.43j ) is again obvious: we are confronted 
with differential equations instead of integral equations [e.g., Roberts and McKellar (1990), Haeri and 
Haeri (1991), Williams et al. (1991), Hawes and Williams (1991), Hadicke (1991)]. Hence we can 
obtain solutions for both spacelike and timelike momenta, which enables a study of the singularity 
structure of the quark propagator. [The infrared dominant k~ 4 model also has this feature in axial gauge 
as discussed in Ball and Zachariasen (1981).] These differential equations take a particularly simple 
form when expressed in terms of <5V(x) = ( (x) / (x + fi 2 (x)) and &s( x ) =o r y(^)/x(x) (x = k 2 /(2D), 
a v (k 2 ) = a v (x)/(2D), a s (k 2 ) = a s (x)/V2D): 

2 

o s \x) = 2[mav(x) —as{x)) , av'{x) = — — [&v{%) (x + 1) + mW s (x) — 1] . (6.55) 

We begin by considering the m = case [in which case Eqs. ( |6.55|) decouple]. Obviously there is the 
trivial \x = solution that corresponds to a Wigner-Weyl realisation of chiral symmetry. In this case 

1 e~ 2x 

C(x) = l- — + C (6.56) 

2x x 



^Deriving these results is easiest if one uses the equations: 

(1- Z) = (DZ 2 )/(s + M 2 ) and (M — mZ) = (2DMZ 2 )/(s + M 2 ) 
which are equivalent to Eqs. fl6.45|). 
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with C an arbitrary constant. The choice C —1/2 yields a propagator with no singularity on the real 
x axis and is required if this solution of the differential equation is to be consistent with the integral 
equation, so that there is no singularity within the Euclidean integration domain [in agreement with 
footnote 3 in Munczek (1986)]. The solution can therefore be written as: 

o r — 1 4- p~ 2x 

Ms) = 2 J > °s(x) = . (6.57) 



There is also a nontrivial solution. In searching for this solution one finds that for m = Eqs. ( |6.52j ) 
and ( |6.53[ ) are scale invariant; i.e., given \i and ( as solutions then 



x -\- A 2 u (2?) 

Hx{x) = V(x) an d Ca(z) = : — tt~- C0*0 (6-58) 

X + jJ, z {X) 

are also solutions. It is therefore sufficient to find the solution for which /^(0) = 1. These scale trans- 
formations correspond to: 

(?s(x) = Xas(x) and a v (x) = &v(%) (6.59) 



and the scale invariance of Eqs. ( |6.52j) and ( |6.53D is a result of the decoupling of Eqs ( |6.55| ). It is, of 



course, easiest to find the solution from Eqs ( |6.55|) which yield ay of Eq. (|6.57|) and 

a s (x) = e~ 2x (6.60) 

which correspond to: 

2x 2 H ,_ 4x 3 + [l + (2x-l)e 2 f 

KX} l + (2x-l)e 2 -' UXj 2xe 2 *[l + (2x-l)e 2 -] 1 ' 

To determine whether chiral symmetry is dynamically broken in this model one must evaluate the 
difference 

V[a v = Eq. W S = 0] - V[a v = Eq. O, a s = Eq. CT) ] 



and hence the chirally symmetric solution is dynamically favoured in this pointed out by 

Munczek (1986). The underlying reason for this is the decoupling of the equations for ay and as- It 
is reasonable to expect that with a model gluon propagator that incorporates asymptotic freedom [i.e., 
one which behaves as ~ [q 2 ln(g 2 )] -1 in the deep spacelike region] these equations will not decouple and 
hence chiral symmetry will be dynamically broken. 

Comparing these two confining models with m = we see that the choice of vertex has a quite dramatic 
qualitative and quantitative effect on the nature of the solution and the realisation of chiral symmetry. 
Qualitatively, the dressed vertex has ensured that the quark propagator is an entire function, whereas 
in the bare vertex case it had a branch point at x = 1/4, and it has also eliminated dynamical chiral 
symmetry breaking from the model. Quantitatively, there is a dramatic enhancement of the propagator 
at large (— x); i.e., the deep timelike region. This is easily seen since, using Eq. (|6.44|) , one finds that in 
the bare vertex case ay(x = k 2 /(2D)) — > 2 and as{x) — ► 1\f—x as x — > — oo which are to be compared 
with Eqs. Q637D and ( ggg ). 
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In the second model with m^O there is no trivial \p = 0] solution and solving Eqs. ( |5.52| ) and ( |6.53| ) 
[or, equivalently, Eqs. (|6.55| )1 is a little more difficult. These equations can be combined to yield: 

d 2 , s ( 1 \ d , . / 1 + m 2 \ , , 4m 

- Mx)+i ( }+ -)- Mx)+ ^ + — r y s(x) = _ (6 . 63 ) 

from which one obtains (x = y 2 ): 

Mv) = exp(-2|/ 2 ) Jx(my) + — / d^K^mOJM)^- (6-64) 

with ci an undetermined constant and J\ and Bessel and modified Bessel functions of order one, 
respectively.! [In the limit m — > Eq. ( |6.64 ) reduces to Eq. ( |6.60 ); the dynamical chiral symmetry 



breaking solution. In the following we set c\ = and concentrate on the behaviour of that part of the 
solution associated with explicit chiral symmetry breaking.] The form of ay follows from Eq. ( |6.55| ) 

(y 2 = x): 

Mv) = ^siv) + ~^ 8 (y)) ■ (6.65) 

It is clear that av(y) and ~<Js{y) are even under y — > —y and, hence, are functions of x = y 2 . Further, 
after a little thought it becomes clear that Ws{x) and oV(x) are entire functions in the complex x 
plane. A combination of numerical and analytic methods establishes that both as and <5V are positive 
definite on the real x-axis. It is of interest to determine the asymptotic behaviour of these functions 
at large spacelike and timelike x = k 2 /(2D). In the deep spacelike region [x ~ +oo] one finds easily 
that Ws{x) = m/x and Wy(x) = 1/x. In the deep timelike region (x ~ — oo) the analysis is a little more 
complicated. Using Laplace's method one finds that at large E = \J —x: 



a s (E) = exp (2[E - m} 2 ) [l + erf (V2[E - m])] (6.66) 



2E 3 

which yields the following asymptotic timelike behaviour: 



a s (E) = J^! e xp(2[E-m] 2 ) , (6.67) 



ME) = V^^n 2 ^-^!, 1 -^-!^) • (6 - 68) 



This again is a dramatic enhancement compared with the bare vertex solutions: oy ~2(1 — mj J\ — x) 
and as — 2(Jj — x — m). 



Hidden in this discussion of as and ay is the behaviour of the functions fi and ( [which are related to 
the mass function, M, and wave function renormalisation, Z: see after Eq. ( |6.53| )1 that can be expressed 
[K(x) = l + a s \x)/{2a s (x)]: 

K x ) = =4~T = t5~T and COc) = = \x + V 2 (x)] o r s {x)K(x) . (6.69) 
ay(x) K[x) m L J 

The propagator is very often discussed and modelled in terms of these functions [equivalent to M and 
Z]. Our analysis leads us to the conclusion that K(x) > for all finite, real x and hence so is /z; i.e., 
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Fig. 6.6. A plot illustrating the m-dependence of the zero of the prop- 
agator denominator in the {Z,M} representation for S(p). 

the mass function has no zeros or singularities for real x. In Fig. 6.6 we have plotted [x + f^ 2 {x)} for 
a range of values of m. The striking observation is that this has a m dependent zero. A pole in the 
propagator is only avoided, therefore, because ((x) also has a zero at this point, Eq. (|6.69|) ; i.e., ( and 
x + /i 2 have coincident zeros. [We remark that a zero in ( with [i finite entails a pole in the functions 
A(x) and B(x) of Eq. flOiP -l 

To summarise : In the first model [the bare vertex] we see that M(p 2 ) diverges in the timelike region 
[Euclidean p 2 < 0] so that a pole in the quark propagator is avoided and hence the theory is confining. 
In the second model [the CP or BC vertex] for m = there is confinement [the quark propagator is 
entire] and the phase with no DCSB [i.e., with M(p 2 ) = 0] is energetically favoured. When m ^ in the 
second model we find DCSB and that confinement occurs through cancelling zeros in the numerator and 
denominator of the quark propagator. It is interesting to speculate that the effective gluon propagator 
is primarily responsible for confinement and that the vertex modifies the particular manifestation for 
the quark propagator. However, other models need to be studied before any such conclusion can be 
drawn. 



Application of a Confinement Test. The confining models discussed above could be studied in the 
timelike region because their simple form allowed them to be written as coupled differential equations. 
In the general case an integral equation must be solved and the extension into the timelike region is 
highly nontrivial. It is still possible to probe the timelike behaviour of the quark propagator obtained by 
numerical solution of a model quark DSE. In order to determine whether the quark propagator obtained 
as a solution to our DSE can represent a confined particle we follow Hollenberg et al. (1992) and adapt 

§In this expression the additional solution of the associated homogeneous equation: 

-^-exp(~2 2/ 2 )y 1 (4m 2/ ) , 
Amy 

x = y 2 ; Y\ a Bessel functions of order one, has been discarded because it is irregular at y = 0. 
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a method commonly used in lattice QCD to estimate bound state masses. We write as before 

asip2) = jtAWY + BW (6 - 70) 



^s(T,x) = I 4^r^ iP4T+f ^s(p 2 ) ■ (6.71) 



and define 

t 
(2^) 4 

This is the scalar part of the Schwinger function of the model quark propagator. If we now define 

A S (T) = J d 3 xA s (T,x) (6.72) 

and, for notational convenience, E(T) = — In [Ag(T)], then it follows that if there is a stable asymptotic 
state with the quantum numbers of this Schwinger function then 

dE(T) 

lim A- ' = m ; 6.73 

where m > is the asymptotic [on-shell] mass of this excitation; i.e., this limit yields the dynamically 
generated quark mass evaluated at the mass pole. If the limit in Eq. ( |6.73| ) exists for a given propagator 
then we can conclude that the propagator has a timelike mass pole, which suggests that the associated 
excitation is not confined. 

For example, one can consider the model of Burden et al. (1992b) discussed above and one finds 

alE{T) 



alT 



kT (6.74) 



where k is a constant and hence the limit in Eq. Q6.73 ) does not exist. In this case the Schwinger function 



does not satisfy the cluster decomposition property and hence the quarks are confined (Roberts et al, 
1992). Alternatively one may say that through self interaction the quark acquires an infinite dynamical- 
mass. This provides another way of understanding the claim that the model of Burden et al. (1992b) 
is confining. 

Another example of this approach is to apply it to the IR vanishing gluon propagator of Eq. (|6.29|) . We 
set a s to a constant here for simplicity and then, for the boson analogue of Ag(T), one finds: 

A(T) a rVexpf-^Vcos^fl -amf^Tl . (6.75) 



by/2 

Note that the Schwinger function in this case is not positive definite, which is an easily identifiable signal 
in A(T) due to the pair of complex conjugate poles. This violates the axiom of reflection positivity. 
It follows from this that Eq. (|6.75|) describes a field with a complex mass spectrum and/or residues 
that are not positive. This is appropriate for particles that decay and forms the basis of the argument 
that such a propagator allows coloured states to exist only for a finite time [of the order of 1/6] before 
hadronising; i.e., that the propagator describes confined gluons (Stingl, 1986; Habel et al. 1990a, 1990b; 
Zwanziger, 1991). 

To illustrate this procedure for an explicit numerical solution we return to the example of the model 



quark DSE discussed in Sec. [Tljand attempt to determine whether the model gluon propagator specified 
by Eqs. ( |6.28| ) and ( |6.29| ) leads to quark confinement; i.e., the absence of free quarks in the QCD 
spectrum (Hawes et al. (1993). In applying this method here it is obvious that numerical evaluation of 
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the Fourier transforms required in using Eq. ( |6.73D will be hindered by numerical noise as T is increased. 
In order to minimise the effect of this noise on the derivative, we fitted A${T) to a form 

Cexp(-mT) (6.76) 

and extracted the derivative from this fit. [Importantly, there is no indication in the results of the 
structure suggested by Eq. ( |6.75|) .1 This was particularly useful with the propagators obtained using 
small values of b which had large dynamical masses [as one would expect since the condensate is large in 
this case] and hence a rapid decline with T. Results were found for the confinement test in the following 
cases: 1) The propagators obtained with lnr = 0.1 and b 2 in the range [0.1,1.0]; 2) The propagator 
obtained with lnr = and b 2 = 0.35 which yields the largest value of — (qq)^ on the (6 2 ,lnr) domain 
considered; 3) Two propagators obtained with (6 2 ,lnr) = (0.1,0.6) - one using the Ball-Chiu vertex 
and another using the Curtis-Pennington addition. The results obtained by fitting the form Eq. ( |6.76| ) 
to our numerical output are presented in Table. K% 



Table 6.2. Asymptotic dressed-quark-mass values for the family of prop- 
agators with lnr = 0.1, the propagator which showed maximal DCSB 
(at lnr = and b 2 = 0.35), and for lnr = 0.6, b 2 = 0.1 with both the 
Ball-Chiu and Curtis-Pennington vertices. These values were arrived at 
by exponential fit and match the limits of the asymptotic mass curves. 



lnr 


b 2 


^free 


C 


Comments 


.1 


.25 


0.410 


0.664 


B-C vertex 


.1 


.35 


0.296 


0.633 


j; 


.1 


.45 


0.176 


0.619 


!) 


.1 


.54 


0.0275 


0.619 


n 


.0 


.35 


0.406 


0.667 


n 


.6 


0.1 


0.210 


0.648 


B-C vertex 


.6 


0.1 


0.210 


0.507 


C-P vertex 



In Fig. 6.7 we present plots of E'(T) for the family of propagators obtained with lnr = 0.1 and this 
clearly illustrates that an unambiguous determination of the dressed-quark-mass is possible. It will 
be observed that the mass decreases with increasing b 2 . This is easily understood in terms of the 
chiral phase transition: as b increases beyond b c there is no DCSB and massless current quarks remain 
massless. Since the behaviour of all the other solutions we obtained was qualitatively the same as that 



described by the results presented in Table 6.2| and Fig. 6.7 we infer that the model considered here 
does not yield a confining quark propagator. We also note that the rainbow approximation studies of 
the fermion DSE in Alkofer and Bender (1993), which address the question of confinement by a direct 
continuation to Minkowski momentum space, found a quark propagator with a pole at timelike p 2 ; i.e., 
a non-confining propagator. We remark that, within numerical noise, the Curtis-Pennington addition 
made no difference to the dressed-quark-mass value extracted in the cases considered and only slightly 
reduced the normalisation constant C . Clearly, the Curtis-Pennington addition leads only to a minor 
quantitative effect in this part of our study too. 

It is worthwhile to recall the suggestion [Habel et al. (1990a,b)] that a model gluon propagator with 
the infrared form of Eq. (|6.29| ) would lead to a confining quark propagator of the form 

ij ■ P + Co 



(17 • p + ci) (17 • p + c\) 



(6.77) 
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Mass curves, ln(x)=0.1, Minimal Ball-Chiu Vertex 
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Fig. 6.7. Dressed-quark-mass curves for the family of propagators with 
the minimal Ball-Chiu vertex and lnr = 0.1. The masses are in units of 
Aqcd- 



where Co and c\ are real and complex constants, respectively; i.e., to a fermion propagator with complex 
conjugate poles just as the gluon propagator has; confinement being realised through the absence of 
poles on the real timelike axis. As Eq. ( 6.75| ) shows, such a propagator has a characteristic signature 
in Ag(T) which we do not see in Fig. 6.7. In this model then it is clear that a fermion propagator of 
the type in Eq. ( |6.77|) does not arise. This does not eliminate the possibility that it can arise in the 
approach of Habel et al. (1990a, 1990b), however, since the rational-polynomial Ansatze employed for 
the vertex functions therein leads to a completely different quark-gluon vertex to that used here. The 
above results do however suggest that a fermion propagator of the type in Eq. ( |6.77p cannot arise if the 
quark-gluon vertex is free of kinematic singularities. 



6.3 Euclidean ±± Minkowski Continuation 

We have discussed this issue in Sec. |2.3| and here we simply put these earlier remarks into context by 



way of a few examples. 

The model of Eqs. ( |6.43| ) and ( |6.44| ) yields a quark propagator which has complex conjugate branch 



points with spacelike real parts. Clearly, therefore, the naive transcription rules of Sec. |2.3| will not 
be suitable in this case. In the modification of this model, Eqs. fl6.5p and ( |6.43| ), one obtains a quark 
propagator that is an entire function with an essential singularity at timelike infinity. In this case it is 
not possible to close the integration contour and hence a "Wick rotation" is simply not possible. This 
latter case illustrates that the inability to perform a Wick rotation appears even in models which go 
beyond the rainbow approximation. 

Detailed analyses of the analyticity properties of quark propagators obtained with other model gluon 
propagators in rainbow approximation have been carried out (Stainsby and Cahill, 1992; Maris and 
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Holties, 1992; Stainsby, 1993, Sec. 3; Maris, 1993, Sec. 3.3). In the studies by Stainsby and Cahill 
(1992) and Stainsby (1993) three forms of gluon propagator were considered: 



with [s = q 2 



D,M) = V ^|r^ (6-78) 



/ x 37TY 2 / S\ ill , „ . 

= Ti^exp -- + 7 ^, (6.79) 



4A 2 V A 



11 In l + e + 



s 



T 

QCD 



where x = 1-14 GeV, A = 0.002 GeV 2 , A QCD = 0.19 GeV and e = 2.0, taken from Praschifka et 
al. (1989) where it was shown that such a gluon propagator is well suited to a description of low energy 
hadronic physics; 

®2(s) = — and otds) = C 2 exp f ) (6.80) 

s + s \ s J 

with Ci « 1, s » (0.0707) 2 GeV 2 and C 2 » 50.0 and s « (0.245) 2 GeV 2 chosen so that the solutions 
obtained for A(p 2 ) and B(p 2 ), with p 2 on the real spacelike axis, were very similar to those found with 
a\. Solving the integral equations directly, with no approximation of the kernel, it was found that in 
each case there were pairs of complex-conjugate, logarithmic branch points with spacelike real part in 
the functions A(p 2 ) and B(p 2 ) and hence in the quark propagator. The position of the branch point 
depended on the form of the propagator but not its nature or existence: 

ai : s B w 0.22 + 0.25 % , a 2 : s B » 0.015 + 0.33 % , a 3 : s B « 0.06 + 0.41 i . (6.81) 



Maris (1993) considered two model forms of gluon propagator 



, ^\ 47ra(g 2 ) 

^(g) = 1 V - I — ^ — (6.82) 



with 



q 1 (g 2 )= / A/ " N and Wg 2 ) ^ ^ 4^) 
lnfr + ^) ^ ^ 2 

\ QCD / 

The first form is a simple extension of the one loop renormalisation group result for the running coupling 
in QCD to small-g 2 and the second form involves an additional "confining term", S A (q). The DSEs 
obtained were analysed using the approximation of Eq. (|4.7|) to obtain differential equations. 



In the case of a,\ there are five singularities, two of which are simply those of the running coupling; 
i.e., at x = —t and x = —r + 1, where x = p 2 /Aq CD , and can therefore be discarded as artifacts of 
the model. Of the three not trivial singularities, two are complex conjugate branch points and one is 
a branch point on the timelike p 2 axis with — r + 1 < x B < 0. This latter singularity might be used to 
define the quark mass shell. 

The singularity structure of the a 2 model is interesting and complicated. The most interesting feature 
is that the singularity on the timelike axis bifurcates into a pair of complex conjugate branch points in 
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a) 





b) 



<- 




iS 

o 
o 

iS 



Fig. 7.1. Shown are two equivalent forms of the BSE; i.e., (a) the proper 
meson-quark vertex [r me son] and (b) the Bethe-Salpeter amplitude for the 
meson [i\ = ST mcson S]. 



the presence of the "confining" term thereby eliminating the quark mass shell. In addition, the branch 
points that were present with a\ alone also bifurcate and move as rf is changed. 



In each of these examples the naive Wick rotation, given by the transcription rules of Eqs. ( |2.83|) , ( P-84Q 
and ( |2.85 ), does not yield the correct continuation of the Euclidean space DSE to Minkowski space. 



This serves to emphasise the points we made in the last subsection of Sec. 2.3 



7 Applications to Hadron Structure 

In this section we will apply the knowledge harvested from the studies of the DSEs in QCD, all of which 
is necessary, to the study of the kinematic and dynamical properties of hadrons. We will work for the 
most part in Euclidean space. 



7.1 Coupled Dyson-Schwinger and Bethe-Salpeter equation Phenomenology 

One direct application of the tower of DSEs to hadron phenomenology is to solve the bound state 
equations using a kernel constructed from the dressed quark and gluon propagators discussed above. In 
Fig. 7.1 we illustrate the BSE for mesons. We are assuming no n- flavour-singlet mesons and so only the 
homogeneous BSE is shown. In this figure we employ the convention that i\ = ST meBoa S and the BSE 
quark-anti-quark kernel is K = (iV)(iD)(iT) + ■ ■ -. Retaining only the first term shown for K yields 
the dressed-ladder approximation for the BSE. 

There have been a number of studies of the covariant BSE for mesons using a kernel constructed from 
a model dressed gluon propagator, whose ultraviolet behaviour is fixed by the QCD renormalisation 
group and whose infrared behaviour is based on the studies described in Sec. |5.1| , and a dressed quark 
propagator obtained as a solution of the rainbow approximation DSE with the same model gluon 
propagator. In a series of studies (Cahill et ai, 1987; Praschifka, 1988; Praschifka et ai, 1989; Stainsby, 
1993, Sec. 4; Stainsby and Cahill, 1993) it was shown that this system of equations provides a very 
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good description of the u- and d-quark sector of the mesonic spectrum with the choice 

( , \ 



dm = s 



q 2 In ( t + 



(7.1) 



V y \ ^QCD / / 

with C = (37r) 3 and r = 3. These model studies incorporated quark confinement in the sense discussed 
in Sec. B73 and will be discussed further in Sec. fOI. 



In addition, there have been phenomenologically successful studies of this coupled system using an 
infrared-finite form of the model gluon propagator subject to the approximation of Eq. ( |4.7|) (Dai et 
ai, 1991; Aoki et ai, 1991) and noncovariant, Coulomb-gauge studies of a pion BSE whose kernel is 
constructed using a dressed quark propagator obtained as a solution of the Coulomb-gauge quark-DSE 
with a bare gluon propagator (Govaerts et ai, 1984) and a model dressed gluon propagator which is 
singular in the infrared (Alkofer and Amundsen, 1988). 

We cannot summarise all of these studies and choose instead to discuss, as a particular example, 
recent studies by Munczek and Jain (1992) and Jain and Munczek (1993) which consider pseudoscalar- 
and vector-g-g bound states of u-, d-, s-, c- and 6-quarks. The treatment necessarily used significant 
approximations and was carried out in Euclidean space, which required extrapolation into the timelike 
region in order to extract the meson masses. The second of these articles is somewhat simpler and more 
systematic and we summarise its arguments and results here. It is important to note that solving the 
BSE with a kernel that is consistent with the solution of the quark DSE is necessary to guarantee the 
Goldstone-boson nature of the pion in the chiral limit (Delbourgo and Scadron, 1979). 

As in the studies summarised briefly above, the principal initial approximation in this model is that the 
quark-gluon vertex is simply given by = 7 M and the q — q scattering kernel is given by an effective 
single-gluon propagator. This is a ladder- type approximation with a modified gluon propagator D^. 
In this approximation, the BSE for a qq bound state of mass m 2 g = p 2 can be written in Minkowski 
space as, 

1fi-7r r d 4 k 

(Sf)- 1 (q + ap) X (p,q)(Sn- 1 (q-(3p) = -%U / -^-^xip^D^ {k - q) , (7.2) 



3 J (2tt) 4 

with a + (3 — 1. The quark propagators are obtained by solving the DSE for a fermion flavor /, in the 
rainbow approximation, 

(Sfy^q) = 4- mf -i™EJ J^ lllS f{k) lv D^{k-q) . (7.3) 

Note that for convenience a s (— (k — q) 2 ) has been absorbed into the definition of the effective gluon 
propagator here [as in Jain and Munczek (1993)]. The Euclidean space form of these equations is [see 



Sees. 2.3 and FT3 



(S f )- 1 (q + ap) X (p,q)(ST 1 (q-m = ~~ J J^T A 1 ,xM lv D , v {k - q) , (7.4) 

1fi7r r d 4 k 

(S-f)-\q) = vy ■ q + m, + J j^l^ {k) lv D, v {k - q) . (7.5) 

The ladder-rainbow approximation has been argued to be most appropriate in Landau gauge (Fomin 
et ai, 1983; Atkinson et ai, 1988c) in which case the model dressed gluon propagator can be written 

^(A) = - ^W) . (7.6) 



N8 



Jain and Munczek (1993) choose the momentum dependence of the gluon propagator taken to be 



D(k 2 



k 2 ln(r + x) 



ln[ln(r + x)} 
ln(r + x) 



+ D IR (k 2 



(7.7) 



where x = k 2 /A 2 QCB and A QC d ^ 200MeV is theQCD scale parameter. In Eq. (|777D, d M = 12/[33 - 2JV>], 
6 = 2fo/0i, Pi = N f /3 - 11/2, /3 2 = 19iV//12 - 51/4 and N f is the number of quark flavors, [taken as 
five in these calculations] and so, at large k 2 , the expression in parentheses agrees with the two- loop 
renormalisation group result for the running coupling in QCD, Eq. (|5.9|) . As in many of the studies 
described above, the parameter r regulates the infrared behaviour of the running coupling. As long 
as its value is larger than ~ 2, the results of the calculations are largely insensitive to it. The value 
adopted by Jain and Munczek is r = 10. 

The infrared effects in the model are dominated by Dm(k 2 ), which for simplicity was taken to have a 
modified Gaussian form 

D IR (k 2 ) = vra k 2 e-^ 2 , (7.8) 

with a = (0.387 GeV)~ 4 and \i = (0.510 GeV)~ 2 . Defining the classical potential by analogue with 
Eq. ( |3.4|) , including the negative sign associated with attraction and multiplying by the factor 167r/3, 
this form corresponds to 

™ = i^<' a - < w «»>($) < 7 - 9 > 

which is approximately linear out to 0.4 fm. We note that Eq. ( J7.8[ ) does not incorporate the infrared 
singularity ~ 5 4 (k) that was associated with quark confinement in Sec. |6.2| . This term was included in 
the earlier study by Munczek and Jain (1992). 



The most general decomposition for the vector, pseudoscalar and scalar bound state wave functions can 
be written as 



Xvfao) = h-£Xvo-h-P£-<lXvi-h-<l£-<lXv2+ e -<lXv3 

+ h ■ e,7 -P}X V 4 + h ■ e >7 • q]x V5 ~b-P,l- <?]e- QX V6 + ^7s7 -*X V 

with = e^apqvPaep, e 2 = -1 and e ■ p = 0; 



(7.10) 



x P (p,q) = 75 



Xj 



h-pXpi -ii-qx P2 - h-pn-q] Xi 



(7.11) 



x s (p» q) = Xso - *7 • pXsi - l i ■ iXs2 - h ■ p, i ■ <f\ x S3 ■ ( 7 - 12 ) 

In these equations, flavor indices are implicit on x Vl , X Pi an d Xsv which are functions of p 2 , q 2 and p- q. 

An important result in this study is that the dominant contributions to the vector and pseudoscalar 
wave functions is provided by X^}- The BSEs for the dominant wave functions, Xjoi are: 



x vo (p ,q ,p-q) 



with 



'vo 



A a A b D 
d A k 



{(q 2 - a(3p 2 + m a m b ) + q ■ p(a - p)}I V o + &Xx 



x P0 (p ,q ,p-q) 



A a A b D 



{(q 2 - a(3p 2 + m a m b ) + (a - (3)q ■ p}I P0 + 



(7.13) 

(7.14) 
(7.15) 
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with 

Ipo = /(04^-^o; (7-16) 

3 

A a ~A b D 

with 



X S0 (P 2 ^ 2 ,P- q) = -j— 7-^{(q 2 ~ app 2 - m a m b ) + (a - p)q ■ p}I S0 + ^X so > ( 7 - 17 ) 



= / (04^(*-9)x OT ; (7-18) 



and where 



/i /i c , ^ /i At o\ B a (q + ap) B b (q-pp) 

A a = A a (q + ap), A b = A b (q - (3p) , m a = — -, m 6 = — — 7.19 

A a (g + ap) A(g - /5p) 

and D = [(g + «p) 2 + m 2 ][(g - /V) + m\\. 



In these equations the Ax J0 contain contributions from several of the wave functions Xn- The technique 
adopted by Jain and Munczek (1993) to solve these equations was to write p • q = pq cos 9 and expand 
the functions Xji m terms of Tschebyshev polynomials, T^ n \cos9): 

Xji (q 2 ,M 2 B , cos 9^j = £ X ^(q\Ml) T< n >(cos0) , (7.20) 

where the subscript J ranges over V, P and 5 and % — 1,2, . . . over the independent wave functions. 
The equation was then projected with the zeroth order Tschebyshev polynomial and the right-hand-side 
approximated by keeping only the leading order term in the Tschebyshev expansion for % > 1, and 
by retaining the three lowest terms, an d X^, for x J0 - 

The vector and pseudoscalar meson equations were reduced to single equations for Xjq by neglecting 
all contributions on the right-hand-side of the equations for x^, xf] anc ^ Xj]>i other than the one 
involving Xjoi thus breaking the self-consistency of the tower of equations that couple the Tschebyshev 
moments. In the scalar meson case, however, the contribution due to xf^ was f° un d to be important. 
Consequently the procedure used in the vector and pseudoscalar case was adapted in order to obtain a 
pair of coupled equations for xfg an d xf^ which were solved self consistently. 

In arriving at the final equations to be studied Jain and Munczek (1993) also employed the approxima- 
tion of expanding the mass functions, m a and m^, and the functions A and B in a Taylor series around 
q 2 + a 2 p 2 or q 2 + j3 2 p 2 , depending on their argument, retaining up to first derivative terms whose contri- 
bution to the calculated meson masses was less than 10-15% . Where necessary the solutions for A and 
B, obtained on the spacelike real axis from the quark DSE, were extrapolated to timelike values of their 
arguments. Electromagnetic effects were incorporated by including the one photon exchange kernel in 
the BSE and, in order ensure that Goldstone's theorem remained valid, in the quark DSE. Isospin mass 
splittings were included by using different masses for the u- and d- quarks. The projected BSEs were 
then solved numerically for bound-state mass-eigenvalues and wave functions. The total mass difference 
between states of the type qu and qd, where q can be any one of the quarks, was obtained as a sum of 
the isospin and electromagnetic mass splittings except in the case of the pion for which it was assumed 
that the pure isospin contribution to the mass difference between 7r° and ir + is negligible. 

There are a number of important conclusions: 1) The contribution to a given meson mass due to XJi>i 
is negligible in all cases as was the contribution of all Tschebyshev moments of order greater than zero. 
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Table 7.1. Calculated pseudoscalar and vector meson masses and pseu- 
doscalar decay constants obtained by Jain and Munczek (1993) using 
renormalisation group invariant quark masses of m u = 8.73 MeV = rhd, 
rh s = 203 MeV; "constituent-quark" masses of M c = 1.54 GeV and 
Mb = 4.74 GeV, defined as the value of the mass function at q 2 = 0, and 
the values of a and /i given after Eq. ( |7.8| ). These tables were adapted 
from Jain and Munczek (1993). 



meson 


71 


K+ 


ss 


Vc 


D s 


D 




B c 


B s 


B+ 


mass (MeV) 


135 


494 




2979 


1969 


1865 








5279 


(exp.) 






















mass (MeV) 


135 


494 


703 


2821 


1872 


1756 


9322 


6126 


5249 


5149 


/m (MeV) 


93 


114 




213 


148 


125 


287 


207 


119 


119 


a 


0.5 


0.35 


0.5 


0.5 


0.25 


0.2 


0.5 


0.27 


0.10 


.09 


meson 


P 


K *+ 




J/V 


D* 

^ s 


n* 
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B* 


B* s 


B* + 


mass (MeV) 


770 


892 


1019 


3097 


2110 


2007 


9460 






5325 


(exp.) 






















mass (MeV) 


760 


976 


1184 


3100 


2187 


1997 


9460 


6277 


5415 


5290 


a 


0.5 


0.65 


0.5 


0.5 


0.3 


0.3 


0.5 


0.2 


.1 


0.07 



The net contribution being less than 10%; 2) There was some dependence of the meson mass on a, the 
relative-momentum distribution parameter, but in all cases it was less than 15%. In quoting results the 
authors settled on a value of a for which the response of the mass to changes in a was minimised; 3) 
The hyperfine splitting was in general too large but the splitting between the ground state and radial 
excitations was in good agreement with experiment. Some of these observations are apparent from 
Table [TT[ 

As remarked by Jain and Munczek (1993), a number of improvements of this study are possible. The fact 
that the splitting between vector and pseudoscalar mesons is larger than experimentally observed may 
indicate the need for an effective scalar interaction in the kernel of the BSE. This must be introduced, 
however, in such a manner as to ensure the preservation of Goldstone's theorem and hence will entail a 
modification of the quark DSE. A natural way to proceed in this direction is to go beyond ladder-rainbow 
approximation, which requires a dressed vertex in the quark DSE. 

Simple improvements in methodology are also possible, for example, one might avoid the Tschebyshev 
expansion, Taylor expansion and extrapolation procedures. In a study of the pion and scalar-diquark 
BSEs using a quark propagator based on the DSE studies of Burden et al. (1992b), which are summarised 
in Sec. |6T2] , and subject to the approximations 

X P ~l5X P0 and Xo+ « Xo+o , (7.21) 

Stainsby and Cahill (1993) solved the integral equations directly by matrix inversion on a discrete 
grid. The fact that the arguments of the quark propagators in the kernel of the BSEs became complex 
in Euclidean space, since P 2 = — M 2 at the solution, did not present a problem because the quark 
propagators were regular on the complete integration domain as a consequence of the quark confinement 
mechanism. 

This approach can also be applied to baryons using the relativistic Faddeev equations (Cahill et al, 
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1989; Burden et ai, 1989; Cahill, 1992). Such studies are facilitated by, but not founded upon, the 
observation that the scalar-diquark, (qq)o+, acts effectively as a resonance in the two-body subchannel 
in which case the Faddeev equation can be written in the form 

*(f,p) = ^/^r 0+ (p + ig + ^p)r 0+ (g + ip+^P) (7.22) 

xS([2a -l]P-p-q) S([l - a]P + q) D 0+ (aP - q) V(q; P) 

where / is the normalisation, T + is the scalar-diquark Bethe-Salpeter amplitude, which, at the same 
level of approximation, is obtained as the solution of (Cahill et ai, 1987) 

r 0+ (p; P) = I -0^ ~IV(p - q) 1 ,S{q - iP)r 0+ (p; P)S(q + \P) lv , (7.23) 

S is the dressed quark propagator, obtained as the solution its DSE with a model dressed gluon prop- 
agator, and D + is the scalar-diquark propagator, which can be extracted from the (qq)o+ (<?<?)o+ 
scattering matrix. The first study of Eq. ( |7.22|) (Burden et ai, 1989) used parametrisations of the quark 
and scalar-diquark propagators based on DSE studies. These lead to a quark-core contribution to the 
baryon mass of ~ 1.2 GeV. Various studies of the effect of the pion-cloud suggest that this will lead to a 
(—200) ~ (—300) MeV contribution and hence this is a good result. However, there was a flaw in this 
study: the parametrisations of the quark and scalar-diquark propagators introduced unphysical poles 
into the integration domain, poles which would not be present in the kernel of Eq. (|7.22| ) if confinement 
is taken into account. Current studies (Cahill, 1993) are addressing this problem. 



7.2 Effective Actions and QCD 

Coupled BSE-DSE systems such as the ones described above arise naturally in the effective action 
approach to field theories. Indeed, the coupled DSE-BSE approach and the effective action approach 
are equivalent; a given effective action serves to specify a certain order of resummation of the Feynman 
diagrams that contribute to a given DSE. The application of effective actions and their interplay with 
the DSE tower has been discussed by McKay and Munczek (1989) and reviewed by Haymaker (1991). 
There are two commonly used effective actions. 



CJT Effective Action. The CJT effective action (Cornwall et ai, 1974) which can be written in 
Euclidean space as: 



T CJT [S] = TrLn S^S 



Tr 



l-S^S +V 2 [S] 



(7.24) 



where S is the dressed fermion propagator, So is the free fermion propagator and V^S 1 ] is the sum of 
all two-particle, irreducible vacuum diagrams (Stam, 1985): 



V 2 [S] = iTrf^n 



■y 2 



d A k 



d 4 q 



(27r) 4 J (2vr) 



D„ u (k - q) tr [ lfl S(q) T u (q, k) S{k)} (7.25) 



where D^ u is the dressed gluon propagator and T u is the dressed quark-gluon vertex. From Eq. (|T 
one recognises this as 



d 4 k 



tr [E(&) S{k)\ . 



(7.26) 



92 



It is therefore obvious that the variational equation 

5 



5S 



-r 



CJTt 



leads to S" 1 = S^ 1 + E 



(7.27) 



which is the quark DSE. 

One can make use of Eq. ([7.27Q to find that, evaluated at the stationary point, 



r CJT [5] = TrLn 
which was used above in Eq. ( |6.48 ). 



S^S 



+ ±Tr [l - S 



(7.28) 



Auxiliary Field Effective Action. The "auxiliary field" effective action, T AF [E], is best introduced by 
way of an example since it cannot be written in a closed form. Consider a field theory defined by the 
generating functional 

Z[J] = jfj VqVq exp [~A%q] + J d 4 xd 4 y J a p(x,y)q a (x)q p (y)^ (7.29) 

with an action 

A[q,q\ = J d A xd 4 y [gjx)^ 1 ^ - y)q a (y) - \q a {y)q p {x) M 5 ^{x,y) q^{x)q 5 {y)] , (7.30) 

where Ai s ^(x, y) = AisaiVi x )> anc ^ normalised such that Z[0] = 1. Equation ( |7.29|) can be rewritten us- 
ing the Hubbard- Stratonovich transformation to introduce hermitian, bilocal auxiliary fields, f3 a p(x,y). 
One multiplies Z by 



/ Il'^'/^«/3 e_7HS ^ _ ^ dixd4y Ja l 3 ( x ' v ^ al3( ' x,y ' } 
a/3 



(7-31) 



YlVfiap exp I- d 4 xd 4 y [/3 al3 (y, x) \ M p "{x, y) /3 jS (x, y) + J a p{x,y)(3 a p{x,y) 

a(3 ^ J 

and renormalises; interchanges the order of integration and performs a change of functional-integration 
variables: 

Pap(x, y) -> Papfa y) + qp{x)q a (y) (7.32) 
which yields a functional integral over the quark fields that is quadratic and hence can be evaluated: 



Z[J] = jn f Yl^Pafi exp -A[(3] - f d 4 xd 4 y J a p(x,y)p a p(x,y) 

a/3 L J 



where 
with 



A[p] 



-TrLn 



Ga(3 fa y) 



HS 



(7.33) 



(7.34) 



Qaji x ^v) = S 1 {x - y)S aP + E af3 (x,y) = S 1 {x - y)S a/3 - Mf 7 {x, y) (5 lS {x, y) . (7.35) 
The auxiliary field effective action, r AF [/3], is obtained from Z via a Legendre transformation 



T AF [(3] + J d 4 xd 4 yJ a p(x,y)Mx,y) = W AF [J] = - \n[Z[J\] 



(7.36) 
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with 

SW AF \J] 5T AF \P] 

; = -Pap(x,y) and — — = J a p(x, y) , (7.37) 

5J a p{x,y) 6f3 af 3{x,y) 

from which it is clear that, at tree level, 

r AF [/3] « A[/3] . (7.38) 

The stationary point equation 

(U ' J ~ yields E aP (x,y) = MfJx,y)g jS (x,y) (7.39) 



which is simply the Hartree approximation to the DSE for the fermion self energy in this theory. [As 
an example one might consider the cases of QED3 and QED4 , which are obtained with 



M%(x,y) = SirfMrisaD^x - y) , (7.40) 

and serve to illustrate this point.] It follows that in Hartree approximation r CJT and r^g e _ level have 
the same stationary point equation and solutions. Indeed, when evaluated at the stationary point, 
r^g e _ level =r CJT . However, r^g e _ level is bounded below whereas r CJT is not (Haymaker, 1991) and so, 
for many applications, r^g e _ level is more useful. 



QCD Phenomenology using the Auxiliary Field Effective Action. The generating functional in equa- 
tion ( |7.29| ) is directly applicable to QED 3 and QED 4 since the gauge boson field can be integrated out 
to yield a quartic four-fermion interaction [see, for example, Roberts and Cahill (1986)]. Here, however, 
we are interested in its application to QCD. A model field theory that has received a good deal of 
attention is defined by the action: 

A[q,q] = J d 4 xq(x)[ 1 -d + M}q(x) + ^Jd 4 xd 4 yj;(x)g 2 D, u (x-y)j a u (y) (7.41) 
where M is the quark mass matrix, j*(x) = q{x)^-j^q(x) and 

D(k 2 )+^) (7-42) 

represents a model quark-quark interaction, realistic constraints on which are discussed in Sec. [5[ [We 
note that a form of the Nambu-Jona-Lasinio model can be obtained by writing D(k 2 ) = constant but 
this does not satisfy the constraints discussed in Sec. |^.] With this current-current interaction one 
identifies 

y J h F a F (lv)s D a D f — J (7.43) 

where D, C, F refer to Dirac-, colour- and flavour-space, respectively. Although such a model field 
theory does not include explicit three- or four-gluon vertices, non-Abelian effects can be incorporated 
through the structure of D{k 2 ), as discussed in Sec. |]. Many independent studies can be correlated 
within this framework; for example, those of McKay and Munczek (1985); Cahill and Roberts (1985); 
Praschifka et al. (1987a); McKay et al. (1988); Barducci et al. (1988); Haymaker (1991) and references 
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d 4 k 
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therein. The nature of the truncations and approximations necessary to obtain Eq. ( [7.41 ) from QCD 
is discussed by Cahill and Roberts (1985) and McKay et al. (1988). 



As mentioned above, there is a large body of work in this area and we choose to summarise it by way 
of a discussion of a series of studies using a form of Eq. ( |7.41| ) that has come to be called the Global 
Colour-symmetry Model [GCM] (Roberts and Cahill, 1985) and which use the following form for the 
model dressed gluon propagator: 

Oik*) = ^fl (7.44) 



with q;i given in Eq. ( 16.791) . This form of the model has been used successfully to calculate the kinematic 
and dynamical properties of mesons and has also been applied to baryons. One ingredient of its success 
is that it incorporates confinement in the sense of ensuring the absence of quark production thresholds 
in any Feynman diagram associated with a physical process, a criterion that is discussed in Sec. B72. 
From the point of view of a connection between the GCM and QCD, Landau gauge, £ = 0, is the 
best choice, as we have seen in the DSE studies described above. However, in developing an efficacious 
phenomenology, a Feynman-like gauge with D^ v {k) = 5^ai(k 2 )/k 2 can be used without penalty. 



The generating functional for the GCM is given by Eqs. ( |7.29[ ) and (7.41) from which the auxiliary-field 
effective action, in the form of Eq. (|7.33|) , follows. In making contact with low energy phenomenology 
in QCD and the tower of DSEs it is useful to express the generating functional in terms of a local 
functional integral. This can be achieved (Roberts and Cahill, 1987) by writing 



f3 a p(x,y) = B aj3 (w = \{x + y),z = x - yj 



e 



a/3 1 



* e 



a/3 



[W,Z 



where {6^(to, z); k = 0, 1, . . .} is a complete set of orthogonal functions, to be specified later, and 
{G Q/ g(2;)} are those solutions of Eq. (|7.39p which actually minimise the action: the translationally 
invariant "vacuum field configuration". Using Eq. ( |7.45| ), the generating functional becomes 



k a/3 



-Am 



d 4 w 



E 

k 



J a/3 ^a/3 1 



w) 



(7.46) 



where A[$\ follows from direct substitution of Eq. (|7.45|) in Eq. (|7.34j ). It is clear that since the 
solutions of Eq. ( |7.39|) have been isolated explicitly the lowest order contribution from *4.[$], in an 
expansion about $ = 0, is of 0($ 2 ). The stationary point equation in the GCM is obviously the 
rainbow-ladder approximation to the quark DSE. Clearly then, of all the combinations (a/?), the only 
nonzero ® a p are the isoscalar, colour-singlet, Dirac- vector and Dirac-scalar combinations, which are 
related to the Fourier transforms of the functions A(p 2 ) and B(p 2 ), much discussed above. It follows 
that, in any expansion of A[<&] about $ = 0, the dressed quark propagator will be an important element. 



In Eq. ( |7.46| ) particular (a/9)-combinations correspond to local fields with the quantum numbers of ir- 
mesons, p-meson, etc., and whose internal structure is described by the associated Q( g-functions. This 
is most easily demonstrated by way of an example. Considering an expansion in the (a/?)-combination 
that corresponds to the isovector, colour- singlet, Dirac-pseudoscalar channel; i.e., the 7r-meson, one 
obtains 



A[n] = - TrLn 



S 1 (x -y)- g 2 D IMU 



x 



v)i. 



(f(x + y),x-y) 7, 



A" 



+ J HS [vf] (7.47) 
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where S(x — y) is the dressed quark propagator and, with i,j isospin labels, 



with 



d 4 qi d 4 q 2 



4 9 2 Dfiu(qi ~ Q2)tr 



mp^i)i^®)(-p^2)iv- 



(2tt) 4 (2tt) z 

where P is the centre-of-mass momentum, which follows from Eqs. (|7.31|) and (|7.43| ). 
At 0(7r 2 ), neglecting the constant vacuum energy, one obtains 



(7.48) 



(7.49) 



Am 



d 4 p 

(2^ 



7r*(P)I*Ag[e](P)7r<(-P) 



where 
with 



ASm(p) = *%fm(p) + ^sm(p) 



r%H l m(p) 



d A k d A q d 4 r 

(2^(2^(27)^ u ^ s,)l J ""' 71 



(7.50) 
(7.51) 

(7.52) 



xtr 



S(k + ±P) 7 ^6*(P, k - q) 7 „£ S(k - |P) 7p ^ e;.(-P, fc - r) 7(7 ^ 



From Eq. (|7.50|) it is natural to identify 7 fAfj[G](P) as the inverse of the matrix propagator for the 
7f fc -field and it is clear that, with the functions Q/ g chosen such the action is minimised, 7 ?Afj[9](P) is 
non-negative for all P 2 > 0. 

Until now the functions {Q h (w, z), k — 0, 1, . . .} have not been specified. If one chooses them such that 

(7.53) 



*Ag[e](P) =0 , k^l and 



^Ag[6](P) 
5Q k (-P,q) 







which entails that the convolution 

rt(Pp) = 

satisfies 

d A q 



^9 2 D^P-q)l^®*{P,(l)l^ 



rt(Pg) + / ^9 2 D» v (p-q)^S( q +lP)TKP,q)S(q-lP)^ =0 



(7.54) 



(7.55) 



which is the ladder BSE in the GCM, then the local field variables introduced via Eq. (|7.45|) represent 
meson fields whose internal structure is described by a ladder-BSE, the kernel of which involves the 
model dressed gluon propagator and the dressed quark propagator obtained with this gluon propagator 
via the quark-DSE. The amplitudes T k can be normalised in the standard fashion (Itzykson and Zuber, 
1980, pp. 482-485). With this interpretation of Eq. (|7.45|) then the superscript k labels excited states 
in a given channel so, in the present example: k = represents the 7r-meson; k — 1, its first excited 
state; etc. 
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This discussion demonstrates that the tree-level effective action for the GCM describes non-pointlike 
meson-fields whose internal structure is given by the solution of the coupled rainbow-DSE-ladder-BSE 
system; i.e., it is the "quark-core" of the hadron. The connection between the auxiliary-field effective 
action approach and the coupled DSE-BSE approach discussed in Sec. |7.1| is therefore clear. The higher 
order terms in the tree- level effective action, 0(7r 4 ) etc., describe couplings between these quark-core 
pions and, although only the pion has been discussed explicitly, the approach can be generalised to 
incorporate other mesons (Praschifka et ai, 1987a) and baryons (Cahill et ai, 1989). 

In order to study low energy properties of QCD it is common to employ a derivative expansion of tree- 
level effective action, Eq. ( [7.47| ). In the GCM one obtains a complete effective action, a small portion 
of the real part of which we illustrate here: 

d 4 x {| (d^TX ■ <9 M 7? + vtl\tx ■ 7?) + \ (d p p v ■ d p p v - ■ d v p v + m 2 ^ ■ p M ) 

-gpmrP» ■ t? x <9 M 7? + («, T], u, a h . . .) + tr N (7 ■ d + m N - -^m N ^f -jfj N + ... } . (7.56) 

Although only shown to second order in tt, the constraints placed by chiral symmetry on the nature 
of pion self-interactions (Roberts et ai, 1988; Roberts et ai, 1994) and on interactions between the 
pion and other hadrons (Roberts et ai, 1989) are properly represented. The full action also has an 
imaginary part which contains the Wess-Zumino term, as any realistic model of QCD must (Praschifka 
et ai, 1987a). In following this procedure one obtains the bound state masses by solving the BSE [or 
Faddeev equation] , as discussed above, and all of the couplings between quark-core mesons and baryons 
as integrals whose integrand contains the dressed quark propagators and Bethe-Salpeter amplitudes. 



For example, if one retains only the dominant functions in the Bethe-Salpeter amplitude for the tt- 



and p-mesons, as discussed in Sec. [7TT| : T p the dominant vector-amplitude and T w = B the dominant 
pseudoscalar-amplitude; then, writing the dressed quark propagator as S(p) = — ry ■ pcv(p) +&s(p), 
one has 



fl = ^ / dxxB 2 (ay-2 [a s a' s + xa v a' v ] 

87T Z JO K 



X 



II 



l \2' 



(7 c 



X 



ll 

a v a v 



I \2 



0" T 



,(7.57) 



which is obtained from the normalisation condition for the pion's Bethe-Salpeter amplitude in ladder 
approximation (Alkofer et al., 1993): 



(7.58) 



Nr. 



d 4 k 



with k, 



a/3 



(2vr) 4 

k + \ap + \f3 q, and 



IVC*; P)S(k .)T n (k; -P) dS }^° +) + T n (k; P) d T n (k; -P)S(k 0+ ) 



8P U 



8P,,. 



Nc_ 
4tt 2 



r 2 

I dxxos 
Jo 



(7.59) 



with the result that rn^f 2 



[m. 



rn d 2 ) (qq)^. [Recall that this expression for f„ differs from that 



in Eq. flb.27p and that the difference is due to neglecting contributions to the RHS of Eq. ( |6.26| ) from 
^hijp' iP)i p' 7^ P- When Z(p 2 ) = 1, both forms are the same.] One also finds (Praschifka et ai, 1987b) 



gp-n 



1 N c 
f P f 2 2vr 2 Jo 



ds s T p {s) A{s) B(s) ^4^2)2 + °K) . 



(7.60) 



with f p given by an expression similar to that in Eq. ( [7.57] ) and obtained in the same way. [The fact 
that r,r = B is a manifestation of Goldstone's theorem in this approach, as it also is in the coupled 
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DSE-BSE approach (Delbourgo and Scadron, 1979). This can be seen heuristically by substituting 
r fc oc i r fc 75 E in Eq. ( [7.55Q and noting that the equation that results is solved by E = B] Because of 
the structure of the model dressed gluon propagator, Eq. (|7.44j) , all of these integrals are finite [except for 



the quark condensate, which diverges as [In Au V /AQ CD ] dM , as it must in a model that is consistent with 
the ultraviolet constraints imposed by QCD]. The phenomenological success of the GCM is illustrated 
in Table W71 . 



One can proceed beyond the tree-level effective action and calculate the effect of meson-loops. The 
introduction of the auxiliary fields, and their subsequent identification with the quark-core; i.e., ladder 
bound states, via the connection with the BSE, Eq. ( |7.55|) , provides a mechanical framework for the 
systematic ordering and summation of diagrams. As examples we note that the model has been used 
successfully to calculate the u-p mass splitting induced by p — > tttt — > p, p — > ujtt — > p, uj — > pn — > u> 
and uj — > 717171 — > uj self energy corrections (Hollenberg et ai, 1992) and to analyse the contribution of 
pion loops to the electromagnetic charge radius of the pion (Alkofer et ai, 1993), which is found to be 
an additive correction to r\ of < 15% at m n = 0.14 GeV. Furthermore, recent calculations of the uj-p 
mixing component of the N-N potential (Goldman et ai, 1992; Krein et ai, 1993; Mitchell et ai, 1994) 
also fit neatly within this framework. The important feature of these meson-loop contributions is that 
they too are finite because the Bethe-Salpeter amplitudes that describe the quark core of the hadrons 
appear in every integral and provide natural momentum cutoffs, thus ensuring convergence. 

This example illustrates the usefulness of the auxiliary- field effective action [or, equivalently, the coupled 
rainbow-DSE-ladder-BSE] approach to QCD phenomenology. A few parameters [two or three in the 
studies described in this section] in the model dressed gluon propagator, chosen in order to parametrise 
the unknown behaviour of the quark-quark interaction in the infrared, are all that is necessary in order 
to construct a model field theory, with confinement and dynamical chiral symmetry breaking, which 
provides a good description of low energy hadronic phenomena. Furthermore, the comparison of the 
calculated results with experimental data allows one to place constraints on the nature of the quark- 
quark interaction at small-g 2 . This phenomenological success provides considerable incentive to resolve 
the questions regarding the formulation of the DSEs in QCD. 



8 Summary, Conclusions and Outlook 

We began with an introduction to the Dyson- Schwinger equation [DSE] formalism with specific reference 
to QED and to the renormalisation procedure in this framework. It was then relatively straightforward 
to extend this discussion to QCD. Following this we presented a survey of studies of strong-coupling 
QED, both because it is an interesting theory in its own right and because it is an ideal pedagogical tool 
for DSE studies and their application to QCD. We reviewed the DSE studies of the infrared behaviour 
of the gluon propagator in QCD and its role in phenomenological studies of dynamical chiral symmetry 
breaking and confinement. Finally, we collected all of this information together and illustrated how it 
can be used in relativistic bound state calculations using the BSE for mesons and the Faddeev equa- 
tion for baryons. This information also provides the foundation for the phenomenologically successful 
effective-action approach to QCD. We illustrated this approach and demonstrated its equivalence, and 
importance, to the coupled Bethe-Salpeter-Dyson-Schwinger equation studies. 

The strengths of the DSE approach to nonperturbative field theory are its manifest covariance and the 
natural way in which dynamical chiral symmetry breaking can be described. It is also straightforward 
to build in the correct perturbative limits for QCD because of asymptotic freedom. From the various 
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Table 7.2. This is an illustrative set of calculated physical quantities, 
obtained in the GCM using the propagator of Eq. (|7.44| ). Here m qq are 
effective diquark-masses and ctj are scattering lengths in 7r-7r scattering. 
The superscripts indicate the reference that a particular result is taken 
from: 1) Roberts et al. (1988); 2) Roberts et al. (1994); 3) Praschifka et 
al. (1987a); 4) Praschifka et al. (1989); 5) Hollenberg et al. (1992); 6) 
Cahill et al. (1987); 7) Cahill (1992). 





Calculated 
Massless u,d 


Experiment 

(where applicable/known) 







(0 if quarks massless) 


f 2 


0.091 GeV 


0.093 GeV 


m u 3 ' 4 


0.745 


0.783 




0.030 


0.013 


1 P 


0.232 


0.154 


3 4 

m fl 13,4 


1.310 


1.283 


3 4 

m fl - m w ^ 


0.565 


0.500 


™qq 4,6 
m o+ 


0.607 




m qq 4,6 


1.170 




qq 6 
iu _ 


0.948 




qq 6 


1.950 




m N 7 


1.20 ~ 1.30 


0.939 


r 1,2 


0.59 fm 


0.66 fm 


0° 2 


0.17 (dimensionless) 


0.20 (dimensionless) 


a 22 


-0.048 


-0.037 


a} 2 


0.030 


0.038 


a° 2 


0.0017 


0.014 


a 22 


-0.0011 





model studies it is clear that once the scale of one physical parameter; for example, f„, is fixed by 
adjusting the infrared behaviour of the quark-quark interaction, then other quantities are automatically 
of approximately the right magnitude. The fine tuning relies on the details of the models, with very 
successful descriptions of a broad range of phenomena being possible with only two or three parameters. 

While much has been accomplished, it is apparent that much more needs to be done. It is important 
to develop a better understanding of the detailed way in which confinement is manifested. Since the 
DSE formalism is not closed at any finite order it will always be important to draw on information 
gathered from other nonperturbative studies. For example, it would be very helpful if we could learn 
more about the infrared behaviour of the gluon propagator and the role of ghosts in QCD from lattice 
gauge theory studies. However, because of present day limitations on lattice studies, this information is 
currently unavailable but it is an area where cross-fertilisation could provide significant benefits. The 
results of recent lattice simulations have shown that such information may be available in the not too 
distant future. 

An important question, in field theory in general, and in the DSE approach in particular, is the connec- 
tion between Euclidean-space and physical Minkowski space. As we demonstrated, a straightforward 
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transcription of the variables in the DSE will almost never be correct. A direct formulation and solution 
of the DSEs in Minkowski space is one possibility but this is presently prevented by the complete lack 
of information about the n-point functions on the timelike axis at all but very large momenta. At the 
present time, both lattice and DSE studies are carried out in Euclidean space and the axioms of field 
theory used to extract values of physical observables. It is not yet clear when or if it will be possible to 
perform DSE studies directly in Minkowski space. 
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A Landau-Khalatnikov Transformations 

As we discussed in Sec. |3.7|, on any Ansatz for the fermion-gauge-boson vertex in Abelian gauge theories, 



one would like to impose the constraint that, when used in a DSE, it ensures local gauge covariance of 
the solutions. The gauge transformation laws which relate the propagators and fermion-gauge-boson 
vertex in QED to their Landau gauge counterparts were first given by Landau and Khalatnikov (1956) 
and Fradkin (1956). These rules are most easily specified in coordinate space and we give below the 
corresponding Euclidean space transformation laws. 

In an arbitrary gauge, the photon propagator is modified from its transverse, Landau gauge, form 
D liv {x' 1 0) by the addition of a longitudinal piece parameterised by an arbitrary function A: 

D^(x; A) = D^(x; 0) + ^A(x). (A.l) 
The corresponding rule for the fermion propagator is 

S(x; A) = S(x; oy 2 l^)-H*)] ? ( A .2) 
where e in the exponent is the gauge coupling constant. The rule for the fermion-photon vertex is 
B,(z;x,y\A) = B,(z; x,y\0)e^- A ^ 

+ S(x-y;0)e e2[Ai0) ^ x ~ y)] ^-[A(x-z)-A(z-y)}, (A.3) 

oz^ 

where is the non-amputated vertex defined in momentum space in terms of the amputated vertex, 
r M , by 

B,(p, q) = S(p)T u (p, q)S(q)D ta/ {p - q). (A.4) 
The transformation rule for the partially amputated vertex 

A fl (p,q) = S(p)T^p,q)S(q) (A.5) 
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follows from Eqs. ([A. 1 ) , (|A.2 ) and ( |A.3|) and is simply: 

A^;x,y|A) = A,(z;x,y\0)e^- A ^. 



(A.6) 



In the usual covariant gauge fixing procedure the photon propagator takes the form 

1 



p[i + n(p)] V k 2 j x /. 1 

which is obtained by taking A in Eq. ( |A.1| ) to be 

d A k e~ lk - x 



A( X ) = -e/ 



(27r) d A; 4 ' 



(A.7) 



(A.8) 



Within this set of gauges, one finds that in QED3 the transformation rule for the fermion propagator, 
Eq. (|A.2| ), becomes 

S( X ;0 = S(x;0)e- e2ilxl/8 \ (A.9) 



The free massless propagator is S 1 (p; 0) = ij ■ p which corresponds to the following function in con- 
figuration space: 



S(x;0) 



7 • x 
47r|x| 3 



(A.10) 



Applying Eq. ( |A.2| ) one obtains the LKF transformed function in an arbitrary covariant gauge: 



7 • X c -eH\x\/8n 

47r|x| 3 



(A.ll) 



For £ > one may evaluate the Fourier amplitude directly to obtain 

,2 



—27 • p 

p2 



1 - - — arctan —— 
8np \ e l i 



(A.12) 



B Quark Condensate 

We discuss here the meaning of the quark condensate (qq). It is convenient to perform a Wick rotation 
and evaluate (qq) in Euclidean space. Then using 



/ 



d 



we can write Eq. (|6.16|) as 



(qq) 



■■) = n 2 dee fit) 

Jo 



3 fA-lv d[2 £ 2 Z{£ 2 )M{£ 2 ) 



(B.l) 



(B.2) 



4tt 2 Jo e 2 + M 2 (e 2 ) 

where we have chosen /j, = Auv a s in Eq. (|6.7|) and Eq. (|6.9|) etc. Now, writing an explicit subscript 
Auv where appropriate, we have 

(qq)a 



3 r A w 



J_ 2 £ 2 Z Auv (£ 2 )M Avv (£ 2 ) ^ 

An 2 Jo ( " £ 2 + M\ (£ 2 ) a uv ^oo An 2 Jo 



d£ 2 M Avv {£ 2 ) . (B.3) 
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The second result follows since the UV dominates for large Auv and since in Landau gauge Za uv (£ 2 ) — 1, 
for large £ 2 . In the case of exact chiral symmetry [i.e., m.A uv = 0] Eq. ( |6.14| ) implies that M\ vv (£ 2 ) = 
(c/£ 2 ) \n{£ 2 1 'Afajj)*"- 1 and hence, from Eq. ([BID , we find 



\n(A 2 uv /A 2 QCD ) 



Auv^oo 47T 2 du 

Since c is independent of the choice of \x then Eq. ( |6.14| ) follows immediately. 



(B.4) 



We see from Eq. ( |6.15| ) that Ma uv (£ 2 ) =m/[| hi(l 2 / k 2 Q CD )] dm in the presence of an ECSB mass, mA uv 7^ 
0. It is clear that Eq. ( |6.17| ) will no longer result from Eq. (|B.3| ) and so the quark condensate, as defined 



in Eq. (|6.16 ), no longer varies with the renormalisation point as (qq) ~ [ln(fi 2 /Aq CD )} 



C Pion Decay Constant 



We derive the expression for f 2 given in Eq. ( |6.27| ). First note that only terms up to O(k) need be kept 
since we will be taking k — > 0. For any function of p 2 , say f(p 2 ), we have 



f((p + kf) = f( P 2 ) + 2k-p ^- 2 {p 2 ) + 0(k 



(CI) 



In Eq. (|6.27| ) the trace is over isospin, spinor, and colour labels. The spinor trace simplifies things 
considerably since [in an obvious shorthand notation] 



tr [S{p) 



k ■ dS(p)}[T^(p,p) + 0(k) l5 ][S(p)][(r n /2)Yl 5 } 
trtk ■ dS(p)T™(p,p)S(p)(T n /2)ri 5 ) + 0(k 2 ) 



(C.2) 



where we have used the fact that tr {S '(jo)75>S ^(p)7 <T 7s} = 0, which follows from Dirac matrix algebra, 
and the assumption that T 5 (p',p) = T 5 (p,p) + O^k)^, where the C(/c)7 5 term is assumed to have 
a purely 75 spinor structure. Now we can use T 5 (p,p) which is known from the Goldberger-Treiman 
relation, Eq. ( |6.21| ). We can now use the above results in Eq. ( |t).26p , operate on both sides with {k a /k 2 ) 
and use tr(r m r n ) = 25 mn to give 



3i 



f= — 



d 4 p 



k 2 J (2tt) 



^ tr S pi nor 



k-dS(p) l5 B(p 2 )S(p)A(p 2 )^ 5 



(C.3) 



where the trace is now over spinor labels only. The factor of three comes from the trace over colour. After 
evaluating (d/dp u )S(p) and taking the spinor trace we find the result of Eq. ( |6.27| ). For completeness 
we also give here the Euclidean-space form of Eq. (|6.27|) , which is 



J 7T 



4vr 2 



dp 



2 p 2 Z(p 2 )M(p 2 ) 
[p 2 + M 2 (p 2 )}< 



M{p 2 ) - 



p 2 dM 
~2df 



(C.4) 
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